Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,684,019 C→T 8.4% R380W (CGG→TGG)  rstB → sensory histidine kinase of RstAB two‑component system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,0190CT 8.4% 334.5 / 16.2 107R380W (CGG→TGG) rstBsensory histidine kinase of RstAB two‑component system
Reads supporting (aligned to +/- strand):  ref base C (64/34);  new base T (3/6);  total (67/40)
Fisher's exact test for biased strand distribution p-value = 7.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-01

Only 100 of 107 total aligned reads displayed.

GGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCC  >  NC_000913/1683879‑1684164
                                                                                                                                            |                                                                                                                                                 
ggTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCg                                                                                                                                            <  2:920439/148‑1 (MQ=255)
ggTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGc                                                                                                                                           <  3:234014/149‑1 (MQ=255)
 gTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgag                                                                                                                                          <  4:110015/149‑3 (MQ=255)
 gTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgag                                                                                                                                          <  4:626682/149‑3 (MQ=255)
   aCTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgagg                                                                                                                                         <  2:462906/148‑4 (MQ=255)
     tGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAc                                                                                                                                      >  3:101098/1‑149 (MQ=255)
      gCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgagg                                                                                                                                         >  2:563854/1‑142 (MQ=255)
        cATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACcg                                                                                                                                    >  1:247327/1‑148 (MQ=255)
          ttCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACcggcg                                                                                                                                 >  4:683997/1‑149 (MQ=255)
               cGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTg                                                                                                                             >  2:769937/1‑148 (MQ=255)
                   tGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCggg                                                                                                                         >  1:374831/1‑148 (MQ=255)
                    gAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTGATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggct                                                                                                                       >  4:787054/1‑149 (MQ=255)
                       aCCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctgg                                                                                                                     >  2:850420/1‑148 (MQ=255)
                         cAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggc                                                                                                                  >  3:475351/1‑149 (MQ=255)
                           gCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggctg                                                                                                                >  3:97471/1‑149 (MQ=255)
                               gCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCa                                                                                                             <  1:443515/148‑1 (MQ=255)
                                 tACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAAt                                                                                                           <  2:247327/148‑1 (MQ=255)
                                            aaTAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCg                                                                                                                                                  >  1:50006/1‑98 (MQ=255)
                                            aaTAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCg                                                                                                                                                  <  2:50006/98‑1 (MQ=255)
                                               agagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAg                                                                                             >  1:648512/1‑148 (MQ=255)
                                                gagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGAt                                                                                                                                              >  4:557720/1‑98 (MQ=255)
                                                gagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATACCAGCCGGGAt                                                                                                                                              <  3:557720/98‑1 (MQ=255)
                                                  gCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACt                                                                                         <  3:755199/149‑1 (MQ=255)
                                                     aCATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGc                                                                                      >  3:594602/1‑149 (MQ=255)
                                                       aTTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCaa                                                                                    >  3:565794/1‑149 (MQ=255)
                                                        ttAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCaa                                                                                    >  1:287383/1‑148 (MQ=255)
                                                        ttAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCCCTCAACCGGCGGCTGCGGGCTGGGGCTGGAAAGTGTCCACTCTATAGCACTGGCaa                                                                                    >  2:404944/1‑148 (MQ=255)
                                                          aaTTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATgg                                                                                 >  3:720163/1‑149 (MQ=255)
                                                           aTTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGAt                                                                                                                                              >  4:610704/1‑87 (MQ=255)
                                                           aTTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGAt                                                                                                                                              <  3:610704/87‑1 (MQ=255)
                                                               ttgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGt                                                                            >  4:948562/1‑149 (MQ=255)
                                                                tgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTa                                                                           <  3:236420/149‑1 (MQ=255)
                                                                 gAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCCCTCTATAGCACTGGCAATGGGCGGTa                                                                           <  2:287383/148‑1 (MQ=255)
                                                                   ggATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACg                                                                         <  1:850420/148‑1 (MQ=255)
                                                                    gatgatGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGt                                                                       <  3:386549/149‑1 (MQ=255)
                                                                     atgatgGCCCAGGGATTGCCCCATAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCg                                                                                                                                                  <  2:289301/73‑1 (MQ=255)
                                                                     atgatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCg                                                                                                                                                  >  1:289301/1‑73 (MQ=255)
                                                                     atgatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGt                                                                       >  2:905001/1‑148 (MQ=255)
                                                                       gatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTaa                                                                    <  4:210318/149‑1 (MQ=255)
                                                                       gatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTaa                                                                    >  3:557872/1‑149 (MQ=255)
                                                                           gCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATtg                                                                 <  2:545613/148‑1 (MQ=255)
                                                                            cccAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATtgt                                                                <  1:48987/148‑1 (MQ=255)
                                                                             ccAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTACTTGTGa                                                              >  3:637709/1‑149 (MQ=255)
                                                                             ccAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGGTAATTGTGa                                                              >  4:875786/1‑149 (MQ=255)
                                                                              cAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGa                                                              >  1:307671/1‑148 (MQ=255)
                                                                                gggATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGacac                                                           >  3:720052/1‑149 (MQ=255)
                                                                                        cccAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGaa                                                    >  2:146985/1‑148 (MQ=255)
                                                                                        cccAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGaa                                                    >  1:740747/1‑148 (MQ=255)
                                                                                         ccAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAAc                                                   >  1:72133/1‑148 (MQ=255)
                                                                                          cAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACt                                                  >  1:150072/1‑148 (MQ=255)
                                                                                            gAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGggctggggctgg                                                                                                               <  3:521319/85‑1 (MQ=255)
                                                                                            gAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGggctggggctgg                                                                                                               >  4:521319/1‑85 (MQ=255)
                                                                                             aaaaCCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGggt                                              >  4:349687/1‑149 (MQ=255)
                                                                                                aCCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGggtgg                                            >  2:964412/1‑148 (MQ=255)
                                                                                                    cgAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGcc                                        >  1:444594/1‑148 (MQ=255)
                                                                                                      aaCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCg                                      >  2:100405/1‑148 (MQ=255)
                                                                                                      aaCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGACAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCg                                      >  2:631638/1‑148 (MQ=255)
                                                                                                         atatCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTtc                                  >  3:130869/1‑149 (MQ=255)
                                                                                                         atatCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTtc                                  >  4:608546/1‑149 (MQ=255)
                                                                                                           atCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTtct                                 <  2:307671/148‑1 (MQ=255)
                                                                                                             cTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCg                               >  2:748682/1‑148 (MQ=255)
                                                                                                             cTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCg                               >  1:595497/1‑148 (MQ=255)
                                                                                                             cTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGt                              <  4:594602/149‑1 (MQ=255)
                                                                                                               ttGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGTCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGtt                             >  1:720250/1‑148 (MQ=255)
                                                                                                               ttGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGtt                             >  1:149299/1‑148 (MQ=255)
                                                                                                               ttGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGtt                             <  2:150072/148‑1 (MQ=255)
                                                                                                                tGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGttt                            <  1:964412/148‑1 (MQ=255)
                                                                                                                tGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGttt                            <  2:648512/148‑1 (MQ=255)
                                                                                                                   aCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCt                        <  3:20700/149‑1 (MQ=255)
                                                                                                                    ccTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTg                       >  3:891418/1‑149 (MQ=255)
                                                                                                                     cTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTgg                      >  3:723259/1‑149 (MQ=255)
                                                                                                                     cTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTg                       >  1:118837/1‑148 (MQ=255)
                                                                                                                      ttttGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTgg                      >  1:487545/1‑148 (MQ=255)
                                                                                                                      ttttGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGc                     >  3:389751/1‑149 (MQ=255)
                                                                                                                         tGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGt                  >  3:316840/1‑149 (MQ=255)
                                                                                                                           ttCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTa                >  3:457459/1‑149 (MQ=255)
                                                                                                                            tCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTa                <  1:84377/148‑1 (MQ=255)
                                                                                                                             cGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTAt               <  1:146985/148‑1 (MQ=255)
                                                                                                                             cGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTAt               <  2:72133/148‑1 (MQ=255)
                                                                                                                             cGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATg              <  4:316840/149‑1 (MQ=255)
                                                                                                                              gCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATgg             >  4:139496/1‑149 (MQ=255)
                                                                                                                               ccTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAAt                                                                   <  1:305677/94‑1 (MQ=255)
                                                                                                                               ccTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAAt                                                                   >  2:305677/1‑94 (MQ=255)
                                                                                                                                cTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATgga            >  2:877900/1‑147 (MQ=255)
                                                                                                                                cTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCa           >  4:459457/1‑149 (MQ=255)
                                                                                                                                agcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATgg             >  4:366750/3‑147 (MQ=255)
                                                                                                                                 tCGANCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCAt          >  3:850405/1‑149 (MQ=255)
                                                                                                                                 gcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATg              <  1:589382/144‑1 (MQ=255)
                                                                                                                                 gcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATg              >  2:589382/2‑145 (MQ=255)
                                                                                                                                    aTCCCAGCCGGGATCGCTCAACCGGCGGCTGCGAGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAAc       >  3:812839/1‑149 (MQ=255)
                                                                                                                                     tCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTg                       <  4:248336/132‑1 (MQ=255)
                                                                                                                                     tCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTg                       >  3:248336/1‑132 (MQ=255)
                                                                                                                                     tCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAAc       <  1:905001/148‑1 (MQ=255)
                                                                                                                                      cccAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACAt     >  4:62265/1‑149 (MQ=255)
                                                                                                                                       ccAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCGGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATc    >  3:942210/1‑149 (MQ=255)
                                                                                                                                       ccAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATc    >  4:445507/1‑149 (MQ=255)
                                                                                                                                       ccAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATc    <  4:648561/149‑1 (MQ=255)
                                                                                                                                         aGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACTCCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATcc   <  2:194107/148‑1 (MQ=255)
                                                                                                                                         aGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATccc  >  4:193558/1‑149 (MQ=255)
                                                                                                                                         aGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATcc   <  2:374831/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
GGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCC  >  NC_000913/1683879‑1684164

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