Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 962,473 A→T D160V (GAT→GTT)  rpsA → 30S ribosomal subunit protein S1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913962,4730AT98.8% 304.7 / ‑3.7 81D160V (GAT→GTT) rpsA30S ribosomal subunit protein S1
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (45/35);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCACTGACTACGGTGCATTCG  >  NC_000913/962327‑962619
                                                                                                                                                  |                                                                                                                                                  
tATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTc                                                                                                                                                  >  2:99060/1‑149 (MQ=255)
tATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTc                                                                                                                                                  >  1:397446/1‑149 (MQ=255)
  tCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGa                                                                                                                                               <  1:173008/150‑1 (MQ=255)
   cAACGGCAAAGTTAAGCGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAg                                                                                                                                             <  1:1031991/151‑1 (MQ=255)
     aCGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCTCCTGGAAGGCAAAGAGCTTGAATTTAAAGTAAGCAAGCTGGTTCAGAAgcg                                                                                                                                           <  2:893138/151‑1 (MQ=255)
     aCGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAgcg                                                                                                                                           >  2:405252/1‑151 (MQ=255)
         cAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGcaac                                                                                                                                       <  1:1133946/151‑1 (MQ=255)
          aaaGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGcaac                                                                                                                                       >  1:533486/1‑150 (MQ=255)
          aaaGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGcaac                                                                                                                                       >  2:666945/1‑150 (MQ=255)
           aaGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGcaaca                                                                                                                                      >  2:612342/1‑150 (MQ=255)
             gTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACg                                                                                                                                   >  1:163144/1‑151 (MQ=255)
             gTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACg                                                                                                                                   <  2:966481/151‑1 (MQ=255)
                aaGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACgttg                                                                                                                                <  2:287002/151‑1 (MQ=255)
                   ggcggcTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACgttgttg                                                                                                                             <  1:748861/151‑1 (MQ=255)
                     cggcTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACgttgttgtt                                                                                                                           >  1:485536/1‑151 (MQ=255)
                      ggcTTCACTGTTGAGCTGAGCGGTATTCGTGCGTTCCTGCCAGGTTCACTGGTAGACGTTCGTCCGTTGCGTGACACTCTGCACCTGGAAGGCAGAGGGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAGGCGCAACACTgatgttgt                                                                                                                            <  1:767705/149‑1 (MQ=255)
                      ggcTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGttt                                                                                                                          <  1:1140811/151‑1 (MQ=255)
                      ggcTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGttt                                                                                                                          >  1:679702/1‑151 (MQ=255)
                         ttCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTtctc                                                                                                                       >  1:747100/1‑151 (MQ=255)
                            aCTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCtcgtc                                                                                                                    <  2:400886/151‑1 (MQ=255)
                               gTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTcgtcgt                                                                                                                  <  1:132616/150‑1 (MQ=255)
                                ttGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGc                                                                                                                >  2:364419/1‑151 (MQ=255)
                                 tGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGc                                                                                                                <  1:1105415/150‑1 (MQ=255)
                                    gCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTTCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGtt                                                                                                            <  1:355997/151‑1 (MQ=255)
                                      tGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTAt                                                                                                          >  1:682088/1‑151 (MQ=255)
                                        aaCGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCg                                                                                                        >  1:357544/1‑151 (MQ=255)
                                            gTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATc                                                                                                    <  2:379794/151‑1 (MQ=255)
                                               ttCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCg                                                                                                  >  2:1040892/1‑150 (MQ=255)
                                                 cGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGaaa                                                                                               <  2:1061469/151‑1 (MQ=255)
                                                 cGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGaaa                                                                                               <  1:396861/151‑1 (MQ=255)
                                                 cGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGaa                                                                                                >  2:1040407/1‑150 (MQ=255)
                                                  gTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGaaa                                                                                               >  2:538762/1‑150 (MQ=255)
                                                   tGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAAc                                                                                             >  2:10270/1‑151 (MQ=255)
                                                      gTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAgc                                                                                          <  2:821826/151‑1 (MQ=255)
                                                      gTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAgc                                                                                          <  2:272760/151‑1 (MQ=255)
                                                        tCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAgcgc                                                                                        >  1:708091/1‑151 (MQ=255)
                                                        tCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAgcgc                                                                                        >  1:580651/1‑151 (MQ=255)
                                                         ccTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAgcgc                                                                                        <  1:221369/150‑1 (MQ=255)
                                                         ccTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCa                                                                                       <  1:970480/151‑1 (MQ=255)
                                                             ccAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAgc                                                                                   >  2:1112510/1‑151 (MQ=255)
                                                              cAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAgcg                                                                                  >  2:499552/1‑151 (MQ=255)
                                                                     tctGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATc                                                                             >  1:953796/1‑149 (MQ=255)
                                                                       tGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAgct                                                                         >  1:144616/1‑151 (MQ=255)
                                                                        ggTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAgctg                                                                        <  2:183002/151‑1 (MQ=255)
                                                                        ggTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAgct                                                                         <  1:250027/150‑1 (MQ=255)
                                                                               gTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGaa                                                                  >  2:459983/1‑150 (MQ=255)
                                                                               gTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGaa                                                                  >  2:122309/1‑150 (MQ=255)
                                                                                     ccGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTg                                                            <  2:583168/150‑1 (MQ=255)
                                                                                      cGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGc                                                           >  1:118419/1‑150 (MQ=255)
                                                                                        gTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAgg                                                        <  2:485561/151‑1 (MQ=255)
                                                                                          gCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGa                                                       <  2:397465/150‑1 (MQ=255)
                                                                                           cGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGa                                                       >  1:1094372/1‑149 (MQ=255)
                                                                                            gTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAgg                                                    >  1:85875/1‑151 (MQ=255)
                                                                                             tGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGCAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGc                                                   >  2:613743/1‑151 (MQ=255)
                                                                                                     tGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGt                                           >  1:1148623/1‑151 (MQ=255)
                                                                                                     tGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGt                                           >  1:724485/1‑151 (MQ=255)
                                                                                                      gCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGtt                                          >  2:987811/1‑151 (MQ=255)
                                                                                                          cttgAAGGCAAAAAGCTTGAAATTAAAGTAATCAAGCTGGTTCAGAAGGGGAACAACCGTGTTTGTTCCCGGCGGGGCGTTATCGAATCCGAAAACCGCGCAGAGGGCGAACCGCCGCTGGAAAACCTGCAGGAAGGCATGGAAGtt                                          <  2:822609/144‑1 (MQ=255)
                                                                                                            ggAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAgg                                     <  2:900993/150‑1 (MQ=255)
                                                                                                              aaGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTAt                                  <  2:1005076/151‑1 (MQ=255)
                                                                                                                ggCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCg                                <  2:565983/151‑1 (MQ=255)
                                                                                                                ggCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCg                                <  2:169009/151‑1 (MQ=255)
                                                                                                                 gCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGt                               <  1:769599/151‑1 (MQ=255)
                                                                                                                  cAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGtt                              >  1:664675/1‑151 (MQ=255)
                                                                                                                  cAAAGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGtt                              <  2:456957/151‑1 (MQ=255)
                                                                                                                    aaGAGCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTa                             >  2:823999/1‑150 (MQ=255)
                                                                                                                       agCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGa                          >  2:207337/1‑150 (MQ=255)
                                                                                                                       agCTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGa                          >  2:591533/1‑150 (MQ=255)
                                                                                                                         cTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAAcc                       >  1:983277/1‑151 (MQ=255)
                                                                                                                         cTTGAATTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAAcc                       >  2:465385/1‑151 (MQ=255)
                                                                                                                             aaTTTAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCAc                   >  1:400954/1‑151 (MQ=255)
                                                                                                                                ttAAAGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCACTg                 >  1:323740/1‑150 (MQ=255)
                                                                                                                                  aaaGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCACTGAc               <  2:441606/150‑1 (MQ=255)
                                                                                                                                    aGTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCACTGACTAc            <  1:612316/151‑1 (MQ=255)
                                                                                                                                     gTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCACTGACTACg           <  2:682120/151‑1 (MQ=255)
                                                                                                                                     gTAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCACTGACTACg           >  2:303449/1‑151 (MQ=255)
                                                                                                                                      tAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCACTGACTACgg          >  1:364203/1‑151 (MQ=255)
                                                                                                                                      tAATCAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCACTGACTACgg          >  1:852969/1‑151 (MQ=255)
                                                                                                                                          cAAGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATTGTTAAGAACCTCACTGACTACGGTGc       >  2:468233/1‑150 (MQ=255)
                                                                                                                                            aGCTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCACTGACTACGGTGCAtt    <  1:185020/151‑1 (MQ=255)
                                                                                                                                              cTGGTTCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCACTGACTACGGTGCATTCg  >  1:846833/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TATCAACGGCAAAGTTAAGGGCGGCTTCACTGTTGAGCTGAACGGTATTCGTGCGTTCCTGCCAGGTTCTCTGGTAGACGTTCGTCCGGTGCGTGACACTCTGCACCTGGAAGGCAAAGAGCTTGAATTTAAAGTAATCAAGCTGGATCAGAAGCGCAACAACGTTGTTGTTTCTCGTCGTGCCGTTATCGAATCCGAAAACAGCGCAGAGCGCGATCAGCTGCTGGAAAACCTGCAGGAAGGCATGGAAGTTAAAGGTATCGTTAAGAACCTCACTGACTACGGTGCATTCG  >  NC_000913/962327‑962619

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