Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,521,282 C→T 5.2% R595H (CGT→CAT)  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,2820CT 5.2% 306.6 / 3.7 96R595H (CGT→CAT) yncDputative iron outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (42/48);  new base T (3/2);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

TATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACG  >  NC_000913/1521133‑1521430
                                                                                                                                                     |                                                                                                                                                    
tataCCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTAc                                                                                                                                                      >  1:1373893/1‑150 (MQ=255)
  taCCCGGAGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGACAAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTACTGTACATAACCTATACACGCAAAGACCATATAACgtg                                                                                                                                                   >  1:1300192/1‑148 (MQ=255)
   aCCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACggg                                                                                                                                                   <  1:1483448/150‑1 (MQ=255)
    cccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGc                                                                                                                                                  >  2:745437/1‑150 (MQ=255)
    cccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGc                                                                                                                                                  >  1:38376/1‑150 (MQ=255)
    cccGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGc                                                                                                                                                  >  2:1616907/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                                                                                                                                    >  2:190167/1‑150 (MQ=255)
                  cgacgaGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATcc                                                                                                                                    <  2:119972/150‑1 (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                                                                                                        <  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:55072/1‑150 (MQ=255)
                                                                        cGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCaa                                                                              >  1:568178/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:1343065/150‑1 (MQ=255)
                                                                                                             aTCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGc                                         <  1:745437/150‑1 (MQ=255)
                                                                                                                  ccGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGGTTACGCGAAAATCTCCTGCGAAAc                                    >  1:123672/1‑150 (MQ=255)
                                                                                                                    ggTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGt                                  <  2:55072/150‑1 (MQ=255)
                                                                                                                    ggTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGt                                  >  1:1621773/1‑150 (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                    <  2:260720/141‑1 (MQ=255)
                                                                                                                                        atcgCCCGTAATATGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAg                    >  1:260720/4‑144 (MQ=255)
                                                                                                                                        atcgCCCGTAATATGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAg                    >  1:1165132/4‑144 (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)
                                                                                                                                                   tACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCATACATTGCTGCGATAGGTCGCATCAAGACAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCTAAACGTTGATCCCATGACAGTTCANCCCNNTnacgac   >  2:227579/1‑150 (MQ=255)
                                                                                                                                                    aCGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGACACGTTGATCCCATGCCAGTTCAGCGCCTTGacgacg  >  1:369653/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCCCATATAACGCGCTTCCGTGCCTGCATACCAACCATCTTCCGGTACATAACCTATCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACG  >  NC_000913/1521133‑1521430

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: