Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,637,295 C→T R21H (CGT→CAT)  der ← GTPase, multicopy suppressor of ftsJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,637,2950CT97.8% 274.9 / ‑2.1 89R21H (CGT→CAT) derGTPase, multicopy suppressor of ftsJ
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (43/44);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

TCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGG  >  NC_000913/2637163‑2637424
                                                                                                                                    |                                                                                                                                 
tcttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaa                                                                                                                             >  1:203091/1‑139 (MQ=255)
tcttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaa                                                                                                                             >  2:1510908/1‑139 (MQ=255)
tcttctGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTaa                                                                                                                             >  1:1245063/1‑139 (MQ=255)
       tgCCATCAATCCCGCGGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGt                                                                                                                      >  2:1131492/1‑139 (MQ=255)
          cATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                                                   >  2:1514589/1‑139 (MQ=255)
            tCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                                                  >  2:1649800/1‑138 (MQ=255)
            tCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                                                  <  1:1649800/138‑1 (MQ=255)
                                   gATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                          >  1:496710/1‑139 (MQ=255)
                                   gATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACa                                                                                          >  2:394051/1‑139 (MQ=255)
                                    aTAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACaa                                                                                         >  2:1345707/1‑139 (MQ=255)
                                    aTAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                                                   >  2:1171621/1‑129 (MQ=255)
                                    aTAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                                                   <  1:1171621/129‑1 (MQ=255)
                                     tAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAg                                                                                        >  1:991233/1‑139 (MQ=255)
                                     tAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAg                                                                                        >  2:715293/1‑139 (MQ=255)
                                     tAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAg                                                                                        >  1:118263/1‑139 (MQ=255)
                                       aaCTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAgcg                                                                                      >  2:1637435/1‑139 (MQ=255)
                                         cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                                                  <  2:926625/109‑1 (MQ=255)
                                         cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                                                  >  1:926625/1‑109 (MQ=255)
                                         cTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                    >  2:404262/1‑139 (MQ=255)
                                           cACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                                                  <  2:203091/139‑1 (MQ=255)
                                           cACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                                                  <  2:496710/139‑1 (MQ=255)
                                            aCGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAcc                                                                                 <  1:1510908/139‑1 (MQ=255)
                                             cGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACca                                                                                <  1:570396/139‑1 (MQ=255)
                                                  ttCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGt                                                                           >  2:400734/1‑139 (MQ=255)
                                                  ttCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGATAGGGCGCCCGACAAGCGCGACCACAGGt                                                                           <  1:1131492/139‑1 (MQ=255)
                                                    cAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                                                   >  1:809520/1‑97 (MQ=255)
                                                    cAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                                                   <  2:809520/97‑1 (MQ=255)
                                                    cAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAc                                                                         <  1:715293/139‑1 (MQ=255)
                                                     aaTTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAcc                                                                        <  2:1245063/139‑1 (MQ=255)
                                                      aTTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCa                                                                       <  1:219095/139‑1 (MQ=255)
                                                                     tACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCAt                                                        <  1:1345707/139‑1 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAAc                                                                                                                        >  1:381670/1‑68 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAAc                                                                                                                        <  2:381670/68‑1 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAAc                                                                                                                        <  2:150371/68‑1 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAAc                                                                                                                        >  1:150371/1‑68 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                                                  >  1:1390301/1‑74 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGAt                                                                                                                  <  2:1390301/74‑1 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAgcg                                                                                      >  2:1569644/1‑102 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAgcg                                                                                      <  1:1569644/102‑1 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                                                  <  1:870055/106‑1 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGAc                                                                                  >  2:870055/1‑106 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTa                                                                          <  2:592984/114‑1 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTa                                                                          >  1:592984/1‑114 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                                 >  2:711313/1‑139 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                                 >  2:1437851/1‑139 (MQ=255)
                                                                            ggTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                                 <  1:711313/139‑1 (MQ=255)
                                                                                aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                    >  1:1534971/1‑100 (MQ=255)
                                                                                aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGcgcg                                                                                    <  2:1534971/100‑1 (MQ=255)
                                                                                aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGt                                                                           >  1:1358172/1‑109 (MQ=255)
                                                                                aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGt                                                                           <  2:1358172/109‑1 (MQ=255)
                                                                                aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                                 >  2:992964/1‑135 (MQ=255)
                                                                                aCGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTtata                                                 <  1:992964/135‑1 (MQ=255)
                                                                                     tCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAg                                        <  1:394051/139‑1 (MQ=255)
                                                                                      cAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGa                                       <  1:404262/139‑1 (MQ=255)
                                                                                       aGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGAc                                      <  1:400734/139‑1 (MQ=255)
                                                                                            cGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgc                                 <  2:118263/139‑1 (MQ=255)
                                                                                             gggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca                                <  2:991233/139‑1 (MQ=255)
                                                                                              ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                                                   <  2:1241694/55‑1 (MQ=255)
                                                                                              ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGa                                                                                                                   >  1:1241694/1‑55 (MQ=255)
                                                                                              ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                              <  1:167391/60‑1 (MQ=255)
                                                                                              ggAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTc                                                                                                              >  2:167391/1‑60 (MQ=255)
                                                                                               gAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgcacg                              <  1:1514589/139‑1 (MQ=255)
                                                                                                   tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACcac                                                                               >  2:929074/1‑86 (MQ=255)
                                                                                                   tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACcac                                                                               <  1:929074/86‑1 (MQ=255)
                                                                                                   tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACa                                     <  2:865042/128‑1 (MQ=255)
                                                                                                   tCCGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACa                                     >  1:865042/1‑128 (MQ=255)
                                                                                                     cGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTa                                                                                                      <  1:600056/61‑1 (MQ=255)
                                                                                                     cGCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTa                                                                                                      >  2:600056/1‑61 (MQ=255)
                                                                                                      gCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACg                                                                                                                       <  2:347085/43‑1 (MQ=255)
                                                                                                      gCAACCAGCGCATCTCGGGTGCGAGTTAGATGGTTAAATAACg                                                                                                                       >  1:347085/1‑43 (MQ=255)
                                                                                                             gcgcATCTCGGGTGCGAGTTAGTTGGTTAAATAACGTGGATTTTCCTACGt                                                                                                        <  1:1049557/51‑1 (MQ=255)
                                                                                                             gcgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGt                                                                                                        >  2:1049557/1‑51 (MQ=255)
                                                                                                             gcgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGcgtcg                <  1:1637435/139‑1 (MQ=255)
                                                                                                              cgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAgcg                                                                                      >  1:982578/1‑68 (MQ=255)
                                                                                                              cgcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCACCCGACCAgcg                                                                                      <  2:982578/68‑1 (MQ=255)
                                                                                                               gcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAg                                        <  2:4810/113‑1 (MQ=255)
                                                                                                               gcATCTCGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAg                                        >  1:4810/1‑113 (MQ=255)
                                                                                                                  tctcGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCgttgtt           <  1:1437851/139‑1 (MQ=255)
                                                                                                                    tcGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGa                                                                                            >  1:1610149/1‑56 (MQ=255)
                                                                                                                    tcGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGa                                                                                            <  2:1610149/56‑1 (MQ=255)
                                                                                                                     cGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTa                                                                          >  2:1412096/1‑73 (MQ=255)
                                                                                                                     cGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTa                                                                          <  1:1412096/73‑1 (MQ=255)
                                                                                                                     cGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca                                <  2:350600/115‑1 (MQ=255)
                                                                                                                     cGGGTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAacgca                                >  1:350600/1‑115 (MQ=255)
                                                                                                                      gggTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                                                   <  2:1496737/47‑1 (MQ=255)
                                                                                                                      gggTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAggg                                                                                                   >  1:1496737/1‑47 (MQ=255)
                                                                                                                        gTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAcc                                                                        <  1:1438520/72‑1 (MQ=255)
                                                                                                                        gTGCGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAcc                                                                        >  2:1438520/1‑72 (MQ=255)
                                                                                                                           cGAGTTAGATGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAgg  >  2:258468/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                 
TCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGG  >  NC_000913/2637163‑2637424

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