Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,971,761 T→C 100% T311A (ACC→GCC)  tar ← methyl‑accepting chemotaxis protein II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,971,7610TC97.7% 275.3 / ‑3.1 84T311A (ACC→GCC) tarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (39/43);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

TCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGT  >  NC_000913/1971664‑1971857
                                                                                                 |                                                                                                
tCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGa                                                                                                >  1:630550/1‑100 (MQ=255)
 ccaccaCTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAg                                                                                               <  1:2968670/100‑1 (MQ=255)
 ccaccaCTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAg                                                                                               <  1:2007186/100‑1 (MQ=255)
  caccacTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAgc                                                                                              <  1:1444869/100‑1 (MQ=255)
   accacTTTGCCGCCGTGCTTGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAgc                                                                                              <  1:747971/99‑1 (MQ=255)
     cacTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGAGGCGAgctgc                                                                                           >  1:2412388/1‑100 (MQ=255)
     cacTTTGCCGCCGGGCCGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAgctt                                                                                            <  2:524719/99‑2 (MQ=255)
          tGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCAt                                                                                      <  2:1834529/100‑1 (MQ=255)
           gccgccGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATg                                                                                     >  1:744749/1‑100 (MQ=255)
            ccgccgTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGc                                                                                    <  1:2407463/100‑1 (MQ=255)
              gccgTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTg                                                                                  >  1:518094/1‑100 (MQ=255)
                cgTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTggc                                                                                >  2:2971505/1‑100 (MQ=255)
                cgTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTggc                                                                                <  2:1576449/100‑1 (MQ=255)
                 gTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGTGGCGAGCTGCTCCATGCTggcn                                                                               <  2:1784182/100‑2 (MQ=255)
                  tGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTggcgg                                                                              >  1:405770/1‑100 (MQ=255)
                   gCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTggcggc                                                                             <  1:261803/100‑1 (MQ=255)
                    cTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTggcgg                                                                              <  2:943269/98‑1 (MQ=255)
                    cTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCa                                                                            >  2:1330444/1‑100 (MQ=255)
                     tGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAg                                                                           >  1:511696/1‑100 (MQ=255)
                      gggCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGt                                                                          <  1:355948/100‑1 (MQ=255)
                       ggCGGTGTCGGAGGCACTTTTCGCCAGTTGCGAGTTCTGGTTGTTGTTATTGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGtt                                                                         <  1:280225/100‑1 (MQ=37)
                        gCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGttt                                                                        >  2:2719784/1‑100 (MQ=255)
                         cGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTttc                                                                       <  1:2100636/100‑1 (MQ=255)
                          ggTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTttct                                                                      >  1:1625506/1‑100 (MQ=255)
                          ggTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTttct                                                                      <  2:535822/100‑1 (MQ=255)
                           gtgtCGGTGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTttctt                                                                     <  2:2852710/100‑1 (MQ=255)
                             gtCGGAGGCACTTTGCGGCAGTTGCGAGGCCTGGCGGGCGTTTTCGGCGTTTTGCTTCACTGTCGCTGCGAGCTGCTCCATGCTGGCGGCAGTttcttc                                                                    <  2:1971863/99‑1 (MQ=37)
                                gggggCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAg                                                                 <  2:2456327/96‑1 (MQ=255)
                                 gAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAgcg                                                               >  1:3027472/1‑100 (MQ=255)
                                  aGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAgcgc                                                              >  1:1353486/1‑100 (MQ=255)
                                   ggCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGcgcg                                                             >  1:296951/1‑100 (MQ=255)
                                    gCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCgg                                                            <  2:3047252/100‑1 (MQ=255)
                                     cACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGa                                                           >  1:1059647/1‑100 (MQ=255)
                                     cACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGa                                                           >  1:1902198/1‑100 (MQ=255)
                                       cTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATg                                                         <  2:1738745/100‑1 (MQ=255)
                                         ttGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGcc                                                       <  1:468398/100‑1 (MQ=255)
                                          tGCTCCAGTTGCGAGGCCTGGCGTGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGcc                                                       <  1:941900/99‑1 (MQ=255)
                                            cgcCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCCGCGCGGATGCctgc                                                    >  2:2052650/1‑100 (MQ=255)
                                            cgcCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCctgc                                                    <  2:982951/100‑1 (MQ=255)
                                            cgcCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCctgc                                                    <  2:1982843/100‑1 (MQ=255)
                                            cgcCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCctgc                                                    >  2:1094038/1‑100 (MQ=255)
                                            cgcCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCctgc                                                    >  2:2790231/1‑100 (MQ=255)
                                            cgcCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCctgc                                                    >  1:2281242/1‑100 (MQ=255)
                                             gcCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCctgct                                                   >  1:1508213/1‑100 (MQ=255)
                                             gcCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTAGCGGCAGTTTCTTCCAGCGCGGATGCctgct                                                   >  1:1985852/1‑100 (MQ=255)
                                              ccAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCctgctg                                                  <  1:236761/100‑1 (MQ=255)
                                                aGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGt                                                 <  1:1439082/99‑1 (MQ=255)
                                                 gTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTc                                               <  2:270542/100‑1 (MQ=255)
                                                 gTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTc                                               <  1:1629973/100‑1 (MQ=255)
                                                   tGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAg                                             >  2:3031497/1‑100 (MQ=255)
                                                    gCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGt                                            >  2:1829445/1‑100 (MQ=255)
                                                    gCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGt                                            >  1:2726062/1‑100 (MQ=255)
                                                         gCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACggg                                       >  1:1088064/1‑100 (MQ=255)
                                                         gCCTGGCGGGCGTTATCGGCGTTGTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACggg                                       >  2:2851856/1‑100 (MQ=255)
                                                          ccTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGgga                                      <  1:1297147/100‑1 (MQ=255)
                                                           cTGGCGGTCGTTATCGGCGTGTTGCTTCACTTTCGCGGCGAGCTGCTCCATGCTGGCGGTAGGTTCTTCCATCGCGGATGCCTGCTGTTCAGTACGggag                                     <  2:1633560/100‑1 (MQ=255)
                                                             ggcgggcgTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGggagga                                   >  1:1898826/1‑100 (MQ=255)
                                                              gcgggcgTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGaa                                  <  1:3026173/100‑1 (MQ=255)
                                                                  gcgTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGAt                              >  2:162576/1‑100 (MQ=255)
                                                                     ttATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCCGCGCGGATGCCTGCTGTTCAGTCCGGGAGGAAAGATc                             >  2:35450/1‑98 (MQ=39)
                                                                     ttATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCgg                           >  1:1578479/1‑100 (MQ=255)
                                                                       aTCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGgtg                         >  1:1896860/1‑100 (MQ=255)
                                                                       aTCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGgtg                         <  1:3031530/100‑1 (MQ=255)
                                                                       aTCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGgtg                         <  2:2249582/100‑1 (MQ=255)
                                                                          ggCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTg                      >  1:582881/1‑100 (MQ=255)
                                                                          ggCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCAGATGCCTGCTGTACAGTACGGGAGGAAAGATCGGTGTTg                      >  1:2956013/1‑100 (MQ=255)
                                                                           gCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGc                     >  1:2295189/1‑100 (MQ=255)
                                                                            cGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGcn                    <  2:1791503/100‑2 (MQ=255)
                                                                            cGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGcc                    <  1:2467003/100‑1 (MQ=255)
                                                                            cGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGcc                    <  1:653197/100‑1 (MQ=255)
                                                                              ttttGCTTCATGGTCGCGGCGAGCTGCTCCATGCTGGCGTCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGcc                    <  2:2789139/98‑1 (MQ=255)
                                                                                 tGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCccgccg               <  1:2719816/100‑1 (MQ=255)
                                                                                  gCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCcgccgc              >  1:1349119/1‑100 (MQ=255)
                                                                                     tCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAAt           >  2:531626/1‑100 (MQ=255)
                                                                                      cACTGTCGCGGCGAGCTGCTCCATGCTGGTGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTGCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAAtt          <  1:440860/100‑1 (MQ=255)
                                                                                      cACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAAtt          <  1:2942548/100‑1 (MQ=255)
                                                                                      cACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAAtt          <  1:1532286/100‑1 (MQ=255)
                                                                                       aCTGTCGCGGCGAGCTGCTTCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAAttt         <  1:1040959/100‑1 (MQ=255)
                                                                                         tGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAAATCGGTGTTGCCCGCCGCAATTTCa       >  1:272864/1‑100 (MQ=255)
                                                                                          gTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCAc      >  2:262779/1‑100 (MQ=255)
                                                                                          gTCGCGGCGAGCTGCTCCATGCTGGCGGCAGGTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTAGCCCGCCGAAATTTCa       >  2:1854834/1‑99 (MQ=255)
                                                                                           tCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACg     <  2:2816689/100‑1 (MQ=255)
                                                                                           tCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACg     <  2:2816672/100‑1 (MQ=255)
                                                                                             gcggcgAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACgg    >  1:2172928/1‑99 (MQ=255)
                                                                                             gcggcgAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACgg    <  1:934648/99‑1 (MQ=255)
                                                                                              cgGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCaagatcggaaga  >  1:927475/1‑90 (MQ=255)
                                                                                                 |                                                                                                
TCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGT  >  NC_000913/1971664‑1971857

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