Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,822,694 G→T 100% G99W (GGG→TGG)  spoT → bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,822,6940GT98.1% 340.7 / ‑3.3 114G99W (GGG→TGG) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base T (58/52);  total (62/52)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

GATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCG  >  NC_000913/3822597‑3822793
                                                                                                 |                                                                                                   
gATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAgtgg                                                                                                   <  1:996267/100‑1 (MQ=255)
 aTGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAgtggg                                                                                                  <  1:1988166/100‑1 (MQ=255)
 aTGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAgtggg                                                                                                  <  1:1275878/100‑1 (MQ=255)
  tGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGgt                                                                                                 >  1:541026/1‑100 (MQ=255)
   ggcggcGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGgtg                                                                                                <  1:668985/100‑1 (MQ=255)
   ggcggcGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGgtg                                                                                                <  1:524383/100‑1 (MQ=255)
   ggcggcGCTGCTGCAGGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGgtg                                                                                                <  1:1499056/100‑1 (MQ=255)
        cgcTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGaa                                                                                           >  2:127167/1‑100 (MQ=255)
        cgcTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTAGTAGAGTGGGTGTCGaa                                                                                           >  2:1432171/1‑100 (MQ=255)
            gctgcATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACtt                                                                                       >  2:133064/1‑100 (MQ=255)
            gctgcATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACtt                                                                                       <  1:299953/100‑1 (MQ=255)
             ctgcATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTg                                                                                      <  1:403346/100‑1 (MQ=255)
              tgcATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGa                                                                                     >  1:1386484/1‑100 (MQ=255)
               gcATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGAt                                                                                    >  2:1326104/1‑100 (MQ=255)
               gcATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGAt                                                                                    <  1:338019/100‑1 (MQ=255)
                 aTGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTTTTTTTTTTTTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGCTa                                                                                   >  2:667492/1‑99 (MQ=255)
                   gACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAAc                                                                                >  2:2486053/1‑100 (MQ=255)
                   gACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAAc                                                                                >  2:354491/1‑100 (MQ=255)
                    ncGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACt                                                                               >  2:1958105/2‑100 (MQ=255)
                    aCGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACt                                                                               <  1:707200/100‑1 (MQ=255)
                     cGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTc                                                                              <  1:2218206/100‑1 (MQ=255)
                     cGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTc                                                                              <  1:2236257/100‑1 (MQ=255)
                     cGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTc                                                                              >  1:233479/1‑100 (MQ=255)
                      gTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCa                                                                             >  2:1071745/1‑100 (MQ=255)
                       tGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTTTTTTTTTTTTCACCGAGATGGTAGAGTGGGTATCTAATCTTGATAAACGCaa                                                                            >  2:1771440/1‑100 (MQ=255)
                       tGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCaa                                                                            <  2:396614/100‑1 (MQ=255)
                       tGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCaa                                                                            >  1:1616591/1‑100 (MQ=255)
                        gATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAg                                                                           >  1:342739/1‑100 (MQ=255)
                            gAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTcc                                                                       >  2:1529099/1‑100 (MQ=255)
                            gAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTcc                                                                       >  1:1502925/1‑100 (MQ=255)
                              aGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCcgc                                                                     >  1:2245457/1‑100 (MQ=255)
                              aGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCcgc                                                                     >  2:2257853/1‑100 (MQ=255)
                              aGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCcgc                                                                     <  2:84835/100‑1 (MQ=255)
                               ggtACTCCCGCCGCCTACCAGGATATGGAACAGCTTGTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCcgc                                                                     <  2:1162766/97‑1 (MQ=255)
                               gATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCcgcg                                                                    <  1:127169/100‑1 (MQ=255)
                                 tACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGAt                                                                  >  2:796933/1‑100 (MQ=255)
                                   cTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGTTGTCGAAACTTGATAAACTCAAGTTCcgcg                                                                    <  2:32555/96‑1 (MQ=255)
                                   cTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATaa                                                                >  1:2483461/1‑100 (MQ=255)
                                   cTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATaa                                                                <  1:450786/100‑1 (MQ=255)
                                    tCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATaag                                                               >  2:2376643/1‑100 (MQ=255)
                                     cccGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATaaga                                                              >  1:2304248/1‑100 (MQ=255)
                                      ccGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATaagaa                                                             >  2:874834/1‑100 (MQ=255)
                                      ccGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATaagaa                                                             >  2:2385843/1‑100 (MQ=255)
                                          cACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCTAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAagag                                                         <  2:1637818/100‑1 (MQ=255)
                                           aCCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGATGGTAGAGTGTGAGTCGAAACGTGATAAACTCAAGTTCCGCGATAAGAAAGAgg                                                        <  2:2071311/100‑1 (MQ=255)
                                            ccTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAagagac                                                       >  2:1556128/1‑98 (MQ=255)
                                             cTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGgcg                                                      >  1:1571547/1‑100 (MQ=255)
                                                ccAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAg                                                   >  2:85850/1‑100 (MQ=255)
                                                    gATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCg                                               <  1:403604/100‑1 (MQ=255)
                                                     atatGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGa                                              >  2:2320703/1‑100 (MQ=255)
                                                     atatGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGa                                              <  1:1393575/100‑1 (MQ=255)
                                                      tatGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGaa                                             >  1:1441215/1‑100 (MQ=255)
                                                      tatGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGaa                                             <  1:575505/100‑1 (MQ=255)
                                                      tatGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGaa                                             >  1:2485720/1‑100 (MQ=255)
                                                      tatGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGCGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGaa                                             >  2:680997/1‑100 (MQ=255)
                                                        tGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGaaaa                                           >  2:1734734/1‑100 (MQ=255)
                                                        tGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGaaaa                                           <  1:1200737/100‑1 (MQ=255)
                                                        tGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGAGATAAGAAAGAGGCGCAGGCCGaaaa                                           >  1:1974990/1‑100 (MQ=255)
                                                         ggAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAAc                                          >  2:1944152/1‑100 (MQ=255)
                                                         ggAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAAc                                          <  2:257905/100‑1 (MQ=255)
                                                         ggAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAAc                                          <  2:562832/100‑1 (MQ=255)
                                                          gAACAACCTTTTGGTAATCGCATCGCCGAGCTGGTAGAGTGCGTGTCGAAACTTTATATTCTCACGTTCCGTCATAATAAAGAGGCGCAGTCCg                                               <  2:1720325/94‑1 (MQ=255)
                                                              aGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTcn                                     <  1:1707600/100‑2 (MQ=255)
                                                              aGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCg                                     <  2:124957/100‑1 (MQ=255)
                                                               gCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCCAGTTGCGCGATAAGAAGGAGCCGCAGGAAGAAAACtt                                        >  2:2528233/1‑96 (MQ=255)
                                                               gCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGc                                    >  1:10689/1‑100 (MQ=255)
                                                                cTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAATTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCa                                   <  2:396770/100‑1 (MQ=255)
                                                                 tttttGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCaa                                  >  1:334594/1‑100 (MQ=255)
                                                                 tatttGGTAAAAGGGTCGCCGAGCTGGTAGAATGGGTGTCGATACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCaa                                  >  1:1730092/3‑100 (MQ=255)
                                                                  ttttGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAg                                 <  1:98990/100‑1 (MQ=255)
                                                                  ttttGGTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAg                                 >  2:1879654/1‑100 (MQ=255)
                                                                  ttttGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAg                                 >  2:922375/1‑100 (MQ=255)
                                                                       gTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAgatgat                            >  2:1883407/1‑100 (MQ=255)
                                                                       gTAAAAGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAgatgat                            <  2:517691/100‑1 (MQ=255)
                                                                         aaaaGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTa                          >  1:2097753/1‑100 (MQ=255)
                                                                          aaaGCGTCGCCGGGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTAt                         >  1:539622/1‑100 (MQ=255)
                                                                          aaaGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTAt                         >  2:1508798/1‑100 (MQ=255)
                                                                            aGCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTGGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATgg                       <  2:1766055/100‑1 (MQ=255)
                                                                             gCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGc                      <  1:1965395/100‑1 (MQ=255)
                                                                             gCGTCGCCGAGCTGGTAGAGTGGGTGTCGAAAATTGATCAACTCCAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGc                      <  2:2059688/100‑1 (MQ=255)
                                                                              cGTCGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCg                     >  1:146834/1‑100 (MQ=255)
                                                                               gTCGCCGAGCTGGCAGAGGGAGGGTCGAAACTTGATAAACTCAAGTTCCGCGGCAAGAGAGATGCGAAGACCAAAAACGTTCACAACATGAACATGGCGa                    >  2:2106766/1‑100 (MQ=255)
                                                                                 cGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGGCGGCCGAAAACTTTCGCAAGATGATTATGGCGATg                  >  1:94989/1‑100 (MQ=255)
                                                                                 cGCCGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATg                  <  1:1167774/100‑1 (MQ=255)
                                                                                   ccGGGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGt                >  1:413953/1‑100 (MQ=255)
                                                                                   ccGAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGt                >  1:1426/1‑100 (MQ=255)
                                                                                     gAGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGc              >  2:847256/1‑100 (MQ=255)
                                                                                      ngCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCa             >  2:1514608/2‑100 (MQ=255)
                                                                                      aGCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCa             <  1:663232/100‑1 (MQ=255)
                                                                                       gCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAg            >  1:1984106/1‑100 (MQ=255)
                                                                                       gCTGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAg            <  2:2403694/100‑1 (MQ=255)
                                                                                         tGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGa          <  1:1071753/100‑1 (MQ=255)
                                                                                         tGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGa          >  1:208297/1‑100 (MQ=255)
                                                                                         tGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGa          <  1:2093942/100‑1 (MQ=255)
                                                                                         tGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGa          <  1:74963/100‑1 (MQ=255)
                                                                                         tGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGa          <  2:2043631/100‑1 (MQ=255)
                                                                                         tGGTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGa          >  2:1846280/1‑100 (MQ=255)
                                                                                          ggTAGAGTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCTGAAAACGTTCGCAAGATGATTATGGCGATGGTGCAGGat         >  2:1800790/1‑100 (MQ=255)
                                                                                             agagTGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATc      <  1:2527826/100‑1 (MQ=255)
                                                                                                 ggggTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCcgcg  <  1:922383/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
GATGGCGGCGCTGCTGCATGACGTGATTGAAGATACTCCCGCCACCTACCAGGATATGGAACAGCTTTTTGGTAAAAGCGTCGCCGAGCTGGTAGAGGGGGTGTCGAAACTTGATAAACTCAAGTTCCGCGATAAGAAAGAGGCGCAGGCCGAAAACTTTCGCAAGATGATTATGGCGATGGTGCAGGATATCCGCG  >  NC_000913/3822597‑3822793

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