Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,111,244 G→C 5.1% P108A (CCA→GCA)  PP_5361 ← protein P47K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,111,2440GC 5.1% 250.7 / 4.3 79P108A (CCA→GCA) PP_5361protein P47K
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base C (2/2);  total (39/40)
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.91e-01

CGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAGCGCGG  >  NC_002947/6111095‑6111392
                                                                                                                                                     |                                                                                                                                                    
cGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTCACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGAGCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTgg                                                                                                                                                      <  1:1590372/150‑1 (MQ=255)
  tttCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCa                                                                                                                                                    <  1:499255/150‑1 (MQ=255)
      gCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTgg                                                                                                                                                <  2:905205/150‑1 (MQ=255)
      gCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTgg                                                                                                                                                <  2:905119/150‑1 (MQ=255)
      gCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTgg                                                                                                                                                <  2:187098/150‑1 (MQ=255)
        cccGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCt                                                                                                                                              >  2:1604813/1‑150 (MQ=255)
          cGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCa                                                                                                                                            >  2:1210209/1‑150 (MQ=255)
             cttgCCAGTCTGTCTGCTGCGTGAAAGTCGCACAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCagag                                                                                                                                         <  2:490360/147‑1 (MQ=39)
              tGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAgaga                                                                                                                                        >  1:1756434/1‑150 (MQ=255)
                  cAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAGGTTCAAGCCAGCGACCACGGTGACCATGGTGTCCAGTCGCGCCAGATCGGACAGACTGCGAGCATCCTCGTCGCGGAAAGTGAAGGTCTCGGCCACTGGCAGGGGCTCAGAGATAcc                                                                                                                                    >  2:656447/1‑150 (MQ=38)
                     gCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGt                                                                                                                                 <  1:1720388/150‑1 (MQ=255)
                     acaccctgctgCCTGATATTCGCTCAGTAAGTTCAAGACATCGCCCACGGTGACCATGGTGTCCAGTCGCACCATAACGGACTGTTTGCCACCCTGGTCGTGGCCGAAAGTGAAGGTCTCATCCACTGTCAGTGTCTCAGAGATACCTGt                                                                                                                                 <  2:585298/145‑1 (MQ=255)
                      cTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGCCAGTGGCTgaga                                                                                                                                          >  1:850724/1‑140 (MQ=255)
                      cTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGCCAGTGGCTgaga                                                                                                                                          <  2:850724/140‑1 (MQ=255)
                        gTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGa                                                                                                                              <  1:1241397/150‑1 (MQ=255)
                           tgctgcCTGATAGTCGCGCAGGAAGTTCGAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTc                                                                                                                           <  2:118888/150‑1 (MQ=255)
                             ctgcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGTGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGa                                                                                                                         <  1:951685/150‑1 (MQ=255)
                              tgcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGCTGACCATGGTGTCAATAATCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGAt                                                                                                                        <  1:1261586/150‑1 (MQ=255)
                                ccTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTgg                                                                                                                                                      <  2:354694/118‑1 (MQ=255)
                                ccTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTgg                                                                                                                                                      >  1:354694/1‑118 (MQ=255)
                                ccTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCa                                                                                                                      <  1:418870/150‑1 (MQ=255)
                                     tAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCaaa                                                                                                                 <  1:1604813/150‑1 (MQ=255)
                                     tAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCaaa                                                                                                                 <  1:148335/150‑1 (MQ=255)
                                         cgcgCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAaataat                                                                                                             <  1:1210209/150‑1 (MQ=255)
                                         cgcgCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAaataat                                                                                                             <  1:320256/150‑1 (MQ=255)
                                          gcgcAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATc                                                                                                            >  2:1158500/1‑150 (MQ=255)
                                          gcgcAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATc                                                                                                            >  1:1512866/1‑150 (MQ=255)
                                                 aaGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGc                                                                                                     <  2:886186/150‑1 (MQ=255)
                                                       aaGCCATCGACCACGGTGACCATTGTGTCCAGTCGCGCCATATCGGACAGACTGCAACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATAAGCAAATAATCGAAACGCCCctct                                                                                               >  1:647562/1‑150 (MQ=255)
                                                        aGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCtt                                                                                              >  2:463608/1‑150 (MQ=255)
                                                         gCCATCGACCACGGTGCCCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTc                                                                                             <  1:187720/150‑1 (MQ=255)
                                                          ccATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCg                                                                                            <  1:908982/150‑1 (MQ=255)
                                                           cATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCgg                                                                                           <  1:604585/150‑1 (MQ=255)
                                                           cATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGCCAGTGGCTgagacc                                                                                                                                        <  1:795187/105‑3 (MQ=255)
                                                           cATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGCCAGTGGCTgagacc                                                                                                                                        >  2:795187/1‑103 (MQ=255)
                                                            aTCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGc                                                                                          >  2:1540243/1‑150 (MQ=255)
                                                             tCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGcc                                                                                         >  2:1168280/1‑150 (MQ=255)
                                                                 ccACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTc                                                                                     >  2:1330900/1‑150 (MQ=255)
                                                                   aCGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTcgc                                                                                   <  2:1167144/150‑1 (MQ=255)
                                                                     ggTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATATCGGTGGACTCGATCAGCAaataat                                                                                                             >  2:517312/1‑122 (MQ=255)
                                                                     ggTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATATCGGTGGACTCGATCAGCAaataat                                                                                                             <  1:517312/122‑1 (MQ=255)
                                                                             tGGTGTCCAGTCGCGCCATATCGGACAGGCTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAAcctcct                                                                         <  2:576760/150‑1 (MQ=255)
                                                                              ggTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACctcctc                                                                        >  2:1362822/1‑150 (MQ=255)
                                                                                 gtCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGac                                                                     >  2:104116/1‑149 (MQ=255)
                                                                                  tCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGc                                                                    >  1:754090/1‑150 (MQ=255)
                                                                                      gTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATc                                                                                                            <  1:1420529/106‑1 (MQ=255)
                                                                                      gTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATc                                                                                                            >  2:1420529/1‑106 (MQ=255)
                                                                                      gTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGt                                                                >  2:345633/1‑150 (MQ=255)
                                                                                               tatCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCa                                                       >  2:313122/1‑150 (MQ=255)
                                                                                               tatCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCa                                                       >  2:1672700/1‑150 (MQ=255)
                                                                                               tatCGGACAGACTGCAAACCTACTCATAACGGACAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCa                                                       >  2:313097/1‑150 (MQ=255)
                                                                                                atCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCaa                                                      <  1:1540243/150‑1 (MQ=255)
                                                                                                   ggACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGt                                                   >  2:178324/1‑150 (MQ=255)
                                                                                                    gACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCTGCCAGTCGCGCAACCTCCTCGGGCATGTCTTCGCGCAACGTg                                                  >  2:969081/1‑150 (MQ=255)
                                                                                                    gACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTg                                                  >  1:1748642/1‑150 (MQ=255)
                                                                                                           tGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCaaa                                           >  2:1250831/1‑150 (MQ=255)
                                                                                                               aCCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCaaa                                           >  2:425854/1‑146 (MQ=255)
                                                                                                               aCCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCaaa                                           <  1:425854/146‑1 (MQ=255)
                                                                                                                cccTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCtt                                                             <  2:192392/127‑1 (MQ=255)
                                                                                                                cccTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCtt                                                             >  1:192392/1‑127 (MQ=255)
                                                                                                                cccTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATNCAg                                      <  2:1194408/150‑1 (MQ=255)
                                                                                                                  cTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTGAGAGCTACCGGTGGACTCGATCAGCAAATAATCGAAACGCACCTCTTCGGCCAGGCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGcc                                    >  1:1470532/1‑150 (MQ=255)
                                                                                                                   tGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGcc                                    <  1:1492864/149‑1 (MQ=255)
                                                                                                                   tGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGcc                                    >  2:1492864/1‑149 (MQ=255)
                                                                                                                     cTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGtt                                 <  1:451073/150‑1 (MQ=255)
                                                                                                                      tcgtcgCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTg                                >  1:1575527/1‑150 (MQ=255)
                                                                                                                           gcgGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCCGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCAt                           >  1:725519/1‑150 (MQ=255)
                                                                                                                              gAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCAtctc                        <  1:463608/150‑1 (MQ=255)
                                                                                                                                 aGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCTCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGAc                     >  1:794413/1‑150 (MQ=255)
                                                                                                                                    gAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTc                                                                                             >  2:1669226/1‑75 (MQ=255)
                                                                                                                                    gAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTc                                                                                             <  1:1669226/75‑1 (MQ=255)
                                                                                                                                    gAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAg                  >  1:1791626/1‑150 (MQ=255)
                                                                                                                                    gAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAg                  >  2:1029159/1‑150 (MQ=255)
                                                                                                                                         tctcAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTtct             >  2:646807/1‑150 (MQ=255)
                                                                                                                                         tctcAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTtct             >  2:1697620/1‑150 (MQ=255)
                                                                                                                                          ctcAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTtctc            <  2:1791626/150‑1 (MQ=255)
                                                                                                                                             aGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTc         <  1:599434/150‑1 (MQ=255)
                                                                                                                                             aGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTc         <  1:313097/150‑1 (MQ=255)
                                                                                                                                             aGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTc         <  1:313122/150‑1 (MQ=255)
                                                                                                                                              gCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCTCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCa        <  2:794413/150‑1 (MQ=255)
                                                                                                                                                   tGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCa                                                                                        >  2:1304799/1‑65 (MQ=255)
                                                                                                                                                   tGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCa                                                                                        <  1:1304799/65‑1 (MQ=255)
                                                                                                                                                    ggCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAGCGCgg  <  2:1748642/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAGCGCGG  >  NC_002947/6111095‑6111392

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: ‑