Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,521,279 T→A 5.4% N596I (AAT→ATT)  yncD ← putative iron outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,521,2790TA 5.4% 284.4 / 3.7 93N596I (AAT→ATT) yncDputative iron outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (40/48);  new base A (3/2);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

TTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACG  >  NC_000913/1521131‑1521427
                                                                                                                                                    |                                                                                                                                                    
tttATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATAtt                                                                                                                                                     <  1:558498/150‑1 (MQ=255)
  tataCCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTAc                                                                                                                                                   >  1:1373893/1‑150 (MQ=255)
    taCCCGGAGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGACAAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTACTGTACATAACCTATACACGCAAAGACCATATAACgtg                                                                                                                                                >  1:1300192/1‑148 (MQ=255)
     aCCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACggg                                                                                                                                                <  1:1483448/150‑1 (MQ=255)
      cccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGc                                                                                                                                               >  1:38376/1‑150 (MQ=255)
      cccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGc                                                                                                                                               >  2:1616907/1‑150 (MQ=255)
      cccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGc                                                                                                                                               >  2:745437/1‑150 (MQ=255)
       ccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCACAGCCCATATTACGGGCg                                                                                                                                              >  1:452976/1‑150 (MQ=255)
       ccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCg                                                                                                                                              >  1:943040/1‑150 (MQ=255)
       ccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCg                                                                                                                                              >  1:1575029/1‑150 (MQ=255)
       ccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCg                                                                                                                                              >  2:415668/1‑150 (MQ=255)
           tGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATcc                                                                                                                                          <  2:38376/150‑1 (MQ=255)
             aaTAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCt                                                                                                                                        >  1:614931/1‑150 (MQ=255)
              aTAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTg                                                                                                                                       >  1:1329628/1‑150 (MQ=255)
              aTAACCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTg                                                                                                                                       <  1:934910/150‑1 (MQ=255)
                    cgacgaGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATcc                                                                                                                                 >  1:1409034/1‑150 (MQ=255)
                    cgacgaGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATcc                                                                                                                                 <  2:119972/150‑1 (MQ=255)
                    cgacgaGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATcc                                                                                                                                 >  2:190167/1‑150 (MQ=255)
                       cgaGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGAt                                                                                                                              >  1:1149718/1‑150 (MQ=255)
                          gagTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTAc                                                                                                                           <  1:1718616/150‑1 (MQ=255)
                            gTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCg                                                                                                                         <  1:1251339/150‑1 (MQ=255)
                                    cGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCAAGATTACGGGCGATCCCTGGCATCCGCTTACCGTTACAATc                                                                                                                 >  2:1341865/1‑150 (MQ=255)
                                     ggCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATcc                                                                                                                >  2:1343065/1‑150 (MQ=255)
                                     ggCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATcc                                                                                                                <  1:1259980/150‑1 (MQ=255)
                                       cgcTTTTGGCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGAACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTg                                                                                                              <  1:186323/150‑1 (MQ=255)
                                       cgcTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTg                                                                                                              <  1:124861/150‑1 (MQ=255)
                                       cgcTTTTGCCGTATTTTCATAATCTGCCATAATATCGCACATATAACTCGCTTCCGTGCCTGCATACCAACCATATTCCGGTACATAACCTATCGACGCAAAGCGCATATTACGGGCGATCCCTGGCATCCTATTACCGTTACAATCCtt                                                                                                              <  1:363705/150‑2 (MQ=255)
                                          ttttGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTTCATACCAACCATCTTCCGGAACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTc                                                                                                           <  1:218256/150‑1 (MQ=255)
                                           tttGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCa                                                                                                          >  1:1647686/1‑150 (MQ=255)
                                             tGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACCATCCTGTTCAtt                                                                                                        >  2:694654/1‑150 (MQ=255)
                                               ccGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGc                                                                                                      >  1:492715/1‑150 (MQ=255)
                                                cGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCa                                                                                                     <  2:1149718/150‑1 (MQ=255)
                                                cGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCa                                                                                                     <  1:1701034/150‑1 (MQ=255)
                                                   aTTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAAc                                                                                                  <  1:1047020/150‑1 (MQ=255)
                                                    ttttCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACa                                                                                                 <  1:921944/150‑1 (MQ=255)
                                                      ttCATCATCTGCCATAATATCGCCCATCTAACGCGCTTCCGGGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACAtt                                                                                               <  2:833790/150‑1 (MQ=255)
                                                       tcatcaTCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATtg                                                                                              <  2:1373893/150‑1 (MQ=255)
                                                        catcatCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATtgc                                                                                             <  1:889159/150‑1 (MQ=255)
                                                            atcTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCg                                                                                         <  1:1083183/150‑1 (MQ=255)
                                                                  cATAATATCGCCCATATCACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGt                                                                                   <  2:1441469/150‑1 (MQ=255)
                                                                         tCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCa                                                                            <  2:1239108/150‑1 (MQ=255)
                                                                          cGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATTTTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCaa                                                                           >  2:307285/1‑150 (MQ=255)
                                                                          cGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCaa                                                                           >  1:568178/1‑150 (MQ=255)
                                                                          cGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCaa                                                                           >  1:55072/1‑150 (MQ=255)
                                                                           gCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATCACCTATCGACGCAACGCCCCTATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAg                                                                          >  1:402951/1‑150 (MQ=255)
                                                                           gCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAg                                                                          <  2:1028515/150‑1 (MQ=255)
                                                                            cccATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGc                                                                         >  1:734889/1‑150 (MQ=255)
                                                                            cccATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGc                                                                         <  2:943040/150‑1 (MQ=255)
                                                                            cccATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTACAAACATTGCTGCGATAGGTCGCATCAAGc                                                                         <  2:1575029/150‑1 (MQ=255)
                                                                             ccATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGcc                                                                        <  1:1616907/150‑1 (MQ=255)
                                                                             ccATATAACGCGCGTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTCTCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGGATCAAGcc                                                                        <  1:180308/150‑1 (MQ=255)
                                                                               atatAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAg                                                                      <  2:614931/150‑1 (MQ=255)
                                                                               aaataACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAg                                                                      <  1:905406/148‑1 (MQ=255)
                                                                                tataACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAgg                                                                     <  2:1697908/150‑1 (MQ=255)
                                                                                tataACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAgg                                                                     <  1:850321/150‑1 (MQ=255)
                                                                                   aaCGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTcc                                                                  >  1:606203/1‑150 (MQ=255)
                                                                                      gcgcTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAgg                                                               <  2:707058/150‑1 (MQ=255)
                                                                                       cgcTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGAATACCGTTACAATCCTGTTCATTGCAACCATTGCTGCGATAGCTCGCATCCAGCCAGGTCCAGGa                                                              >  1:1144806/1‑150 (MQ=255)
                                                                                           tCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGAcgcg                                                          <  2:366242/150‑1 (MQ=255)
                                                                                             cGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGtt                                                        >  2:1401984/1‑150 (MQ=255)
                                                                                             cGTGCCTGCAGCCCAGCCATCTTCCGGTACATCACCTATCGACGCCAAGCCCATATTACGGGCGATCCCGGGCATCCGCTTACCGTTACAGTACTGTACATTGCAAACATTGCTGCGATAGGTCGCAACCAGCCAGGTCCAGGACGCGtt                                                        >  1:47944/1‑150 (MQ=255)
                                                                                               tGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTa                                                      <  2:452976/150‑1 (MQ=255)
                                                                                                   tGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTAcgcc                                                  >  2:1700840/1‑149 (MQ=255)
                                                                                                     cATACCAACCATCTTCCGTTACATAACCTATCGACGCAAAGCCCATATTAAGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGaa                                                <  2:1345741/150‑1 (MQ=255)
                                                                                                     cATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGATTACGCGaa                                                >  2:1349852/1‑150 (MQ=255)
                                                                                                         ccaaccaTCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCATACATTGCTGGGATAGGTCGCATCAAGCGAGGTCCAGGACGCGTTTACGCGATAAtc                                            >  1:85751/1‑150 (MQ=255)
                                                                                                          caaccaTCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCa                                                                       >  2:1337001/1‑122 (MQ=255)
                                                                                                          caaccaTCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCa                                                                       <  1:1337001/122‑1 (MQ=255)
                                                                                                               aTCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGc                                      <  1:745437/150‑1 (MQ=255)
                                                                                                               aTCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGc                                      <  1:1343065/150‑1 (MQ=255)
                                                                                                                    ccGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGGTTACGCGAAAATCTCCTGCGAAAc                                 >  1:123672/1‑150 (MQ=255)
                                                                                                                      ggTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGt                               >  1:1621773/1‑150 (MQ=255)
                                                                                                                      ggTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGt                               <  2:55072/150‑1 (MQ=255)
                                                                                                                       gTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGtt                              <  2:1695458/150‑1 (MQ=255)
                                                                                                                       gTACATAACCTATCGACGCAAAGCACATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGtt                              <  2:1300192/150‑1 (MQ=255)
                                                                                                                          cATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGCAAATCTCCTGCGAAACGTTGAt                           >  2:1389651/1‑150 (MQ=255)
                                                                                                                          cATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGAt                           >  2:964442/1‑150 (MQ=255)
                                                                                                                           aTAACCTATCGACGCTAAGACCATATTCCGGGAGATCCCTTGCAGCCGCTTACCGTTACCATCCTGTTCATTGCAAACTTTGCTGAGATAGGTCGCATCAAGCCAGGTCCAGGACTCGTTTACGCGAAAATCTCCTGCGAAACGTTGATc                          <  2:1065918/150‑1 (MQ=255)
                                                                                                                           aTAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATc                          <  1:964442/150‑1 (MQ=255)
                                                                                                                             aaCCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATccc                        <  1:1619038/150‑1 (MQ=255)
                                                                                                                                cTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATg                     <  2:1329628/150‑1 (MQ=255)
                                                                                                                                 tATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGc                    <  2:237029/150‑1 (MQ=255)
                                                                                                                                  aTCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGcc                   <  2:1061351/150‑1 (MQ=255)
                                                                                                                                     gACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTAGATCCCATGCCAGt                >  2:757118/1‑150 (MQ=255)
                                                                                                                                       cGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTc              <  1:190167/150‑1 (MQ=255)
                                                                                                                                       cGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTc              >  1:407561/1‑150 (MQ=255)
                                                                                                                                        gCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCa             >  2:1398361/1‑150 (MQ=255)
                                                                                                                                          atcgCCCGTAATATGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAg                 >  1:1165132/4‑144 (MQ=255)
                                                                                                                                          atcgCCCGTAATATGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAg                 >  1:260720/4‑144 (MQ=255)
                                                                                                                                          atcgCCCGTAATATGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAg                 <  2:260720/141‑1 (MQ=255)
                                                                                                                                          atcgCCCGTAATATGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAg                 <  2:1165132/141‑1 (MQ=255)
                                                                                                                                          atcgCCCGTAATATGGGCGATCCCTGGCAGCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAATATCTCCGGCGAAACGTTGATCCCCTGCCAGTTCCgc           >  1:89277/4‑150 (MQ=255)
                                                                                                                                            aGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAgcgc         >  1:1336810/1‑150 (MQ=255)
                                                                                                                                               ccATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCtt      <  2:1409034/150‑1 (MQ=255)
                                                                                                                                                 atatTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTga    <  1:415668/150‑1 (MQ=255)
                                                                                                                                                   atTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTgacg  >  2:989556/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACG  >  NC_000913/1521131‑1521427

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