Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,308,636 Δ1 bp 11.4% coding (289/885 nt) nlpI ← lipoprotein NlpI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,308,6360C.11.4% 243.8 / 26.1 79coding (289/885 nt)nlpIlipoprotein NlpI
Reads supporting (aligned to +/- strand):  ref base C (35/34);  new base . (4/5);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

CACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCA  >  NC_000913/3308489‑3308785
                                                                                                                                                    |                                                                                                                                                     
cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                                     >  1:21482/1‑150 (MQ=255)
cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                                     >  2:1148367/1‑150 (MQ=255)
cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                                     >  1:65951/1‑150 (MQ=255)
 aCGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAgg                                                                                                                                                    <  2:858843/150‑1 (MQ=255)
 aCGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAgg                                                                                                                                                    <  2:720013/150‑1 (MQ=255)
   gACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCa                                                                                                                                                  <  1:574883/150‑1 (MQ=255)
   gACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCa                                                                                                                                                  <  1:1240671/150‑1 (MQ=255)
     ccgccgTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCata                                                                                                                                                <  1:550843/150‑1 (MQ=255)
        ccgTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCatatat                                                                                                                                             <  2:1328677/150‑3 (MQ=255)
            aaTATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTc                                                                                                                                         <  2:831540/150‑1 (MQ=255)
            aaTATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTc                                                                                                                                         <  1:229694/150‑1 (MQ=255)
             atatAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCatatatttaac                                                                                                                                        <  1:1273598/150‑8 (MQ=38)
                taATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGAt                                                                                                                                     <  2:653981/150‑1 (MQ=255)
                 aaTGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGAtt                                                                                                                                    >  2:537723/1‑150 (MQ=255)
                   tGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGc                                                                                                                                  >  1:209536/1‑150 (MQ=255)
                    gCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGTATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGcc                                                                                                                                 <  2:1149856/150‑1 (MQ=255)
                       ttCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTTGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAgc                                                                                                                              <  2:1236224/149‑1 (MQ=255)
                         cccGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAgcgc                                                                                                                            >  2:1207549/1‑150 (MQ=255)
                             cgATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGc                                                                                                                        >  2:423414/1‑150 (MQ=255)
                                ttCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGaa                                                                                                                     >  2:787509/1‑150 (MQ=255)
                                        cgcgTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGtt                                                                                                             <  1:1033988/150‑1 (MQ=255)
                                        cgcgTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGtt                                                                                                             <  1:890761/150‑1 (MQ=255)
                                         gagTAGTTCTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTa                                                                                                            >  2:772051/3‑150 (MQ=255)
                                          cgTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTAc                                                                                                           <  2:201655/150‑1 (MQ=255)
                                          cgTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTAc                                                                                                           <  1:1197631/150‑1 (MQ=255)
                                               ttGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCt                                                                                                      <  2:181539/150‑1 (MQ=255)
                                                tGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTa                                                                                                     >  2:750235/1‑150 (MQ=255)
                                                        ggATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCata                                                                                                                                                <  2:1176308/99‑1 (MQ=25)
                                                        ggATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCata                                                                                                                                                >  1:1176308/1‑99 (MQ=25)
                                                            cAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGAc                                                                                         >  2:101300/1‑150 (MQ=255)
                                                              aGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCt                                                                                                      >  1:318781/1‑135 (MQ=255)
                                                              aGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCt                                                                                                      <  2:318781/135‑1 (MQ=255)
                                                              aGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTAAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCg                                                                                       >  2:70413/1‑150 (MQ=255)
                                                               gCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAGAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                                      >  2:720041/1‑150 (MQ=255)
                                                                cTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCgag                                                                                     <  2:21482/150‑1 (MQ=255)
                                                                cTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCgag                                                                                     >  1:830199/1‑150 (MQ=255)
                                                                  cAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGAc                                                                                   <  1:1153076/150‑1 (MQ=255)
                                                                  cAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGAc                                                                                   <  1:1148367/150‑1 (MQ=255)
                                                                   aaGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                                      >  2:577227/1‑146 (MQ=255)
                                                                   aaGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCga                                                                                      <  1:577227/146‑1 (MQ=255)
                                                                    aGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTa                                                                                 >  1:1391230/1‑150 (MQ=255)
                                                                      tACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATc                                                                               >  2:930378/1‑150 (MQ=255)
                                                                       aCAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCa                                                                              <  1:1029311/150‑1 (MQ=255)
                                                                       aCAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCa                                                                              >  2:1172567/1‑150 (MQ=255)
                                                                          gAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATac                                                                           >  1:106886/1‑150 (MQ=255)
                                                                          gAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATac                                                                           <  1:500671/150‑1 (MQ=255)
                                                                                aaCGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCACACCGAGACTATCATACAACACt                                                                     >  2:1247689/1‑150 (MQ=255)
                                                                                    cTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCcgc                                                                 <  2:127114/150‑1 (MQ=255)
                                                                                    cTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCcgc                                                                 <  1:1207549/150‑1 (MQ=255)
                                                                                    cTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCcgc                                                                 >  2:64101/1‑150 (MQ=255)
                                                                                      tCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCgcgc                                                               >  1:93406/1‑150 (MQ=255)
                                                                                       cATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCt                                                              <  1:1172567/150‑1 (MQ=255)
                                                                                             cagcaTCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCAtata                                                        <  2:65951/150‑1 (MQ=255)
                                                                                                caTCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATaaaa                                                     >  1:44436/1‑150 (MQ=255)
                                                                                                caTCAAAATTGCCTGCCCGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCa                                                            >  1:958308/1‑143 (MQ=255)
                                                                                                caTCAAAATTGCCTGCCCGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCa                                                            <  2:958308/143‑1 (MQ=255)
                                                                                                  tCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGc                                                   >  1:1363698/1‑150 (MQ=255)
                                                                                                      aaTTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCtgtg                                               <  1:423414/150‑1 (MQ=255)
                                                                                                        ttGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTgcg                                             >  2:506137/1‑150 (MQ=255)
                                                                                                             tgcctgcGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTc                                        >  2:987357/1‑150 (MQ=255)
                                                                                                              gcctgcGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCg                                       <  1:950428/150‑1 (MQ=255)
                                                                                                                    gTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCg                                 >  2:1275505/1‑150 (MQ=255)
                                                                                                                    gTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCg                                 >  2:1019511/1‑150 (MQ=255)
                                                                                                                       ccggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAg                                                    >  1:600035/5‑128 (MQ=37)
                                                                                                                       ccggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAg                                                    <  2:600035/124‑1 (MQ=37)
                                                                                                                        cggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTa                             >  2:1090129/4‑150 (MQ=38)
                                                                                                                        cggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTa                             >  2:1048397/4‑150 (MQ=38)
                                                                                                                        aaTATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTa                             <  1:772051/150‑1 (MQ=255)
                                                                                                                              tGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGccc                       >  2:1391701/1‑150 (MQ=255)
                                                                                                                                 cT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGAc                    >  2:1343525/1‑150 (MQ=255)
                                                                                                                                       tAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCa               >  2:276070/1‑150 (MQ=255)
                                                                                                                                         aTTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAg             >  1:1273462/1‑150 (MQ=255)
                                                                                                                                             aaTACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGAtt         >  1:140009/1‑150 (MQ=255)
                                                                                                                                               tACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTg       <  1:787509/150‑1 (MQ=255)
                                                                                                                                               tACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTg       <  1:70413/150‑1 (MQ=255)
                                                                                                                                               tACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTg       <  2:1273462/150‑1 (MQ=255)
                                                                                                                                                 cTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGtt     <  1:537723/150‑1 (MQ=255)
                                                                                                                                                 cTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGtt     >  2:1399408/1‑150 (MQ=255)
                                                                                                                                                  ttCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTc    <  1:1399408/150‑1 (MQ=255)
                                                                                                                                                    cAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGt      >  2:616152/1‑146 (MQ=255)
                                                                                                                                                    cAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGt      <  1:616152/146‑1 (MQ=255)
                                                                                                                                                    cAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCa  <  1:54032/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCA  >  NC_000913/3308489‑3308785

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