Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,521,276 A→G 5.3% M597T (ATG→ACG)  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,2760AG 5.3% 276.2 / 4.7 94M597T (ATG→ACG) yncDputative iron outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (39/49);  new base G (3/2);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

TGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGA  >  NC_000913/1521127‑1521425
                                                                                                                                                     |                                                                                                                                                     
tGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCa                                                                                                                                                       >  1:1697908/1‑150 (MQ=255)
tGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCa                                                                                                                                                       <  1:1601953/150‑1 (MQ=255)
  tATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCata                                                                                                                                                     >  2:1701034/1‑150 (MQ=255)
    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                                                                                                                                             >  2:1616907/1‑150 (MQ=255)
          cccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGc                                                                                                                                             >  1:38376/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                                                                                                                                            >  2:415668/1‑150 (MQ=255)
           ccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCg                                                                                                                                            >  1:1575029/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                                                                                                                               <  2:119972/150‑1 (MQ=255)
                        cgacgaGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATcc                                                                                                                               >  2:190167/1‑150 (MQ=255)
                        cgacgaGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATcc                                                                                                                               >  1:1409034/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                                                                                                   <  1:1701034/150‑1 (MQ=255)
                                                    cGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCa                                                                                                   <  2:1149718/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                                                                       <  2:943040/150‑1 (MQ=255)
                                                                                cccATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGc                                                                       >  1:734889/1‑150 (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                                                                   <  1:850321/150‑1 (MQ=255)
                                                                                    tataACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAgg                                                                   <  2:1697908/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:1343065/150‑1 (MQ=255)
                                                                                                                   aTCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGc                                    <  1:745437/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:1165132/141‑1 (MQ=255)
                                                                                                                                              atcgCCCGTAATATGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAg               <  2:260720/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)
                                                                                                                                                     |                                                                                                                                                     
TGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGA  >  NC_000913/1521127‑1521425

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: