Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,308,619:1 +G 11.1% coding (306/885 nt) nlpI ← lipoprotein NlpI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,308,6191.G11.1% 299.9 / 37.3 81coding (306/885 nt)nlpIlipoprotein NlpI
Reads supporting (aligned to +/- strand):  ref base . (37/35);  new base G (4/5);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

TCTTGCGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGAC  >  NC_000913/3308473‑3308767
                                                                                                                                                   |                                                                                                                                                    
tcTTGCGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAg                                                                                                                                                   <  2:1407262/150‑1 (MQ=255)
tcTTGCGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACCCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑aat                                                                                                                                                   >  2:753789/1‑149 (MQ=255)
 cTTGCGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGt                                                                                                                                                  >  2:229694/1‑150 (MQ=255)
 cTTGCGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGt                                                                                                                                                  <  2:198231/150‑1 (MQ=255)
   tGCGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTaa                                                                                                                                                >  1:341945/1‑150 (MQ=255)
     cgcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAAtt                                                                                                                                              <  2:485767/150‑1 (MQ=255)
      gcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTg                                                                                                                                             <  2:233199/150‑1 (MQ=255)
      gcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTg                                                                                                                                             >  2:890761/1‑150 (MQ=255)
      gcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTg                                                                                                                                             <  2:889566/150‑1 (MQ=255)
      gcTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTg                                                                                                                                             <  1:1265798/150‑1 (MQ=255)
       cTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGa                                                                                                                                            >  2:1029311/1‑150 (MQ=255)
             tGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACtt                                                                                                                                      >  2:950428/1‑150 (MQ=255)
                cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                   >  2:1148367/1‑150 (MQ=255)
                cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                   >  1:65951/1‑150 (MQ=255)
                cACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAg                                                                                                                                   >  1:21482/1‑150 (MQ=255)
                 aCGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAgg                                                                                                                                  <  2:720013/150‑1 (MQ=255)
                 aCGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAgg                                                                                                                                  <  2:858843/150‑1 (MQ=255)
                   gACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCa                                                                                                                                <  1:1240671/150‑1 (MQ=255)
                   gACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCa                                                                                                                                <  1:574883/150‑1 (MQ=255)
                     ccgccgTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCata                                                                                                                              <  1:550843/150‑1 (MQ=255)
                        ccgTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCTGAAGTATTCAATTACTT‑AGGCatatat                                                                                                                           <  2:1328677/150‑3 (MQ=255)
                            aaTATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTc                                                                                                                       <  1:229694/150‑1 (MQ=255)
                            aaTATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTc                                                                                                                       <  2:831540/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:890761/150‑1 (MQ=255)
                                                        cgcgTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGtt                                                                                           <  1:1033988/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                                                                                    <  2:318781/135‑1 (MQ=255)
                                                                              aGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCt                                                                                    >  1:318781/1‑135 (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                                                            >  2:1172567/1‑150 (MQ=255)
                                                                                       aCAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCa                                                            <  1:1029311/150‑1 (MQ=255)
                                                                                          gAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATac                                                         <  1:500671/150‑1 (MQ=255)
                                                                                          gAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATac                                                         >  1:106886/1‑150 (MQ=255)
                                                                                                aaCGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCACACCGAGACTATCATACAACACt                                                   >  2:1247689/1‑150 (MQ=255)
                                                                                                    cTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCcgc                                               <  1:1207549/150‑1 (MQ=255)
                                                                                                    cTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCcgc                                               >  2:64101/1‑150 (MQ=255)
                                                                                                    cTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCcgc                                               <  2:127114/150‑1 (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:1019511/1‑150 (MQ=255)
                                                                                                                                    gTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCg               >  2:1275505/1‑150 (MQ=255)
                                                                                                                                       ccggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAg                                  >  1:600035/5‑128 (MQ=37)
                                                                                                                                       ccggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAg                                  <  2:600035/124‑1 (MQ=37)
                                                                                                                                        cggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTa           >  2:1048397/4‑150 (MQ=38)
                                                                                                                                        cggatatGCCTGAAGTATTCAATTACTT‑AGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTa           >  2:1090129/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)
                                                                                                                                                   |                                                                                                                                                    
TCTTGCGCTAACTTGTCACGACCGCCGTAATATAATGCGATCCCGCGATTCAAGTGCGCGTAGTTGTAAGTTGGATCAAGCTCAAGTACAGAATCAAACGCTTCATAGGCAGCATCAAAATTGCCTGCCTGCGTTAAATATATGCCT‑AAGTAATTGAATACTTCAGGCATATCCGGTCGGATTGCCAGCGCTTGCGAAAAATCGTTACGCGCTAATGCCCTCAGACCGAGACTATCATACAACACTCCGCGCTCATATAAAAGCTGTGCGCGTTCGTCATCGGTTAAAGCCCGAC  >  NC_000913/3308473‑3308767

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: ‑