Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,821,953 T→A 6.1% L173Q (CTG→CAG)  metW → methionine biosynthesis protein MetW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,821,9530TA 6.1% 249.4 / 6.6 84L173Q (CTG→CAG) metWmethionine biosynthesis protein MetW
Reads supporting (aligned to +/- strand):  ref base T (45/33);  new base A (3/2);  total (48/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGC  >  NC_002947/5821813‑5822093
                                                                                                                                              |                                                                                                                                            
gTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa                                                                                                                                     <  1:698864/150‑1 (MQ=255)
     tGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACt                                                                                                                                >  2:1175619/1‑150 (MQ=255)
     tGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACt                                                                                                                                >  1:454806/1‑150 (MQ=255)
     tGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACt                                                                                                                                >  1:454650/1‑150 (MQ=255)
      ggCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACtt                                                                                                                               <  2:982379/150‑1 (MQ=255)
      ggCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCAACTTCGCCGTATAAGAAGAACCTTGTCCCGAACTCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACtt                                                                                                                               <  1:515426/150‑1 (MQ=255)
       gCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTg                                                                                                                              >  2:652012/1‑150 (MQ=255)
         gACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACATCACTCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCa                                                                                                                            >  1:1079417/1‑150 (MQ=255)
           ccAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCAcc                                                                                                                          >  2:179950/1‑150 (MQ=255)
            cAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATTCCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCg                                                                                                                         >  1:1094974/1‑150 (MQ=255)
                  ccGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgg                                                                                                                   >  1:378993/1‑150 (MQ=255)
                   cGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCgg                                                                                                                                         <  2:1050022/127‑1 (MQ=255)
                   cGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCAAAGGTACTCGACCGCCTGGCgg                                                                                                                                         >  1:1050022/1‑127 (MQ=255)
                     cATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtg                                                                                                                >  2:1487283/1‑150 (MQ=255)
                       tGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtggg                                                                                                              <  1:652012/150‑1 (MQ=255)
                           ggTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGggtggggtgg                                                                                                          >  1:1228260/1‑150 (MQ=255)
                           ggTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGggtggggtgg                                                                                                          <  2:454806/150‑1 (MQ=255)
                           ggTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGggtggggtgg                                                                                                          <  2:454650/150‑1 (MQ=255)
                              aTCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCg                                                                                                       >  1:624056/1‑150 (MQ=255)
                                 ggACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGGCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCt                                                                                                    >  1:550068/1‑150 (MQ=255)
                                  gACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTa                                                                                                   >  1:1314633/1‑150 (MQ=255)
                                   aCTTCATGCCGTATGCCTTGTACAACACGCCGAACATCCACTTCTGCACCTTAGCAAACGCAGAAAGTATGTGGCACACACCCCGGGCCAAGGTACTCGACCGCCTGGAGGTCGCCCACTTGCACCGCAACGGGTGGGGTGGCCGGCTAt                                                                                                  <  1:1572243/150‑1 (MQ=255)
                                      tCATGCCGTATACCTGGTACAACACGCCGAACACCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGc                                                                                               <  2:1094974/150‑1 (MQ=255)
                                       cATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGGCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGcc                                                                                              <  2:550068/150‑1 (MQ=255)
                                           ccGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGc                                                                                                                                           >  2:815663/1‑101 (MQ=255)
                                           ccGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGc                                                                                                                                           <  1:815663/101‑1 (MQ=255)
                                             gTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAAtct                                                                                        >  1:173600/1‑150 (MQ=255)
                                             gTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAAtct                                                                                        >  1:1184552/1‑150 (MQ=255)
                                              tataCCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGTCTTCGAAGAACTTGTGCTCGCACGCCGGGCCAAGGTACTCGCCCGCCCGGCGGTCGACCCCTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAAtctt                                                                                       <  2:1023392/150‑1 (MQ=255)
                                                   cTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAgg                                                                                  >  1:311765/1‑150 (MQ=255)
                                                   cTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAgg                                                                                  >  1:1793446/1‑150 (MQ=255)
                                                    tGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGt                                                                                 <  2:1489346/150‑1 (MQ=255)
                                                     ggTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTg                                                                                >  1:1236286/1‑150 (MQ=255)
                                                       tACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTgag                                                                              >  2:1691377/1‑150 (MQ=255)
                                                            cacGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCgg                                                                         <  2:1236286/150‑1 (MQ=255)
                                                               gCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCg                                                                                                                                      <  1:718587/86‑1 (MQ=255)
                                                               gCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCg                                                                                                                                      >  2:718587/1‑86 (MQ=255)
                                                                ccGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATc                                                                     >  2:793103/1‑150 (MQ=255)
                                                                ccGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATc                                                                     >  1:1745968/1‑150 (MQ=255)
                                                                      aTCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCAACTATGGCCTAATCTGCTAGGTGAGATCGGCATCTACcgc                                                               >  1:422748/1‑150 (MQ=255)
                                                                       tCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGc                                                                        <  2:534604/140‑1 (MQ=255)
                                                                       tCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGc                                                                        >  1:534604/1‑140 (MQ=255)
                                                                        ccACTTCGGCACCTTCGCCGACTTCGAAGAGCTGGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCAGTCGGCCGCGTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGt                                                             <  2:1314633/150‑1 (MQ=255)
                                                                          aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTca                                                           <  2:378993/150‑1 (MQ=255)
                                                                          aCTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGAACACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTca                                                           >  1:125014/1‑150 (MQ=255)
                                                                            ttCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCAcc                                                                                                                          <  2:752386/85‑1 (MQ=255)
                                                                            ttCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCAcc                                                                                                                          >  1:752386/1‑85 (MQ=255)
                                                                             tCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGTTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTcagca                                                        >  2:1094149/1‑150 (MQ=255)
                                                                             tCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTcagca                                                        <  1:1733577/150‑1 (MQ=255)
                                                                                   ccTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCgg                                                  <  1:179950/150‑1 (MQ=255)
                                                                                     ttCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACGGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGc                                                <  1:1859963/150‑1 (MQ=255)
                                                                                          cGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCa                                           <  2:125014/150‑1 (MQ=255)
                                                                                              ttCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCGCTTGCCCCGCAACGGGTGGGGTGGCCGGCTAGGGCCTAATCTTGTAGGTGAGATGGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAagag                                       >  1:678730/1‑150 (MQ=255)
                                                                                                 gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTc                                                                                                                                       <  2:1711375/51‑1 (MQ=255)
                                                                                                 gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTc                                                                                                                                       >  1:1711375/1‑51 (MQ=255)
                                                                                                 gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCTGGCCTGCAAGAGCAt                                    <  1:22441/150‑1 (MQ=255)
                                                                                                 gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCAt                                    >  1:307152/1‑150 (MQ=255)
                                                                                                  aagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAgg                                                                                  >  2:458293/1‑103 (MQ=255)
                                                                                                  aagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAgg                                                                                  <  1:458293/103‑1 (MQ=255)
                                                                                                  aagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATc                                   >  1:286754/1‑150 (MQ=255)
                                                                                                  aagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATc                                   >  2:1897050/1‑150 (MQ=255)
                                                                                                   agaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCa                                  >  2:576615/1‑150 (MQ=255)
                                                                                                     aaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGc                                <  2:1106141/150‑1 (MQ=255)
                                                                                                         ttGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTcg                             >  2:388495/1‑150 (MQ=255)
                                                                                                              gCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGt                         >  2:945705/1‑150 (MQ=255)
                                                                                                                aCGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCt                       <  1:1034016/150‑1 (MQ=255)
                                                                                                                  gAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGa                     >  1:198436/1‑150 (MQ=255)
                                                                                                                      gCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATg                                                                                                 >  1:1269824/1‑70 (MQ=255)
                                                                                                                      gCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATg                                                                                                 <  2:1269824/70‑1 (MQ=255)
                                                                                                                             cAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCgga          <  2:1184552/150‑1 (MQ=255)
                                                                                                                               aGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCgg                          <  2:1715297/132‑1 (MQ=255)
                                                                                                                               aGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCgg                          >  1:1715297/1‑132 (MQ=255)
                                                                                                                                ggTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACg                   <  1:58887/138‑1 (MQ=255)
                                                                                                                                ggTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACg                   >  2:58887/1‑138 (MQ=255)
                                                                                                                                ggTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCggagga       <  1:945705/150‑1 (MQ=255)
                                                                                                                                ggTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCggagga       >  2:1685339/1‑150 (MQ=255)
                                                                                                                                  tacgccaCCGCCTGGCGGTCGACCGCTTGCACCGCACCGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGATATCGGCATCTACCGCGTCAGCAGTCCGCGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAAt     <  1:507333/144‑1 (MQ=255)
                                                                                                                                     tCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGc  >  1:1821868/1‑150 (MQ=255)
                                                                                                                                     tCGACCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTcagca                                                        <  1:1785974/96‑1 (MQ=255)
                                                                                                                                     tCGACCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTcagca                                                        >  2:1785974/1‑96 (MQ=255)
                                                                                                                                     tCGACCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGc  >  1:1775403/1‑150 (MQ=255)
                                                                                                                                     tCGACCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGc  >  1:159284/1‑150 (MQ=255)
                                                                                                                                       gACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTc                                                     >  2:669558/1‑97 (MQ=255)
                                                                                                                                       gACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTc                                                     <  1:669558/97‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGC  >  NC_002947/5821813‑5822093

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