Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,684,019 C→T 5.6% 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 5.6% 276.8 / 4.5 89R380W (CGG→TGG) rstBsensory histidine kinase of RstAB two‑component system
Reads supporting (aligned to +/- strand):  ref base C (47/36);  new base T (2/3);  total (50/39)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

CCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCC  >  NC_000913/1683873‑1684164
                                                                                                                                                  |                                                                                                                                                 
cccTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCggg                                                                                                                                                <  2:724096/150‑1 (MQ=255)
    gCGGTACTGCCATTCCACGGGTGAAACCACCCTGCTACTGTCGGGGAATAGAGCGTCATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCg                                                                                                                                            <  2:92766/150‑1 (MQ=255)
     cGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGc                                                                                                                                           <  2:713754/150‑1 (MQ=255)
      ggTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGAGGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCt                                                                                                                                          <  1:1074712/150‑1 (MQ=255)
       gTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgagg                                                                                                                                         <  2:459876/150‑4 (MQ=255)
        tACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCa                                                                                                                                        <  2:214458/150‑1 (MQ=255)
         aCTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATAAATTGATGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggcg                                                                                                                                       <  2:547362/150‑6 (MQ=255)
         aCTGCCATTCAACGGCTGAAACCAGCCTGCTACTGGCGGGGAATAGAGCGACATTAATTGTTGAGGAGGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCGGCTGGGATCgaggcg                                                                                                                                       <  1:24175/150‑6 (MQ=255)
          cTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGGACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAc                                                                                                                                      >  2:105545/1‑150 (MQ=255)
          cTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGCACCTTTTGTTCGCCTCGATCCCAGCCGGGATCCCTCAAc                                                                                                                                      >  2:1702759/1‑150 (MQ=255)
           tGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAAcc                                                                                                                                     <  2:1039158/150‑1 (MQ=255)
           tGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACAGATCTTTGAACCTTTTGTTCGCCTCGATCCCCGCCGGGATCGCTCAAcc                                                                                                                                     >  1:944877/1‑150 (MQ=255)
              cATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGCCATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGACCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACcggc                                                                                                                                  >  1:1696661/1‑150 (MQ=255)
               aTTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCCACcggcg                                                                                                                                 >  2:279621/1‑150 (MQ=255)
                ttCAACGGTTGAAACCAGCCTGCTACTGTCTGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACcggcgg                                                                                                                                <  1:570099/150‑1 (MQ=255)
                  cAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAACCCCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAGCCggcggcg                                                                                                                              >  2:120319/1‑149 (MQ=255)
                  cAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCt                                                                                                                              >  1:1149523/1‑150 (MQ=255)
                  cAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTGTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCggcggcg                                                                                                                              >  1:1188715/1‑149 (MQ=255)
                       gTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCCACCGGCGGCTGCggg                                                                                                                         >  2:1787108/1‑150 (MQ=255)
                        ttGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggc                                                                                                                        >  2:611732/1‑150 (MQ=255)
                         tgcaaCCAGCCTGCTACTGTCGGGGAATAGAGCGCCATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTACGCCTCGATCCCAGCCGGGATCTCTCAACAGGCGGCTGCgggcg                                                                                                                       >  2:1125836/4‑149 (MQ=255)
                          gAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctg                                                                                                                      <  2:1600563/150‑1 (MQ=255)
                          gAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctg                                                                                                                      <  1:1810457/150‑1 (MQ=255)
                           aaaCCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctgg                                                                                                                     >  2:1033893/1‑150 (MQ=255)
                             aCCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGGTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctgggg                                                                                                                   <  2:1543161/150‑1 (MQ=255)
                               cAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggct                                                                                                                 >  2:415304/1‑150 (MQ=255)
                                aGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctggggctg                                                                                                                >  1:1259489/1‑150 (MQ=255)
                                  ccTACTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGGTCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTAGGGCTGGc                                                                                                              <  1:902217/150‑1 (MQ=255)
                                   cTGCTACTGTCGGCGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCGCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCa                                                                                                             <  1:146904/150‑1 (MQ=255)
                                      cTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAAtt                                                                                                          <  2:1227008/150‑1 (MQ=255)
                                       tACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTg                                                                                                         <  2:839571/150‑1 (MQ=255)
                                          ggggggggAATAGAGCGACATTAATTGTTGAGGATGATGGCCGAGGGATTGCCCCTGAATACCGCGAACATAGATTTGAACCTTTTGTTCTCCTCGATCCCAGCCGGGATCGCTCAACCTGCGGCTGCGGGCTTGGGCTGGCAATTGTcc                                                                                                      <  2:544947/146‑1 (MQ=255)
                                            gcGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCAc                                                                                                    <  1:710186/149‑1 (MQ=255)
                                             cGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCAct                                                                                                   <  1:529107/150‑1 (MQ=255)
                                               gggAATAGCGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCCTTTGTTCGCCTCGATCCCAGCAGGGATCGATCAACCGGCGGCTGCGGGCTGGTGCTGGCAATTGTACActct                                                                                                 >  1:1389339/1‑150 (MQ=255)
                                                ggAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCta                                                                                                >  1:478456/1‑150 (MQ=255)
                                                       agCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACt                                                                                         >  2:1848513/1‑150 (MQ=255)
                                                       agCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACt                                                                                         <  1:1733095/150‑1 (MQ=255)
                                                          gaaaTTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGc                                                                                      <  2:1521650/147‑1 (MQ=255)
                                                          gACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTAAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGc                                                                                      <  2:1577414/150‑1 (MQ=255)
                                                             aTTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAAt                                                                                   >  1:477817/1‑150 (MQ=255)
                                                              ttAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATg                                                                                  >  2:578797/1‑150 (MQ=255)
                                                               tAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATgg                                                                                 <  2:1566152/150‑1 (MQ=255)
                                                                   tgttgAGGATGATGGCCCAGGGATTGCCCCAGAAACCCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCgg                                                                             <  1:1459998/150‑1 (MQ=255)
                                                                     ttgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCACCCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCCCTCTATAGCACTGGCAATGGGCGGTa                                                                           >  1:933408/1‑150 (MQ=255)
                                                                     ttgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTa                                                                           >  2:1367296/1‑150 (MQ=255)
                                                                          gatgatGGCCCAGGGATTGCCCCAGAAAACCGCGATCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCCACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGtt                                                                      <  1:1162084/150‑1 (MQ=255)
                                                                             gatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCCGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTACACTCTATAGCACTGGCAATGGGCGGTACGGTTAAt                                                                   >  2:1376542/1‑150 (MQ=255)
                                                                              atgGCCCAGGGATAGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGCTCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACTGTTAAtt                                                                  <  2:422278/150‑1 (MQ=255)
                                                                                 gCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATATCACTGGCAATGGGCGGTACGGTTAATtgtg                                                               <  2:202903/150‑1 (MQ=255)
                                                                                      gggATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACAcc                                                          >  2:608420/1‑150 (MQ=255)
                                                                                            gCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGCCACCAGCGaa                                                    >  1:382881/1‑150 (MQ=255)
                                                                                               ccAGAAAACCGCGAACATATCTTTGAAACTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCACCCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTTATTGTGACACCAGCTACCTg                                                 >  1:806484/1‑150 (MQ=255)
                                                                                                 aGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGTCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTggg                                               <  1:850265/150‑1 (MQ=255)
                                                                                                  gAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGggt                                              >  1:1820867/1‑150 (MQ=255)
                                                                                                       ccGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCACCCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTACTTGTGCCACCAGCGAGCTGGGTGGTGc                                         >  2:1796376/1‑150 (MQ=255)
                                                                                                         gcgAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGccc                                       >  2:973831/1‑150 (MQ=255)
                                                                                                         gcgAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGccc                                       >  2:565608/1‑150 (MQ=255)
                                                                                                           gAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCAACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGc                                     >  1:1337209/1‑150 (MQ=255)
                                                                                                             aCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGGTCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCtt                                   >  2:1196347/1‑150 (MQ=255)
                                                                                                              cATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGAGCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTtc                                  <  2:944877/150‑1 (MQ=255)
                                                                                                                 atCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCg                               >  1:399498/1‑150 (MQ=255)
                                                                                                                 atCTTTGAACCTTTTGTACGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCg                               <  2:253809/150‑1 (MQ=255)
                                                                                                                  tcgatGAGCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGATCAACCGTCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGt                              >  2:1003771/5‑150 (MQ=255)
                                                                                                                   cTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGtt                             >  2:946044/1‑150 (MQ=255)
                                                                                                                    tgtGAACCATTTGCTCGCCGCGATCCCAGCCCGGAGCGCTCAACCGGCGGCTGCGGGCGGGTGCTGGCAATTGTCCACTCTATCGCACTGGCAATGGGCGGTACGTTTAATTGTGACACCAGCGAACTGGGTGGGGCCCGCTTCTCGttt                            <  1:1774512/148‑1 (MQ=255)
                                                                                                                     ttGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTa                           <  1:41999/150‑1 (MQ=255)
                                                                                                                       gAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGc                         >  2:236698/1‑150 (MQ=255)
                                                                                                                        aaCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCt                        >  2:1576174/1‑150 (MQ=255)
                                                                                                                        aaCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCt                        <  2:1089536/150‑1 (MQ=255)
                                                                                                                        aaCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCGGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCt                        >  2:203099/1‑150 (MQ=255)
                                                                                                                         aCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTg                       >  1:820328/1‑150 (MQ=255)
                                                                                                                            ttttGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGcc                    <  1:1565358/150‑1 (MQ=255)
                                                                                                                             tttGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCg                   <  1:415304/150‑1 (MQ=255)
                                                                                                                                 ttCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTAt               <  2:477817/150‑1 (MQ=255)
                                                                                                                                 ttCGCCGCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATAGTCCACTCTATAGCACTGCCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTAt               <  2:1149523/150‑1 (MQ=255)
                                                                                                                                   cGCCTCGATCCCCGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACCCCAGCGAACTGGGTGGTGCCCGCGGCGCGGTTAGCTGGCCGTTTTgg             >  2:891262/1‑150 (MQ=255)
                                                                                                                                    gCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGc            >  1:930362/1‑150 (MQ=255)
                                                                                                                                    gCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGc            >  1:858174/1‑150 (MQ=255)
                                                                                                                                    gCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGc            >  1:826884/1‑150 (MQ=255)
                                                                                                                                     gagcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACCCCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCa           >  1:448875/4‑150 (MQ=255)
                                                                                                                                      agcGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGTCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTGTAGATGGCCGTTATGGCAt          >  1:697086/3‑150 (MQ=255)
                                                                                                                                         gATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAAc       >  2:357419/1‑150 (MQ=255)
                                                                                                                                         gATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAAc       <  1:1848604/150‑1 (MQ=255)
                                                                                                                                         gATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAAc       >  1:1630718/1‑150 (MQ=255)
                                                                                                                                         gATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCAGCGTCTCGTTTAGCTGGACGTTATGGCATAAc       >  2:138718/1‑150 (MQ=255)
                                                                                                                                          aTCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACa      >  2:733494/1‑150 (MQ=255)
                                                                                                                                           tCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACAt     >  1:220539/1‑150 (MQ=255)
                                                                                                                                              cAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTGAGCTGGCCGTTATGGCAGGACATccc  >  2:83716/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCC  >  NC_000913/1683873‑1684164

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: