Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,434 G→A 5.5% R34* (CGA→TGA)  pyrL ← pyrBI operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,4340GA 5.5% 482.2 / 11.8 146R34* (CGA→TGA) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (68/70);  new base A (3/5);  total (71/75)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

Only 100 of 146 total aligned reads displayed.

GTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGG  >  NC_000913/4472287‑4472581
                                                                                                                                                    |                                                                                                                                                   
gTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGa                                                                                                                                                    >  2:1039093/1‑149 (MQ=255)
  tGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGAtt                                                                                                                                                  <  1:986987/149‑1 (MQ=255)
  tGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGAt                                                                                                                                                   <  4:778626/148‑1 (MQ=255)
     gggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAAGCCCCTCGATTGa                                                                                                                                                >  1:1034152/1‑149 (MQ=255)
     gggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGa                                                                                                                                                <  3:677819/148‑1 (MQ=255)
     gggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAg                                                                                                                                               <  1:970494/149‑1 (MQ=255)
      ggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAgg                                                                                                                                              >  2:628388/1‑149 (MQ=255)
       gTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAggg                                                                                                                                             >  1:544857/1‑149 (MQ=255)
        tttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAgg                                                                                                                                              >  1:171354/1‑147 (MQ=255)
        tttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAgg                                                                                                                                              <  2:171354/147‑1 (MQ=255)
        tttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATGGAggg                                                                                                                                             >  3:719288/1‑148 (MQ=255)
         ttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCAATCGAGGGGc                                                                                                                                           <  2:292988/149‑1 (MQ=255)
           gCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCt                                                                                                                                          <  3:507254/148‑1 (MQ=255)
           gCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCAATCGAGGGGCtt                                                                                                                                         <  2:479122/149‑2 (MQ=255)
            cTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTg                                                                                                                                         <  4:111264/148‑1 (MQ=255)
             tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCt                                                                                                                                          >  2:941504/1‑146 (MQ=255)
              ttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCgggg                                                                                                                                       >  3:736442/1‑148 (MQ=255)
               tCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGa                                                                                                                                      <  4:479250/148‑1 (MQ=255)
               tCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCAATCGAGGGGCttttt                                                                                                                                      <  4:716164/148‑5 (MQ=255)
                 aGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATg                                                                                                                                   >  1:224004/1‑149 (MQ=255)
                  gTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGGCGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                                                                                  <  2:605124/149‑1 (MQ=255)
                  gTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                                                                                  >  1:1077526/1‑149 (MQ=255)
                   tttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGGCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgt                                                                                                                                  >  3:1037217/1‑147 (MQ=255)
                   tttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAgggg                                                                                                                                            >  3:1031889/1‑138 (MQ=255)
                   tttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAgggg                                                                                                                                            <  4:1031889/138‑1 (MQ=255)
                    ttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGaa                                                                                                                                     <  4:31110/144‑1 (MQ=255)
                    ttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGaa                                                                                                                                     >  3:31110/1‑144 (MQ=255)
                       gccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAt                                                                                                                             <  1:576050/149‑1 (MQ=255)
                        ctgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAAAGATATGTTTCTGATAGAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAt                                                                                                                             >  4:658310/3‑148 (MQ=255)
                        ccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAt                                                                                                                             <  3:939844/148‑1 (MQ=255)
                                gCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                                                                                    <  2:377254/149‑1 (MQ=255)
                                 ccAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaag                                                                                                                   >  2:817655/1‑149 (MQ=255)
                                      aCCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaa                                                                                                               <  3:749906/148‑1 (MQ=255)
                                       ccAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaaa                                                                                                              <  3:299446/148‑1 (MQ=255)
                                           aTTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGCATGGGTGATCAACGGGAAGAAAAACGGCa                                                                                                         <  2:288066/149‑1 (MQ=255)
                                                ggTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCa                                                                                                    >  2:561542/1‑149 (MQ=255)
                                                   cATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGc                                                                                                  <  2:1034152/149‑1 (MQ=255)
                                                     tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTc                                                                                               >  2:349291/1‑149 (MQ=255)
                                                      cgcgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTc                                                                                               >  4:361494/1‑148 (MQ=255)
                                                           cTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGACCAGCGTCttttt                                                                                          <  3:76528/148‑1 (MQ=255)
                                                           cTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑‑AAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                                                                        <  1:48799/149‑1 (MQ=255)
                                                            tAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttttt                                                                                         >  4:638770/1‑148 (MQ=255)
                                                            tAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                                                                        <  1:243249/149‑1 (MQ=255)
                                                                gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGa                                                                                     <  3:577470/148‑1 (MQ=255)
                                                                 tCGTTTATGGAAATGATATGTTTCTGATATAGCGTATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAc                                                                                    >  4:67309/1‑148 (MQ=255)
                                                                    tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgcg                                                                                 >  3:559904/1‑148 (MQ=255)
                                                                    tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCgg                                                                                <  2:132734/149‑1 (MQ=255)
                                                                      tATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTa                                                                              >  1:893164/1‑149 (MQ=255)
                                                                       aTGGAAATGATATGTTTCTGATATAGCGTATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTa                                                                              <  3:67309/148‑1 (MQ=255)
                                                                       aTGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTa                                                                              <  4:313354/148‑1 (MQ=255)
                                                                         ggAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGa                                                                           >  1:147952/1‑149 (MQ=255)
                                                                          gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                                                                               <  2:1052407/144‑1 (MQ=255)
                                                                          gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                                                                               >  1:1052407/1‑144 (MQ=255)
                                                                             aTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCGCGTGACGGCTGGGC‑AAGAAAAAGCCCCTCGATTGAGGGGCTGGGACTGGGAGATCAACGGGAAGAAAAACGGTAGGGCAGCGTCTTTTTTCAGACGCGGTAAGACaa                                                                        >  4:652787/1‑148 (MQ=255)
                                                                              tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaaa                                                                      <  1:1039093/149‑1 (MQ=255)
                                                                              tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaa                                                                       <  3:1787/148‑1 (MQ=255)
                                                                              tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaa                                                                       <  3:330033/148‑1 (MQ=255)
                                                                               gATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaaa                                                                      >  4:207658/1‑148 (MQ=255)
                                                                                atatGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAGAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAATACAAAAta                                                                    >  2:878565/1‑148 (MQ=255)
                                                                                  atGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                                                                   >  3:129834/1‑148 (MQ=255)
                                                                                  atGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                                                                   <  3:318954/148‑1 (MQ=255)
                                                                                  atGTTTCTGATATAGCGGATTAGCCATCGTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                                                                  <  1:789821/149‑1 (MQ=255)
                                                                                     tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                                                                  >  1:35947/1‑146 (MQ=255)
                                                                                     tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                                                                  <  2:35947/146‑1 (MQ=255)
                                                                                     tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                                                               <  1:614446/149‑1 (MQ=255)
                                                                                      ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgcg                                                                                 <  3:587602/130‑1 (MQ=255)
                                                                                      ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAg                                                                            <  3:259165/135‑1 (MQ=255)
                                                                                       tCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                                                             >  1:170699/1‑149 (MQ=255)
                                                                                           atatAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                                                         >  2:535826/1‑149 (MQ=255)
                                                                                           atatAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                                                         >  1:210306/1‑149 (MQ=255)
                                                                                            tataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca                                                           <  2:390873/146‑1 (MQ=255)
                                                                                              taGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTTATCAACGGGAAGAAAAACGTCAGGCCAGCGACTTCTCTCAGACGCGGTAATCCAAAATGTCCAACACACTGa                                                       <  4:760087/148‑1 (MQ=255)
                                                                                               aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGaa                                                      <  4:716802/148‑1 (MQ=255)
                                                                                               aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAAc                                                     >  1:550132/1‑149 (MQ=255)
                                                                                                 cGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCa                                                   <  1:535826/149‑1 (MQ=255)
                                                                                                  ggATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCa                                                   <  3:361494/148‑1 (MQ=255)
                                                                                                         ccATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcct                                           >  1:965387/1‑149 (MQ=255)
                                                                                                            tCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCt                                                                                              >  1:294200/1‑95 (MQ=255)
                                                                                                             cTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGga                                                                                                                     <  3:794914/71‑1 (MQ=255)
                                                                                                             cTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTccc                                        <  4:590207/148‑1 (MQ=255)
                                                                                                                   tctcCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCaaa                                  <  4:437677/148‑1 (MQ=255)
                                                                                                                    ctcCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCgg                                                                                >  1:302053/1‑101 (MQ=255)
                                                                                                                    ctcCTGACGACTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCgg                                                                                <  2:302053/101‑1 (MQ=255)
                                                                                                                     tcCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTg                               >  1:214304/1‑149 (MQ=255)
                                                                                                                     tcCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTg                               <  2:704534/149‑1 (MQ=255)
                                                                                                                       cTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGt                              <  3:963470/148‑1 (MQ=255)
                                                                                                                          aCGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCg                           >  3:617593/1‑148 (MQ=255)
                                                                                                                                gggC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATta                     <  4:559904/148‑1 (MQ=255)
                                                                                                                                 ggC‑AAAAAAAAGCCCCTCGATGGAGGGGCCGGGAATGGGTGCTCAACGGGAAGAAAAACGGCCGGCCAGCCTCTTTTTTCAGACGCGGTAAGACAAAATGTCGACCACACTGAACCAGACATCCTCCCGGCAAATTGTCCGGCATtat                    <  3:652787/148‑1 (MQ=255)
                                                                                                                                   c‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATAc                                                >  2:178165/1‑118 (MQ=255)
                                                                                                                                    aaaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACt                 >  1:858404/2‑149 (MQ=255)
                                                                                                                                      aaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTAAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCa               >  2:891403/1‑149 (MQ=255)
                                                                                                                                          ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGt           >  2:914481/3‑149 (MQ=255)
                                                                                                                                          ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGt           >  2:978268/3‑149 (MQ=255)
                                                                                                                                          aaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGt           <  2:296788/149‑1 (MQ=255)
                                                                                                                                             gCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCa         >  4:195058/1‑148 (MQ=255)
                                                                                                                                             gCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCa         >  4:482537/1‑148 (MQ=255)
                                                                                                                                               cccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGACCGGCATTATACTCATCGTCAGaa      <  1:913849/149‑1 (MQ=255)
                                                                                                                                                 cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAg     >  4:685956/1‑148 (MQ=255)
                                                                                                                                                  tCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc    >  3:269896/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGG  >  NC_000913/4472287‑4472581

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑