Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,002,402 G→A 7.3% intergenic (+26/+17) corA → / ← yigF Ni(2(+))/Co(2(+))/Mg(2(+)) transporter/DUF2628 domain‑containing protein YigF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,002,4020GA 7.3% 551.8 / 21.0 178intergenic (+26/+17)corA/yigFNi(2(+))/Co(2(+))/Mg(2(+)) transporter/DUF2628 domain‑containing protein YigF
Reads supporting (aligned to +/- strand):  ref base G (76/89);  new base A (6/7);  total (82/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

Only 100 of 185 total aligned reads displayed.

TGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAG  >  NC_000913/4002253‑4002550
                                                                                                                                                     |                                                                                                                                                    
tGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTgg                                                                                                                                                      >  2:319994/1‑150 (MQ=255)
 ggcaggcaCTGTGAGTTGATGCCAGAACTGAAGTGGCGCTTCGGCTACCCTGGCGCGATTATCTTTATGAGCCTCGCGGGCCGGGCACCGGAGCTGTACGTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGGGGCTTAGTCTGGc                                                                                                                                                     <  1:1810536/150‑1 (MQ=255)
 ggCATGAACTTTGAGTTTATGCCATAACTGAAGTGGAGCTTCGGCTACGCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGc                                                                                                                                                     <  1:1048017/150‑1 (MQ=255)
 ggCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCTGCTACCCTGGCGAGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGc                                                                                                                                                     <  1:1819303/150‑1 (MQ=255)
 ggCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTggg                                                                                                                                                     >  1:2433783/1‑149 (MQ=255)
  gCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCt                                                                                                                                                    <  1:45364/150‑1 (MQ=255)
   cATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTa                                                                                                                                                   >  2:2165352/1‑150 (MQ=255)
    aTGAACTTTGAGTTTATGCCAGAACTGAAGAGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTaa                                                                                                                                                  >  1:754051/1‑150 (MQ=255)
         cTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCAc                                                                                                                                             >  2:905857/1‑150 (MQ=255)
         cTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGGGGCTTAGCCAGACTAAGCCAc                                                                                                                                             <  2:1016223/150‑1 (MQ=255)
          tttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTTTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGATCGGTGGCTTAGTCTGGCTAAGCCAcc                                                                                                                                            <  2:1835929/150‑1 (MQ=255)
          tttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCAcc                                                                                                                                            >  2:1682857/1‑150 (MQ=255)
           ttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAccg                                                                                                                                           <  1:351492/150‑2 (MQ=255)
           ttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAccg                                                                                                                                           <  1:342450/150‑2 (MQ=255)
           ttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAccg                                                                                                                                           <  2:1590632/150‑2 (MQ=255)
             gAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGt                                                                                                                                         >  1:1096685/1‑150 (MQ=255)
             gAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGt                                                                                                                                         <  2:855844/150‑1 (MQ=255)
              aGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGtt                                                                                                                                        >  1:423890/1‑150 (MQ=255)
               gTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTa                                                                                                                                       >  2:1560058/1‑150 (MQ=255)
               gTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATACTAGCGGGCCTGGCCCCGTATCTGTACTTTAAGCCGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTa                                                                                                                                       >  2:1331102/1‑150 (MQ=255)
                   aTGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCa                                                                                                                                   <  1:2271736/150‑1 (MQ=255)
                    tgacAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCaa                                                                                                                                  >  2:899283/4‑150 (MQ=255)
                    tgacAGAACTGAAGTGGAGCTTCGGCTACACTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCaa                                                                                                                                  >  2:890302/4‑150 (MQ=255)
                    tGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCaa                                                                                                                                  >  2:1852263/1‑150 (MQ=255)
                     gCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCaaa                                                                                                                                 >  1:1884671/1‑150 (MQ=255)
                     gCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGACACCTGTTATTCaaa                                                                                                                                 >  2:2190269/1‑150 (MQ=255)
                      ccAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAg                                                                                                                                <  2:279423/150‑1 (MQ=255)
                      ccAGAACTGAAGTGGAGCTTCGGCTAACCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTCGCACCGTATCTCTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAg                                                                                                                                >  2:882084/1‑150 (MQ=255)
                           aCTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTcc                                                                                                                           >  2:1180764/1‑150 (MQ=255)
                            cTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCa                                                                                                                          <  2:2352507/150‑1 (MQ=255)
                             tGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAg                                                                                                                         >  1:1341099/1‑150 (MQ=255)
                               aaGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGt                                                                                                                       <  2:1047322/150‑1 (MQ=255)
                                  tGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTAtt                                                                                                                    <  1:2302390/150‑1 (MQ=255)
                                  tGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTAtt                                                                                                                    <  2:1300427/150‑1 (MQ=255)
                                      gCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAAc                                                                                                                >  2:1350288/1‑150 (MQ=255)
                                        ttCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAAccc                                                                                                              <  2:1889853/150‑1 (MQ=255)
                                        ttCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAAccc                                                                                                              >  1:365025/1‑150 (MQ=255)
                                         tCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCt                                                                                                             >  1:820409/1‑150 (MQ=255)
                                           ggCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCttt                                                                                                           >  1:2094354/1‑150 (MQ=255)
                                               aCCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTAcc                                                                                                       <  1:723935/150‑1 (MQ=255)
                                               aCCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTAcc                                                                                                       <  1:1475552/150‑1 (MQ=255)
                                               aCCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTAcc                                                                                                       >  2:415081/1‑150 (MQ=255)
                                                 ccTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACctc                                                                                                     >  1:1572162/1‑150 (MQ=255)
                                                 ccTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACctc                                                                                                     <  2:456286/150‑1 (MQ=255)
                                                      cgcgATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTtct                                                                                                >  2:2020977/1‑150 (MQ=255)
                                                        cgATTATCTTTATGATCCTCGCGGGCCTGTCACCGTATCTGTACTTTAAGAGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCATGTATTTAACCCTTTTACCTCTTTCTCa                                                                                              >  2:1044760/1‑150 (MQ=255)
                                                          tttATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGATAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATa                                                                                            <  2:536269/149‑1 (MQ=255)
                                                          aTTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATa                                                                                            <  1:1724812/150‑1 (MQ=255)
                                                               cTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAAcc                                                                                       <  2:2357819/150‑1 (MQ=255)
                                                                 ttATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCAt                                                                                     >  1:570058/1‑150 (MQ=255)
                                                                      aTCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGtt                                                                                <  1:652694/150‑1 (MQ=255)
                                                                        ccTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCg                                                                              <  2:1970875/150‑1 (MQ=255)
                                                                            gcgGGCCTGGCACCTTATCTGTTCTTTAAGCGGAAGAACTTGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGTGTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGtt                                                                          <  1:874484/150‑1 (MQ=255)
                                                                               ggCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAACCGCGAGCGGTGGCTTCGTCTGGCTACGCCCCCGGTTATGCAACGGGGCAAGGTATTTAACCCTTTTCCAGATTTCGCATAGACCCATTTGTTCGTTTTAAc                                                                       >  1:126871/1‑150 (MQ=255)
                                                                               ggCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAAc                                                                       <  2:1227626/150‑1 (MQ=255)
                                                                               ggCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAAc                                                                       <  1:2221621/150‑1 (MQ=255)
                                                                                gCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACa                                                                      <  1:880350/150‑1 (MQ=255)
                                                                                 ccTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAg                                                                     >  2:1091444/1‑150 (MQ=255)
                                                                                   tGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCa                                                                   >  2:1868156/1‑150 (MQ=255)
                                                                                     gCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAAc                                                                 <  1:795655/150‑1 (MQ=255)
                                                                                      cagcGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTGATAGAACCATTTGTTCGTGTTAACAGCAACa                                                                <  1:2010044/147‑1 (MQ=255)
                                                                                           tATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgc                                                           <  1:983385/150‑1 (MQ=255)
                                                                                              cTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctgc                                                        >  2:2420074/1‑150 (MQ=255)
                                                                                               tGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctgct                                                       >  1:1106390/1‑150 (MQ=255)
                                                                                               tGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGgctgct                                                       >  2:178803/1‑150 (MQ=255)
                                                                                                    tttAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGc                                                  <  2:7233/150‑1 (MQ=255)
                                                                                                     ttAAGCGGAAGAAATGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCa                                                 <  2:1082738/150‑1 (MQ=255)
                                                                                                        aGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAAtt                                              >  2:1959099/1‑150 (MQ=255)
                                                                                                        aGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAAtt                                              >  1:165145/1‑150 (MQ=255)
                                                                                                              agaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGa                                        <  1:1852263/150‑1 (MQ=255)
                                                                                                              agaaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGa                                        <  2:2305020/150‑1 (MQ=255)
                                                                                                                aaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAAt                                      >  1:2328424/1‑150 (MQ=255)
                                                                                                                aaCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAAt                                      <  1:890302/150‑1 (MQ=255)
                                                                                                                 aCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATg                                     >  1:256994/1‑150 (MQ=255)
                                                                                                                 aCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATg                                     <  1:27904/150‑1 (MQ=255)
                                                                                                                 aCTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATg                                     <  2:1221223/150‑1 (MQ=255)
                                                                                                                   tGGTTGTAAAAAACGAGAGCGTTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTGCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGtt                                   <  2:754051/150‑1 (MQ=255)
                                                                                                                    ggTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTTCCAGAATGTTa                                  >  1:1907422/1‑150 (MQ=255)
                                                                                                                      ttGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAAc                                <  1:1180764/150‑1 (MQ=255)
                                                                                                                           aaaaaCGCGAGCGGTGGCTTAGTCTGGCTAGGCCACCTGGGATGCAAAGGCTCCAGGTATGTAACCCGTTGAGGTCTTTCTCATAGAGCCATTTGTTCGTGTGAACCGCCAGATCGGCTGCTAGGGCAATTCACAGAATGTGAACGCCaa                           >  2:1138440/1‑150 (MQ=255)
                                                                                                                           aaaaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTTCGGCAATTCCCAGAATGTTAACGCCaa                           >  1:1188139/1‑150 (MQ=255)
                                                                                                                            aaaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAAt                          >  1:7479/1‑150 (MQ=255)
                                                                                                                            aaaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCAATTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAAt                          >  2:2223804/1‑150 (MQ=255)
                                                                                                                             aaaCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAAtt                         <  2:1875541/150‑1 (MQ=255)
                                                                                                                               aCGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAg                       <  1:2190269/150‑1 (MQ=255)
                                                                                                                                  agagCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGc                    <  1:953160/150‑1 (MQ=255)
                                                                                                                                  agagCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGc                    <  1:894702/150‑1 (MQ=255)
                                                                                                                                      cGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACca                >  2:586793/1‑150 (MQ=255)
                                                                                                                                      aggTGGCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACca                >  1:1373578/2‑150 (MQ=255)
                                                                                                                                        gTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaac              >  1:2025205/1‑150 (MQ=255)
                                                                                                                                        gTGGCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaac              >  2:325667/1‑150 (MQ=255)
                                                                                                                                        gTGGCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaac              >  2:67335/1‑150 (MQ=255)
                                                                                                                                         tGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaaca             <  1:1878060/150‑1 (MQ=255)
                                                                                                                                          ggCTTAGCCAGACTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACcaacaa            >  2:2459458/1‑150 (MQ=255)
                                                                                                                                             ttAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATa         <  2:1941713/150‑1 (MQ=255)
                                                                                                                                                  cTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCCCACAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGAt    <  1:319994/150‑1 (MQ=255)
                                                                                                                                                    ggCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAg  <  1:1331102/150‑1 (MQ=255)
                                                                                                                                                    ggCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAg  >  1:195408/1‑150 (MQ=255)
                                                                                                                                                    ggCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAg  <  2:2387069/150‑1 (MQ=255)
                                                                                                                                                    ggCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAg  <  2:1240482/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTCAAAGGCTCCAGGTATTTAACCCTTTTACCTCTTTCTCATAGAACCATTTGTTCGTGTTAACAGCAACATAGGCTGCTACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAG  >  NC_000913/4002253‑4002550

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: