Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,270,013 A→G 5.5% A200A (GCA→GCG aphA → acid phosphatase/phosphotransferase, class B, non‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,270,0130AG 5.5% 375.0 / 8.4 108A200A (GCA→GCGaphAacid phosphatase/phosphotransferase, class B, non‑specific
Reads supporting (aligned to +/- strand):  ref base A (52/50);  new base G (3/3);  total (55/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

Only 100 of 114 total aligned reads displayed.

TTTCATATTCCTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTC  >  NC_000913/4269876‑4270150
                                                                                                                                         |                                                                                                                                         
tttCATATTCCTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCgg                                                                                                                                <  4:1096166/149‑1 (MQ=255)
      atTCCTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCg                                                                                                                          >  1:555594/1‑149 (MQ=255)
          cTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGt                                                                                                                      >  2:835961/1‑149 (MQ=255)
              caccaACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATcc                                                                                                                  >  3:125260/1‑149 (MQ=255)
                ccaACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGc                                                                                                                >  3:291791/1‑149 (MQ=255)
                   aCATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGTTAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCAtt                                                                                                             >  4:1309655/1‑149 (MQ=255)
                   aCATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGACg                                                                                                                                 >  4:1152512/1‑129 (MQ=255)
                     aTGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCt                                                                                                           >  3:1351471/1‑149 (MQ=255)
                             ggTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCt                                                                                                   <  2:555594/149‑1 (MQ=255)
                              gTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGa                                                                                                                                      >  1:1300585/1‑113 (MQ=255)
                              gTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGa                                                                                                                                      <  2:1300585/113‑1 (MQ=255)
                              gTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTc                                                                                                  >  4:752909/1‑149 (MQ=255)
                                 aTCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCaa                                                                                               >  4:1745747/1‑149 (MQ=255)
                                  tCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAAc                                                                                              >  4:1634050/1‑149 (MQ=255)
                                    tttGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTg                                                                                                                        >  2:482349/1‑121 (MQ=255)
                                    tttGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTg                                                                                                                        <  1:482349/121‑1 (MQ=255)
                                           gCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCAcgc                                                                                                                                        <  1:20837/98‑1 (MQ=255)
                                           gCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCAcgc                                                                                                                                        >  2:20837/1‑98 (MQ=255)
                                                 aaCCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAAccc                                                                               >  4:1686558/1‑149 (MQ=255)
                                                   ccAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCtt                                                                             >  3:694233/1‑149 (MQ=255)
                                                   ccAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCtt                                                                             >  2:717377/1‑149 (MQ=255)
                                                    cAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTg                                                                            <  1:717377/149‑1 (MQ=255)
                                                     aGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTctacc                                                                                        >  4:773915/1‑136 (MQ=255)
                                                          aaaaTACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGgc                                                                                                                               <  4:622100/92‑1 (MQ=255)
                                                          aaaaTACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGgc                                                                                                                               >  3:622100/1‑92 (MQ=255)
                                                          aaaaTACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACaa                                                                      >  1:453411/1‑149 (MQ=255)
                                                              tACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCgg                                                                  <  2:309545/149‑1 (MQ=255)
                                                              tACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCgg                                                                  <  2:453411/149‑1 (MQ=255)
                                                              tACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCgg                                                                  >  2:1841757/1‑149 (MQ=255)
                                                               aCAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCggg                                                                 <  3:1745747/149‑1 (MQ=255)
                                                                cAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGt                                                                >  1:1848546/1‑149 (MQ=255)
                                                                cAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGt                                                                <  4:125260/149‑1 (MQ=255)
                                                                 aaaaTCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTg                                                               >  3:1054471/1‑149 (MQ=255)
                                                                    aTCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCaa                                                                                               <  1:143235/114‑1 (MQ=255)
                                                                    aTCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCaa                                                                                               >  2:143235/1‑114 (MQ=255)
                                                                      cGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGttt                                                          <  4:694233/149‑1 (MQ=255)
                                                                        cAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTgg                                                        >  4:1520070/1‑149 (MQ=255)
                                                                          aTGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGc                                                              <  4:1887176/141‑1 (MQ=255)
                                                                              cTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAaga                                                  <  3:752909/149‑1 (MQ=255)
                                                                               tGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTg                                                                                                          <  1:255688/92‑1 (MQ=255)
                                                                               tGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTg                                                                                                          >  2:255688/1‑92 (MQ=255)
                                                                               tGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGc                                                              >  1:1206875/1‑136 (MQ=255)
                                                                               tGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGc                                                              <  2:1206875/136‑1 (MQ=255)
                                                                                 cAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGt                                               <  3:1520070/149‑1 (MQ=255)
                                                                                   ggATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGa                                             <  4:1351471/149‑1 (MQ=255)
                                                                                    gATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGAt                                            >  4:1630338/1‑149 (MQ=255)
                                                                                    gATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGAt                                            >  3:750499/1‑149 (MQ=255)
                                                                                    gATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGAt                                            <  3:1634050/149‑1 (MQ=255)
                                                                                     aTAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTcc                                                                                                 <  1:1301235/95‑1 (MQ=255)
                                                                                     aTAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATc                                           <  4:291791/149‑1 (MQ=255)
                                                                                     aTAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATc                                           <  3:1686558/149‑1 (MQ=255)
                                                                                       aaaaaTATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTcta                                                                                          <  4:1736167/100‑1 (MQ=255)
                                                                                       aaaaaTATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTcta                                                                                          >  3:1736167/1‑100 (MQ=255)
                                                                                              tCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCa                                  <  4:1054471/149‑1 (MQ=255)
                                                                                                 gAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCTTGGTATCCGCATTCTGCGCGCCTCCAACTctacc                                                                                        >  2:1531654/1‑92 (MQ=255)
                                                                                                 gAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTctacc                                                                                        <  1:1531654/92‑1 (MQ=255)
                                                                                                 gAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAGCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAAt                                     <  4:534404/143‑1 (MQ=255)
                                                                                                 gAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAAt                                     >  3:534404/1‑143 (MQ=255)
                                                                                                 gAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGaa                               >  3:1835347/1‑149 (MQ=255)
                                                                                                 gAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGaa                               >  3:1384919/1‑149 (MQ=255)
                                                                                                 gAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGaa                               <  3:1630338/149‑1 (MQ=255)
                                                                                                    ttttttATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCg                                                                   <  4:1227108/110‑1 (MQ=255)
                                                                                                    ttttttATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCg                                                                   >  3:1227108/1‑110 (MQ=255)
                                                                                                    ttttttATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCa                                  <  2:749458/143‑1 (MQ=255)
                                                                                                    ttttttATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCa                                  >  1:749458/1‑143 (MQ=255)
                                                                                                    ttttttATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATAc                            >  2:1502622/1‑149 (MQ=255)
                                                                                                      ttttATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTg                          <  4:750499/149‑1 (MQ=255)
                                                                                                       tttATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCCCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGa                         <  3:1309655/149‑1 (MQ=255)
                                                                                                            ggCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACagag                    >  1:1367469/1‑149 (MQ=255)
                                                                                                             gCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGc                   >  2:986516/1‑149 (MQ=255)
                                                                                                               gATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCgg                 <  1:1841757/149‑1 (MQ=255)
                                                                                                                 ttCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGga               <  2:1367469/149‑1 (MQ=255)
                                                                                                                   cTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGgaga             <  1:1502622/149‑1 (MQ=255)
                                                                                                                     gATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCgg                 <  3:1693110/143‑1 (MQ=255)
                                                                                                                     gATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCgg                 >  4:1693110/1‑143 (MQ=255)
                                                                                                                        aaTGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGc                                                                           <  2:593026/82‑1 (MQ=255)
                                                                                                                        aaTGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGc                                                                           >  1:593026/1‑82 (MQ=255)
                                                                                                                         aTGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTg                                                         >  1:11703/1‑99 (MQ=255)
                                                                                                                         aTGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTg                                                         <  2:11703/99‑1 (MQ=255)
                                                                                                                         aTGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGa       <  4:1384919/149‑1 (MQ=255)
                                                                                                                          tGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTc                                   <  1:850524/120‑1 (MQ=255)
                                                                                                                          tGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTc                                   >  2:850524/1‑120 (MQ=255)
                                                                                                                          tGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCa                                  >  3:893000/1‑121 (MQ=255)
                                                                                                                          tGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGga               >  1:585302/1‑140 (MQ=255)
                                                                                                                          tGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGAt      <  1:986516/149‑1 (MQ=255)
                                                                                                                           gATATTACCGCCGCGCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCgg                 >  3:1651461/1‑137 (MQ=255)
                                                                                                                           gATATTACCGCCGCGCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCgg                 <  3:1710154/137‑1 (MQ=255)
                                                                                                                           gATATTACCGCCGCGCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCgg                 >  4:1710154/1‑137 (MQ=255)
                                                                                                                           gATATTACCGCCGCGCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCgg                 >  1:1162722/1‑137 (MQ=255)
                                                                                                                           gATATTACCGCCGCGCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCgg                 <  2:1162722/137‑1 (MQ=255)
                                                                                                                             tcttACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTCCCTACAAACCCta                                                                             <  4:1638587/73‑2 (MQ=255)
                                                                                                                              atTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGAt      <  3:420838/145‑1 (MQ=255)
                                                                                                                              atTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGctc  >  2:423533/1‑149 (MQ=255)
                                                                                                                                   cgccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTgcgc                                                                                                       >  2:1443986/1‑43 (MQ=255)
                                                                                                                                         aCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCgg                                                                  >  1:1792485/1‑74 (MQ=255)
                                                                                                                                         aCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCgg                                                                  <  2:1792485/74‑1 (MQ=255)
                                                                                                                                         aCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGaa                               <  4:810357/109‑1 (MQ=255)
                                                                                                                                         aCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGaa                               >  3:810357/1‑109 (MQ=255)
                                                                                                                                         aCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATg     <  1:1254949/135‑1 (MQ=255)
                                                                                                                                         aCGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATg     >  2:1254949/1‑135 (MQ=255)
                                                                                                                                         |                                                                                                                                         
TTTCATATTCCTGCCACCAACATGAATCCGGTGATCTTTGCGGGCGATAAACCAGGGCAAAATACAAAATCGCAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTC  >  NC_000913/4269876‑4270150

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