Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,732,811 G→T S1242I (AGC→ATC)  lhr → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,732,8110GT98.8% 328.2 / ‑3.5 86S1242I (AGC→ATC) lhrputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (44/41);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

GATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGTTCCAACGCACGCAC  >  NC_000913/1732663‑1732958
                                                                                                                                                    |                                                                                                                                                   
gATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATcc                                                                                                                                                   <  2:865772/151‑1 (MQ=255)
gATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATcc                                                                                                                                                   >  1:422073/1‑151 (MQ=255)
gATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGGTCTGGCTGACGGAGGAGCCTGGCTTGCACAACAAATCATcc                                                                                                                                                   >  1:279507/1‑151 (MQ=255)
 aTGACGGCCTGGTGGCCCTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCa                                                                                                                                                  <  2:840444/151‑1 (MQ=255)
 aTGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCa                                                                                                                                                  <  1:38709/151‑1 (MQ=255)
  tGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCa                                                                                                                                                  >  2:614527/1‑150 (MQ=255)
    aCGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACg                                                                                                                                               <  1:265484/151‑1 (MQ=255)
    aCGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACg                                                                                                                                               >  2:385223/1‑151 (MQ=255)
     cGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACgg                                                                                                                                              <  2:454092/151‑1 (MQ=255)
     cGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACgg                                                                                                                                              >  1:183122/1‑151 (MQ=255)
     cGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACg                                                                                                                                               <  2:1118535/150‑1 (MQ=255)
      ggCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACgg                                                                                                                                              >  2:421855/1‑150 (MQ=255)
      ggCCGGGTGGCACTGCATCGACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAAc                                                                                                                                                <  2:840436/148‑1 (MQ=255)
       gCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGAt                                                                                                                                            <  1:1257552/151‑1 (MQ=255)
       gCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGAt                                                                                                                                            <  2:432355/151‑1 (MQ=255)
         cTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATAc                                                                                                                                          >  1:287927/1‑151 (MQ=255)
              ggCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGAc                                                                                                                                     <  1:632086/151‑1 (MQ=255)
              ggCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGAc                                                                                                                                     <  1:507189/151‑1 (MQ=255)
                cACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACaa                                                                                                                                   <  1:659704/151‑1 (MQ=255)
                cACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACaa                                                                                                                                   <  1:539735/151‑1 (MQ=255)
                  cTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACaaa                                                                                                                                  >  2:564682/1‑150 (MQ=255)
                   tGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAAt                                                                                                                                >  1:1281943/1‑151 (MQ=255)
                    gCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATc                                                                                                                               <  2:206537/151‑1 (MQ=255)
                    gCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAAt                                                                                                                                <  2:868225/150‑1 (MQ=255)
                      aTCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCgg                                                                                                                             >  2:1132372/1‑151 (MQ=255)
                         tACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGCAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCg                                                                                                                           >  2:1071104/1‑150 (MQ=255)
                          aCAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGc                                                                                                                            <  1:299246/148‑1 (MQ=255)
                          aCAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGaa                                                                                                                         <  2:1239329/151‑1 (MQ=255)
                          aCAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGa                                                                                                                          <  1:947071/150‑1 (MQ=255)
                          aCAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGa                                                                                                                          >  2:196226/1‑150 (MQ=255)
                              gAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCgg                                                                                                                     >  2:464489/1‑151 (MQ=255)
                              gAATATGCCGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACTAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCgg                                                                                                                     <  1:962251/151‑1 (MQ=255)
                                 tatGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTg                                                                                                                  >  2:675633/1‑151 (MQ=255)
                                  atGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTg                                                                                                                  >  2:537030/1‑150 (MQ=255)
                                  atGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGa                                                                                                                 <  2:690537/151‑1 (MQ=255)
                                     ctgcAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGAtc                                                                                                               <  1:1131333/150‑1 (MQ=255)
                                      tgcAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGAtct                                                                                                              >  2:503703/1‑150 (MQ=255)
                                       gcAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGAtct                                                                                                              <  2:285210/149‑1 (MQ=255)
                                           aaTCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGc                                                                                                        <  1:598858/151‑1 (MQ=255)
                                              cGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGccc                                                                                                      <  2:894144/150‑1 (MQ=255)
                                              cGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCt                                                                                                     >  1:389393/1‑151 (MQ=255)
                                              cGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCt                                                                                                     >  2:719150/1‑151 (MQ=255)
                                               gTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCt                                                                                                     >  2:845945/1‑150 (MQ=255)
                                                  cACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCa                                                                                                  >  1:1311272/1‑150 (MQ=255)
                                                   aCTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCa                                                                                                  >  2:316350/1‑149 (MQ=255)
                                                      cccGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAgg                                                                                             >  1:1328330/1‑151 (MQ=255)
                                                         gccgAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGt                                                                                          >  1:214322/1‑151 (MQ=255)
                                                           cgAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGtt                                                                                         >  1:1044766/1‑150 (MQ=255)
                                                           cgAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGtt                                                                                         >  2:300504/1‑150 (MQ=255)
                                                           cgAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGtt                                                                                         >  1:440817/1‑150 (MQ=255)
                                                               gCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATgg                                                                                     <  2:437592/150‑1 (MQ=255)
                                                                cGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGgc                                                                                   <  1:376559/151‑1 (MQ=255)
                                                                 ggATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGgcg                                                                                  >  1:699439/1‑151 (MQ=255)
                                                                        gCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGCCAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCt                                                                           <  1:1071697/151‑1 (MQ=255)
                                                                          gAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTgg                                                                         <  1:188056/151‑1 (MQ=255)
                                                                           aaTCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGc                                                                        >  1:387095/1‑151 (MQ=255)
                                                                            aTCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCa                                                                       <  2:314181/151‑1 (MQ=255)
                                                                             tCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCa                                                                       >  2:842832/1‑150 (MQ=255)
                                                                              cGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCa                                                                       >  1:642179/1‑149 (MQ=255)
                                                                                ttCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGc                                                                   >  1:1042078/1‑151 (MQ=255)
                                                                                      gcTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGtcatc                                                             <  2:187465/151‑1 (MQ=255)
                                                                                        tGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCAc                                                           >  2:164123/1‑151 (MQ=255)
                                                                                           aacaaGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCAcc                                                          >  2:110596/1‑149 (MQ=255)
                                                                                                gCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACa                                                   <  1:130745/151‑1 (MQ=255)
                                                                                                  gATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACAtt                                                 <  1:514891/151‑1 (MQ=255)
                                                                                                         gCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCAc                                          <  2:438400/151‑1 (MQ=255)
                                                                                                               cTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTAc                                    >  1:259506/1‑151 (MQ=255)
                                                                                                                  gCTGACGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTAcgcg                                 >  2:471493/1‑151 (MQ=255)
                                                                                                                      aCGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGcc                               <  1:68379/149‑1 (MQ=255)
                                                                                                                       cGGAGGAGCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCt                             <  2:389382/150‑1 (MQ=255)
                                                                                                                         gaggagCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGAGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCAc                          >  2:409643/1‑151 (MQ=255)
                                                                                                                          aggagCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCAcc                         >  1:1289835/1‑151 (MQ=255)
                                                                                                                           ggagCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCAccc                        <  1:935401/151‑1 (MQ=255)
                                                                                                                           ggagCCTGGTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCAccc                        <  2:598748/151‑1 (MQ=255)
                                                                                                                                  ggTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCgcagca                  <  1:264672/150‑1 (MQ=255)
                                                                                                                                  ggTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCgcagca                  <  1:657546/150‑1 (MQ=255)
                                                                                                                                  ggTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCgcagca                  <  2:513664/150‑1 (MQ=255)
                                                                                                                                   gTTTGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGt                <  2:1302704/151‑1 (MQ=255)
                                                                                                                                    tttGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGtt               >  2:217750/1‑151 (MQ=255)
                                                                                                                                     ttGCACAACAAATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGTTc              >  2:78247/1‑151 (MQ=255)
                                                                                                                                            acaaATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGTTCCAacgca       <  1:1068911/151‑1 (MQ=255)
                                                                                                                                            acaaATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGTTCCAacgc        >  1:1292236/1‑150 (MQ=255)
                                                                                                                                             caaATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGTTCCAacgcac      >  1:780904/1‑151 (MQ=255)
                                                                                                                                             caaATCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGTTCCAacgcac      >  2:1004372/1‑151 (MQ=255)
                                                                                                                                                 tCATCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGTTCCAACgcacgcac  >  1:1270050/1‑151 (MQ=255)
                                                                                                                                                 tCATCCAACGGATACACGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGTTCCAACgcacgcac  >  1:501341/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATGACGGCCTGGTGGCACTGCATCTACAGGAATATGCTGCAGAATCGTTCACTCCCGCCGAAGCGGATCAGGCGAATCGTTCGGCGCTGCAACAAGCGATAGTCGCTGTTCTGGCTGACGGAGGAGCCTGGTTTGCACAACAAATCAGCCAACGGATACGCGACAAAATCGGCGAATCGGTTGATCTCTCTGCCCTGCAAGAGGCGTTATGGGCGCTGGTCTGGCAAGGCGTCATCACCAGCGACATTTGGGCACCGTTACGCGCCCTCACCCGCAGCAGTTCCAACGCACGCAC  >  NC_000913/1732663‑1732958

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