Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,733,986 T→C D63G (GAT→GGT)  clpB ← protein disaggregation chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,733,9860TC99.0% 374.5 / NA 97D63G (GAT→GGT) clpBprotein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (62/34);  total (63/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATATCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTGGGCATCGGCA  >  NC_000913/2733843‑2734132
                                                                                                                                               |                                                                                                                                                  
acAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTgcg                                                                                                                                             >  4:371055/1‑151 (MQ=255)
 cAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTgcgc                                                                                                                                            >  3:253247/1‑151 (MQ=255)
  aGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTgcgc                                                                                                                                            >  2:293844/1‑150 (MQ=255)
  aGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTgcgc                                                                                                                                            >  2:441925/1‑150 (MQ=255)
    ttCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAAc                                                                                                                                         >  2:1020752/1‑151 (MQ=255)
     tCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACt                                                                                                                                        >  4:589610/1‑151 (MQ=255)
       tGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTgg                                                                                                                                      >  2:1585306/1‑151 (MQ=255)
        gACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGc                                                                                                                                     >  4:333336/1‑151 (MQ=255)
           gagaTAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGcc                                                                                                                                    >  4:155051/1‑149 (MQ=255)
           gagaTAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGcc                                                                                                                                    >  4:1029234/1‑149 (MQ=255)
           gagaTAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGcc                                                                                                                                    >  2:1467668/1‑149 (MQ=255)
               tAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCAtt                                                                                                                              >  4:1225360/1‑151 (MQ=255)
               tAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCAtt                                                                                                                              >  2:1484627/1‑151 (MQ=255)
                aaaGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCAttt                                                                                                                             >  2:280470/1‑151 (MQ=255)
                    ttATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGc                                                                                                                         >  4:1043191/1‑151 (MQ=255)
                    ttATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGc                                                                                                                         >  1:245195/1‑151 (MQ=255)
                      aTCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGcc                                                                                                                        >  1:330138/1‑150 (MQ=255)
                      aTCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGcc                                                                                                                        <  4:118117/150‑1 (MQ=255)
                         accacGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGCAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGc                                                                                                                     >  2:136636/1‑150 (MQ=255)
                          ccacGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCgg                                                                                                                   >  3:263349/1‑151 (MQ=255)
                          ccacGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCgg                                                                                                                   >  2:1211498/1‑151 (MQ=255)
                              gTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATg                                                                                                                >  2:1429855/1‑150 (MQ=255)
                              gTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGt                                                                                                               >  1:463244/1‑151 (MQ=255)
                                ttttGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGtt                                                                                                              <  2:963923/150‑1 (MQ=255)
                                ttttGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGAGGtt                                                                                                              >  1:102801/1‑150 (MQ=255)
                                gattGCGCCAGCGTGGGGCAAGGAAAAAGAGGGGGCACCAGATGTTGTGATGGCTGGGCACCGCCACCATTAAAAAAAAACTGACGTAAACGATTAAATGCCTGATTGATATCTGTGCGCAACTGGCAAGCATTTATGCCAGCGGATg                                                                                                                >  2:628080/3‑148 (MQ=255)
                                 tttGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTta                                                                                                             >  4:1595996/1‑150 (MQ=255)
                                   tGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTtaat                                                                                                           >  2:988784/1‑150 (MQ=255)
                                   tGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTtaat                                                                                                           >  2:664244/1‑150 (MQ=255)
                                    gcgcCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTtaata                                                                                                          <  3:49922/150‑1 (MQ=255)
                                      gcCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTGAATGCCTGATTGATACCTGTGCTCAAAGGGCCAGCATTTATGGCAGCGGATGTtaataa                                                                                                         <  4:1597007/149‑1 (MQ=255)
                                          gCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACt                                                                                                   >  3:438493/1‑151 (MQ=255)
                                           cTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACt                                                                                                   >  3:251073/1‑150 (MQ=255)
                                            ttGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTa                                                                                                  >  2:512816/1‑150 (MQ=255)
                                             tGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTaa                                                                                                 <  2:602221/150‑1 (MQ=255)
                                             tGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTa                                                                                                  <  1:300775/149‑1 (MQ=255)
                                               tCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCg                                                                                              >  2:568850/1‑151 (MQ=255)
                                               tCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCg                                                                                              >  2:66517/1‑151 (MQ=255)
                                               tCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCg                                                                                              >  2:847485/1‑151 (MQ=255)
                                               tCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCg                                                                                              >  2:1230268/1‑151 (MQ=255)
                                                cGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCg                                                                                              <  2:588794/150‑1 (MQ=255)
                                                   aaaGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAAcc                                                                                          >  4:572935/1‑151 (MQ=255)
                                                     aGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAAcccc                                                                                        >  2:1254704/1‑151 (MQ=255)
                                                     aGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAAcccc                                                                                        >  1:76134/1‑151 (MQ=255)
                                                     aGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAAcccc                                                                                        <  4:886171/151‑1 (MQ=255)
                                                       aTTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGAc                                                                                                    <  3:151565/137‑1 (MQ=255)
                                                       aTTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGAc                                                                                                    >  1:151551/1‑137 (MQ=255)
                                                       aTTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCt                                                                                      <  4:394385/151‑1 (MQ=255)
                                                            gAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTg                                                                                 <  4:1215065/151‑1 (MQ=255)
                                                            gAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTg                                                                                 <  2:25167/151‑1 (MQ=255)
                                                            gAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTg                                                                                 <  1:378/151‑1 (MQ=255)
                                                              aCGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGAt                                                                               >  1:327748/1‑151 (MQ=255)
                                                               cgcgCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGAtt                                                                              >  4:836836/1‑151 (MQ=255)
                                                                  gcACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTc                                                                             >  1:258928/1‑149 (MQ=255)
                                                                        gATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGgc                                                                     >  4:295207/1‑151 (MQ=255)
                                                                           ccTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCt                                                                  <  1:334788/151‑1 (MQ=255)
                                                                           ccTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCt                                                                  >  4:429798/1‑151 (MQ=255)
                                                                              gtgATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCAt                                                               <  2:709777/151‑1 (MQ=255)
                                                                                gATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCAtt                                                              >  4:496131/1‑150 (MQ=255)
                                                                                gATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCAtt                                                              <  4:684069/150‑1 (MQ=255)
                                                                                gATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCAtt                                                              <  2:774782/150‑1 (MQ=255)
                                                                                 aTGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTaa                                                            <  2:215887/151‑1 (MQ=255)
                                                                                   ggCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAAt                                                          <  2:1502792/151‑1 (MQ=255)
                                                                                   ggCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAAt                                                          >  2:1671229/1‑151 (MQ=255)
                                                                                     cTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATg                                                         >  2:1146854/1‑150 (MQ=255)
                                                                                       ggACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGa                                                        <  2:935368/149‑1 (MQ=255)
                                                                                         aCATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTg                                                    <  4:955979/151‑1 (MQ=255)
                                                                                         aCATCACCACCAGTAACTTCAACCTGCGGTAAACGATTTAATGCATGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCGAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCTTGATTGAGTAGGTCGCTCATTAAATGAAGt                                                     >  4:1260240/1‑150 (MQ=255)
                                                                                          cATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTgg                                                   <  2:370502/151‑1 (MQ=255)
                                                                                          cATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTgg                                                   >  4:1615142/1‑151 (MQ=255)
                                                                                           aTCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGt                                                  <  4:1275414/151‑1 (MQ=255)
                                                                                           aTCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGt                                                  >  1:442281/1‑151 (MQ=255)
                                                                                              accaccAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCg                                               >  1:163539/1‑151 (MQ=255)
                                                                                              accaccAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCg                                               >  4:294287/1‑151 (MQ=255)
                                                                                                  ccaGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATaa                                           <  2:1285426/151‑1 (MQ=255)
                                                                                                           tCAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGt                                  >  1:58757/1‑151 (MQ=255)
                                                                                                            cAACCTGCTGTAAACGATTGAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTg                                   <  2:1175087/149‑1 (MQ=255)
                                                                                                            cAACCTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGt                                  <  4:257069/150‑1 (MQ=255)
                                                                                                               ccTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTg                              >  2:1564899/1‑151 (MQ=255)
                                                                                                               ccTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTg                              >  4:440386/1‑151 (MQ=255)
                                                                                                               ccTGCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTg                              >  2:143465/1‑151 (MQ=255)
                                                                                                                  gCGGTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGccc                           >  2:1177962/1‑151 (MQ=255)
                                                                                                                    ggTAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGa                         <  4:670428/151‑1 (MQ=255)
                                                                                                                      tAAACGATTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGt                       >  4:260493/1‑151 (MQ=255)
                                                                                                                            aTTTAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGt                 <  4:280916/151‑1 (MQ=255)
                                                                                                                              ttAATGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTg                >  1:273226/1‑150 (MQ=255)
                                                                                                                                aaTGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGAtt             <  1:436142/151‑1 (MQ=255)
                                                                                                                                 aTGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTg            <  4:1182897/151‑1 (MQ=255)
                                                                                                                                 aTGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTg            <  2:154788/151‑1 (MQ=255)
                                                                                                                                 aTGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTg            >  3:333731/1‑151 (MQ=255)
                                                                                                                                  tGCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTgg           >  2:268218/1‑151 (MQ=255)
                                                                                                                                   gCCTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTgg           <  2:471314/150‑1 (MQ=255)
                                                                                                                                    ccTGATTGATACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTGGGc         <  4:1439445/151‑1 (MQ=255)
                                                                                                                                       gattgatACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTGGGCa        <  4:136866/149‑1 (MQ=255)
                                                                                                                                       gattgatACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTGGGCATc      >  1:467428/1‑151 (MQ=255)
                                                                                                                                          tgatACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTGGGCATCGGc   >  2:1185381/1‑151 (MQ=255)
                                                                                                                                           gatACCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTGGGCATCGGc   >  4:1443152/1‑150 (MQ=255)
                                                                                                                                            atacCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTGGGCATCGGCa  >  2:1568031/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
ACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTAATGCCTGATTGATATCTGTGCGCAACTGGCCAGCATTTATGCCAGCGGATGTTAATAAAGGACTAACCGAACCCCCTTCCTGATTCAGCAGGGCGCTCATTAAATGAAGTGGTTCGATAAATTGGTTGTCGTGCCCGAGTGCAAGTGATTGGGCATCGGCA  >  NC_000913/2733843‑2734132

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