Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,971,761 T→C 16.0% 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,7610TC16.0% 173.6 / 25.9 81T311A (ACC→GCC) tarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base C (9/4);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

ATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTAC  >  NC_000913/1971663‑1971859
                                                                                                  |                                                                                                  
aTCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGTCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTg                                                                                                   <  2:710613/100‑1 (MQ=255)
aTCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTg                                                                                                   <  1:2312538/100‑1 (MQ=255)
  ccaccaCTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAg                                                                                                 <  2:1989048/100‑1 (MQ=255)
   caccacTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAgc                                                                                                >  2:1184550/1‑100 (MQ=255)
   caccacTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAgc                                                                                                <  1:2370953/100‑1 (MQ=255)
    accacTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAgct                                                                                               <  1:776405/100‑1 (MQ=255)
    accacTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAgct                                                                                               >  1:1178006/1‑100 (MQ=255)
       acTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAgctgct                                                                                            <  1:1351597/100‑1 (MQ=255)
        cTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTc                                                                                           <  1:1042600/100‑1 (MQ=255)
            gccgccGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCagg                                                                                       <  1:1622503/100‑3 (MQ=255)
             ccgccgTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGc                                                                                      >  1:457721/1‑100 (MQ=255)
              cgccgTGCTGGGCGGTGTCGGAGGCCGTTTGCGCTAGTTGCGAGGTCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATg                                                                                       <  2:1017519/98‑1 (MQ=255)
                 cgTGCTGGGCGGTGTCGGAGGCACTTTGCGTCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGc                                                                                      <  2:675982/96‑1 (MQ=255)
                 cgTGCTGGGCGGTGACGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGTCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTggc                                                                                  <  2:1633993/100‑1 (MQ=255)
                   tggtGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGGGGCCTGGCGGGGGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTggcgg                                                                                <  2:2189240/97‑1 (MQ=255)
                   tGCTGGGCGGTGTCGGTGGCACTTTGCGCCAGTTGGGAGGCGTGGCGGGCGTTATCTGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGc                                                                                      <  2:72390/94‑1 (MQ=255)
                   tGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTggcgg                                                                                >  2:530750/1‑100 (MQ=255)
                   tGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTggcgg                                                                                >  1:1603361/1‑100 (MQ=255)
                       gggCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGt                                                                            >  1:1205407/1‑100 (MQ=255)
                        ggCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGtt                                                                           <  2:1128699/100‑1 (MQ=255)
                          cGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTttc                                                                         <  1:1339214/100‑1 (MQ=255)
                          cGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTttc                                                                         <  2:1623042/100‑1 (MQ=255)
                             tgtCGGAGGCACGTTGCGCCAGTTGCGAGGCTTGGCGGGCTTGATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTTCATGCTGGCGGCAGTttctt                                                                       <  2:931997/99‑1 (MQ=38)
                              gtCGGAGGCACTTTGCGCCAGTTGCGAAGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTttctttc                                                                     >  2:900150/1‑98 (MQ=255)
                                cGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAg                                                                   <  2:1758284/100‑1 (MQ=255)
                                 ggAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAgc                                                                  >  2:664752/1‑100 (MQ=255)
                                      cACTTTGCGCCGGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCgg                                                              >  2:1551709/1‑99 (MQ=255)
                                      cACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGa                                                             <  2:1205396/100‑1 (MQ=255)
                                      cACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGa                                                             >  1:2380387/1‑100 (MQ=255)
                                       aCTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCGTCCAGCGCgg                                                              >  2:1959948/1‑98 (MQ=255)
                                           tGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGcc                                                         >  2:501266/1‑99 (MQ=255)
                                             cgcCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCctgc                                                      <  1:208875/100‑1 (MQ=255)
                                              gcCAGTTGCGAGGCCTAACGGGCGTTAACGGCGTATTGATTCACGGTCGCCGTGAGCTGCTCCATGCTGGCGGGAGTTTCTTCCAGCGCGGATTCcggct                                                     >  2:793227/1‑100 (MQ=25)
                                                cAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGt                                                   <  2:478565/100‑1 (MQ=255)
                                                   ttGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCa                                                <  1:1277973/100‑1 (MQ=255)
                                                   ttGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCa                                                <  2:2402664/100‑1 (MQ=255)
                                                    tGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAg                                               >  1:63802/1‑100 (MQ=255)
                                                    tGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAg                                               <  1:223710/100‑1 (MQ=255)
                                                    tGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAg                                               >  2:2335478/1‑100 (MQ=255)
                                                     gCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGt                                              >  1:1642892/1‑100 (MQ=255)
                                                      cGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTa                                             >  1:2218825/1‑100 (MQ=255)
                                                       gAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTAc                                            >  1:1520349/1‑100 (MQ=255)
                                                            cTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGggag                                       <  2:1886425/100‑1 (MQ=255)
                                                             tGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGggagg                                      >  2:2228039/1‑100 (MQ=255)
                                                              ggcgggcgTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCCGCGCGGATGCCTGCTGTTCAGTACGagagga                                     >  1:767730/1‑100 (MQ=255)
                                                                 gggcgTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAg                                  <  1:95650/100‑1 (MQ=255)
                                                                 gggcgTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAg                                  >  2:2380748/1‑100 (MQ=255)
                                                                 gggcgTTATCGGAGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAg                                  >  1:1576602/1‑100 (MQ=255)
                                                                  ggcgTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGa                                 >  2:1359688/1‑100 (MQ=255)
                                                                  ggcgTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGa                                 <  1:1623953/100‑1 (MQ=255)
                                                                   gcgTTATCGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGAt                                >  1:153734/1‑100 (MQ=255)
                                                                    cgTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATc                               <  2:542207/100‑1 (MQ=255)
                                                                       tATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGgt                            <  2:359644/100‑1 (MQ=255)
                                                                        aTCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGgtg                           >  2:1600681/1‑100 (MQ=255)
                                                                        aTCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGgtg                           >  2:689330/1‑100 (MQ=255)
                                                                         tCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGgtgt                          <  2:1265303/100‑1 (MQ=255)
                                                                         tCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGgtgt                          >  1:1355893/1‑100 (MQ=255)
                                                                         tCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGGTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGgtgt                          <  2:1666554/100‑1 (MQ=255)
                                                                         tCGGCGTTTTGATTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGCTGCCTGCTGTTCAGTACGGGAGGAAAGATCGgtg                           <  1:434928/99‑1 (MQ=255)
                                                                          cGGCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGtt                         >  2:2088378/1‑100 (MQ=255)
                                                                           ggCGTTTTGCTTCACTGTCGCGGTGAGGTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGtt                         <  2:1662663/99‑1 (MQ=255)
                                                                           ggCGTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTg                        <  2:12387/100‑1 (MQ=255)
                                                                             cGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGcc                      >  1:1382303/1‑100 (MQ=255)
                                                                              gTTTTGCTTCACTGTCGCGGCGAGCTGCTCCATGTTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATAGGTGTTGccc                     >  1:140866/1‑100 (MQ=255)
                                                                                tttGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGCTGccccc                   >  2:1148853/1‑98 (MQ=255)
                                                                                 ttGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCccgcc                  >  1:2354337/1‑100 (MQ=255)
                                                                                 ttGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCccgcc                  <  2:248085/100‑1 (MQ=255)
                                                                                  tGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCccgccg                 <  2:50369/100‑1 (MQ=255)
                                                                                   gCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCcgccgc                >  1:2331812/1‑100 (MQ=255)
                                                                                   gCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCcgccgc                >  2:1775686/1‑100 (MQ=255)
                                                                                    cTTCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGGTGCCTGCTGTTCAGTACGGGGGGAAAGATCGGTGTTGCCCGCCGCa               >  1:1827201/1‑100 (MQ=255)
                                                                                      tCACTGTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAAt             <  1:1107690/100‑1 (MQ=255)
                                                                                        aCTTTCGCGCTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAAttt           <  2:932472/100‑1 (MQ=255)
                                                                                           gTCGCGGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCAc        >  1:2212994/1‑100 (MQ=255)
                                                                                             cgcgGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACgg      >  2:2044980/1‑100 (MQ=255)
                                                                                             cgcgGCGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACgg      >  2:1371410/1‑100 (MQ=255)
                                                                                              gcgGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACggg     >  2:1247214/1‑100 (MQ=255)
                                                                                              gcgGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACggg     <  1:2326257/100‑1 (MQ=255)
                                                                                                 gTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTAc  >  1:2115638/1‑100 (MQ=255)
                                                                                                 gTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTAc  >  1:1530334/1‑100 (MQ=255)
                                                                                                 gTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTAc  >  1:1329372/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
ATCCACCACTTTGCCGCCGTGCTGGGCGGTGTCGGAGGCACTTTGCGCCAGTTGCGAGGCCTGGCGGGCGTTATCGGCGTTTTGCTTCACTGTCGCGGTGAGCTGCTCCATGCTGGCGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTAC  >  NC_000913/1971663‑1971859

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: