Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,947,095 G→A R28R (CGC→CGT mltA ← membrane‑bound lytic murein transglycosylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,0950GA99.1% 430.7 / NA 113R28R (CGC→CGTmltAmembrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (53/59);  minor base C (1/0);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 114 total aligned reads displayed.

TCTGGTTGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCGCGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGATAACAAACCGCCGCGTTTAATGA  >  NC_000913/2946947‑2947241
                                                                                                                                                    |                                                                                                                                                  
tctGGTTGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCcacg                                                                                                                                                  >  2:1145710/1‑151 (MQ=255)
  tGGTTGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGa                                                                                                                                                 >  1:336218/1‑150 (MQ=255)
  tGGTTGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGAt                                                                                                                                                >  4:342699/1‑151 (MQ=255)
   ggTTGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGa                                                                                                                                                 <  2:1514134/149‑1 (MQ=255)
   ggTTGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGa                                                                                                                                                 <  2:1442633/149‑1 (MQ=255)
   ggTTGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGa                                                                                                                                                 <  2:736512/149‑1 (MQ=255)
   ggTTGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGa                                                                                                                                                 <  4:1403816/149‑1 (MQ=255)
   ggTTGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGAt                                                                                                                                                <  3:424819/150‑1 (MQ=255)
     ttGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCgg                                                                                                                                             >  2:316748/1‑151 (MQ=255)
      tGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCggt                                                                                                                                            >  2:317080/1‑151 (MQ=255)
         cATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCggttgg                                                                                                                                         >  4:1583168/1‑151 (MQ=255)
          aTACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCggttggt                                                                                                                                        <  4:914907/151‑1 (MQ=255)
            aCAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGttt                                                                                                                                      <  4:111117/151‑1 (MQ=255)
               gACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTgg                                                                                                                                    >  4:653870/1‑150 (MQ=255)
                 cgcgGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGa                                                                                                                                   >  4:1625406/1‑149 (MQ=255)
                  gcgGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAaga                                                                                                                                <  4:1088976/151‑1 (MQ=255)
                    ggTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGc                                                                                                                              >  4:304452/1‑151 (MQ=255)
                       gACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAgg                                                                                                                           >  4:1423728/1‑151 (MQ=255)
                       gACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAg                                                                                                                            <  2:886481/150‑1 (MQ=255)
                       gACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAg                                                                                                                            <  2:582331/150‑1 (MQ=255)
                       gACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGGTGGTTTGGAAGAGCAgg                                                                                                                           >  2:1143926/1‑151 (MQ=255)
                        aCGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAggc                                                                                                                          <  4:316945/151‑1 (MQ=255)
                          gAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAggcgg                                                                                                                        >  4:1239047/1‑151 (MQ=255)
                           aaTTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAggcggc                                                                                                                       <  4:758955/151‑1 (MQ=255)
                              tACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCa                                                                                                                      >  4:1125414/1‑149 (MQ=255)
                               aCGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGCTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGc                                                                                                                   <  2:255867/151‑1 (MQ=255)
                                 ggATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCAt                                                                                                                 <  4:932387/151‑1 (MQ=255)
                                 ggATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCAt                                                                                                                 >  1:110787/1‑151 (MQ=255)
                                 ggATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATCGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCAt                                                                                                                 >  1:114535/1‑151 (MQ=255)
                                   atatGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTg                                                                                                               >  4:1279998/1‑151 (MQ=255)
                                     atGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGcc                                                                                                             <  4:1086863/151‑1 (MQ=255)
                                       ggTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCcac                                                                                                           >  1:464277/1‑151 (MQ=255)
                                        gTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCcaca                                                                                                          <  2:745522/151‑1 (MQ=255)
                                        gTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCcaca                                                                                                          <  4:216762/151‑1 (MQ=255)
                                           aaTTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAAc                                                                                                        <  2:447145/150‑1 (MQ=255)
                                             tttGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGt                                                                                                     >  4:1381345/1‑151 (MQ=255)
                                              ttGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGt                                                                                                     <  2:722893/150‑1 (MQ=255)
                                               tGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGc                                                                                                   >  4:1037801/1‑151 (MQ=255)
                                                 cTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGccc                                                                                                 >  2:1637373/1‑151 (MQ=255)
                                                   cGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGccc                                                                                                 <  3:349589/149‑1 (MQ=255)
                                                    ggCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGGTCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCAt                                                                                               >  4:1272638/1‑150 (MQ=255)
                                                       aaaaTCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATaaga                                                                                           <  1:17794/151‑1 (MQ=255)
                                                       aaaaTCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATaaga                                                                                           <  2:1059187/151‑1 (MQ=255)
                                                       aaaaTCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATa                                                                                              >  4:817522/1‑148 (MQ=255)
                                                       aaaaTCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATa                                                                                              >  2:282787/1‑148 (MQ=255)
                                                         aaTCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATaagaag                                                                                         <  4:1000329/151‑1 (MQ=255)
                                                            cACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGt                                                                                       <  4:584622/150‑1 (MQ=255)
                                                             aCCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTa                                                                                      >  2:13821/1‑150 (MQ=255)
                                                             aCCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTAc                                                                                     <  2:380226/151‑1 (MQ=255)
                                                               cGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGATGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACtt                                                                                   >  4:1469260/1‑151 (MQ=255)
                                                                ggCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACttt                                                                                  <  4:1151414/151‑1 (MQ=255)
                                                                 gCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTa                                                                                 <  2:219287/151‑1 (MQ=255)
                                                                         cGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACg                                                                          <  2:1011605/150‑1 (MQ=255)
                                                                         cGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGt                                                                         >  1:419128/1‑151 (MQ=255)
                                                                              cgcCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCttt                                                                    >  2:154136/1‑151 (MQ=255)
                                                                                ccAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCa                                                                  <  3:425741/151‑1 (MQ=255)
                                                                                  aaCGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCAtt                                                                <  4:43430/151‑1 (MQ=255)
                                                                                   aCGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTg                                                               >  2:44712/1‑151 (MQ=255)
                                                                                      gCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTtc                                                            <  2:652988/151‑1 (MQ=255)
                                                                                        aTCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTtctc                                                          >  1:73713/1‑151 (MQ=255)
                                                                                            ggctggTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCa                                                                                                >  4:485086/1‑109 (MQ=255)
                                                                                            ggctggTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCa                                                                                                <  2:484405/109‑1 (MQ=255)
                                                                                            ggctggTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCtcttc                                                       <  2:490634/150‑1 (MQ=255)
                                                                                                ggTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCttcttc                                                    <  2:720367/149‑1 (MQ=255)
                                                                                                  ttCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGt                                                >  4:766027/1‑151 (MQ=255)
                                                                                                   tCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGt                                                <  4:13849/150‑1 (MQ=255)
                                                                                                    cACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTa                                              <  2:1421948/151‑1 (MQ=255)
                                                                                                      ccAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTaa                                             <  1:138173/150‑1 (MQ=255)
                                                                                                      ccAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTaa                                             >  4:37953/1‑150 (MQ=255)
                                                                                                           gagaAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATa                                       <  3:306892/151‑1 (MQ=255)
                                                                                                            agaAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATa                                       >  2:1315229/1‑150 (MQ=255)
                                                                                                             gaAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATa                                       >  4:1576237/1‑149 (MQ=255)
                                                                                                             gaAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAg                                     <  4:1006502/151‑1 (MQ=255)
                                                                                                               aaGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGgc                                   <  2:571448/151‑1 (MQ=255)
                                                                                                               aaGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGgc                                   <  4:661414/151‑1 (MQ=255)
                                                                                                               aaGGCTGGGAAAATTTCCCGCTTTTTTATTGCTGTCCACGATCGGTGGGTTTGGAAGACCAGGGGCAAAGAATTGCACCAACCGTGCCAATAGGAGGGAACTTTACCCAACGTCTTTTCATTGTTCTCTTCTTCAGGTAAAAAATTAAg                                     >  1:444207/1‑149 (MQ=255)
                                                                                                                aGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGgcg                                  <  2:1568799/151‑1 (MQ=255)
                                                                                                                 ggCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTaaaaa                                          >  3:432612/1‑142 (MQ=255)
                                                                                                                 ggCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTaaaaa                                          <  1:432594/142‑1 (MQ=255)
                                                                                                                 ggCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGgcgc                                 <  4:654077/151‑1 (MQ=255)
                                                                                                                  gCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCa                                <  4:1234320/151‑1 (MQ=255)
                                                                                                                    tGGGTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAAc                              >  1:158869/1‑151 (MQ=255)
                                                                                                                      ggTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACg                             >  4:1611901/1‑150 (MQ=255)
                                                                                                                      ggTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGa                            <  2:367418/151‑1 (MQ=255)
                                                                                                                       gTAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGa                            >  2:690353/1‑150 (MQ=255)
                                                                                                                        tAAATTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGa                            >  4:748900/1‑149 (MQ=255)
                                                                                                                         aaaTTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAg                          <  1:156641/150‑1 (MQ=255)
                                                                                                                         aaaTTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAg                          <  2:1379625/150‑1 (MQ=255)
                                                                                                                          aaTTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGa                            >  1:436207/1‑147 (MQ=255)
                                                                                                                          aaTTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGAt                        <  4:1653511/151‑1 (MQ=255)
                                                                                                                           aTTTCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCAAAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGAt                        >  4:926710/1‑150 (MQ=255)
                                                                                                                             ttCCCGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGATaa                      <  2:1581216/150‑1 (MQ=255)
                                                                                                                                ccGTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGATAACaa                   >  4:663926/1‑150 (MQ=255)
                                                                                                                                  gTCTTTATATTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGATAACAAAcc                <  4:317223/151‑1 (MQ=255)
                                                                                                                                   tCTTTATATTGCTGTCTACGATCGGATGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGATAACAAAccg               >  2:1384654/1‑151 (MQ=255)
                                                                                                                                       tataTTGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGATAACAAACcgccgc           >  2:966938/1‑151 (MQ=255)
                                                                                                                                            tGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGATAACAAACCGCCGCGttt       >  4:1377347/1‑150 (MQ=255)
                                                                                                                                            tGCTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGATAACAAACCGCCGCGttt       <  4:342525/150‑1 (MQ=255)
                                                                                                                                              cTGTCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGATAACAAACCGCCGCGTTTaat    >  4:1452121/1‑151 (MQ=255)
                                                                                                                                                 tCCACGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGATAACAAACCGCCGCGTTTaatga  >  1:287634/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCTGGTTGCCATACAGACGCGGTGACGAATTACGGATATGGTTAATTTGCTCGGCAAAATCACCGGCGTTAATCGGCGCGCCAACGGCATCTGGCTGGTTCACCAGAGAGAAAGGCTGGGTAAATTTCCCGTCTTTATATTGCTGTCCGCGATCGGTTGGTTTGGAAGAGCAGGCGGCAAGCATTGCCACAACCGTGCCCATAAGAAGGTACTTTACCCAACGTCCTTTCATTGTTCTCTTCTTCAGGTTAAAAAATAAGGCGCAACGAAGATAACAAACCGCCGCGTTTAATGA  >  NC_000913/2946947‑2947241

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