Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,101,776 G→A 5.6% G291D (GGC→GAC)  uvrD → DNA‑dependent ATPase I/helicase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,101,7760GA 5.6% 213.2 / 2.1 71G291D (GGC→GAC) uvrDDNA‑dependent ATPase I/helicase II
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base A (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CGGTGGCGACAGCCTGATGGCGGTAGGCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGTCGTGGAAACCATC  >  NC_002947/6101627‑6101924
                                                                                                                                                      |                                                                                                                                                    
cgGTGGCGACAGCCTGATGGCGGTAGGCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACcggcgg                                                                                                                                                      >  2:1615205/1‑150 (MQ=255)
 ggTGGCGACAGCCTGATGGCGGTAGGCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATAAGCCGCCTGGAGCAGAACTACCGCTCCACCggcggc                                                                                                                                                     >  2:1124275/1‑150 (MQ=255)
     gCGACAGCCTGATGGCGGTAGGCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATcc                                                                                                                                                 <  1:538716/150‑1 (MQ=255)
       gACAGCCTGATGGCGGTAGGCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTc                                                                                                                                               >  2:1433765/1‑150 (MQ=255)
         cAGCCTGATGGCGGTAGGCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCaa                                                                                                                                             <  2:404814/150‑1 (MQ=255)
          aGCCTGATGGCGGTAGGCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAg                                                                                                                                            <  2:229724/150‑1 (MQ=255)
          aGCCTGATGGCGGTAGGCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGCAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAg                                                                                                                                            <  1:1124275/150‑1 (MQ=255)
           gCCTGATGGCGGTAGGCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAgg                                                                                                                                           >  2:86727/1‑150 (MQ=255)
                aTGGCGGTAGGCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGACAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGcc                                                                                                                                      >  1:931867/1‑150 (MQ=255)
                         ggCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTGACAAGTACACCGCCGACTTCCCCGCCGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTg                                                                                                                             >  1:347259/1‑150 (MQ=255)
                            gATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATc                                                                                                                          >  2:414540/1‑150 (MQ=255)
                              tGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGc                                                                                                                        >  2:357898/1‑150 (MQ=255)
                               gacgacCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTCGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGcc                                                                                                                       <  2:1718600/150‑1 (MQ=255)
                                  gacCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGGGCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCcaac                                                                                                                    >  1:1182021/1‑150 (MQ=255)
                                  gacCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCcaac                                                                                                                    >  2:982421/1‑150 (MQ=255)
                                    ccAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCTGCATCCTCAAGGCGGCCAACGCGCTGATCGCcaacaa                                                                                                                  >  2:1130662/1‑150 (MQ=255)
                                     cAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCaacaac                                                                                                                 >  2:623299/1‑150 (MQ=255)
                                       gTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAg                                                                                                               >  2:1316041/1‑150 (MQ=255)
                                       gTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCACCAACAACAg                                                                                                               <  2:1559821/150‑1 (MQ=255)
                                           aTCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGACATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCggg                                                                                                           >  2:331687/1‑150 (MQ=255)
                                           aTCTACGGCTGGCGCGGCGCCAAGATCAGGAACGTTCACGAGTACACCGCCGACTTCCCCGACAGCGTTTGACCCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCggg                                                                                                           <  1:1790375/150‑1 (MQ=255)
                                              tACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGgcgc                                                                                                        >  1:1448400/1‑150 (MQ=255)
                                                cGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGACATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCt                                                                                                      >  2:323882/1‑150 (MQ=255)
                                                   cTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTgtg                                                                                                   <  2:1448400/150‑3 (MQ=255)
                                                      ccgcggcgcCAAAATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCaa                                                                                                <  1:930840/149‑1 (MQ=255)
                                                        gcggcgcCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAgg                                                                                              >  2:821722/1‑150 (MQ=255)
                                                          ggcgcCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGaa                                                                                            <  1:623299/150‑1 (MQ=255)
                                                          ggcgcCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGaa                                                                                            <  1:1297481/150‑1 (MQ=255)
                                                             gcCAAGATCGAGAGCATTCACCAGTACACCGCCGACTTCCCCGAGCTGGGAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACtg                                                                                         >  2:274140/1‑150 (MQ=255)
                                                              ccAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACtgt                                                                                        <  1:1433765/150‑1 (MQ=255)
                                                               cAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACtgtg                                                                                       <  1:414540/150‑1 (MQ=255)
                                                                aaGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTgg                                                                                      >  1:1259725/1‑150 (MQ=255)
                                                                 aGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGa                                                                                     <  2:1780324/150‑1 (MQ=255)
                                                                 aGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGACATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGa                                                                                     <  1:323882/150‑1 (MQ=255)
                                                                  gATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCTCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGCTCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGAc                                                                                    <  1:1854977/150‑1 (MQ=255)
                                                                   aTCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGAcc                                                                                   >  2:1242735/1‑150 (MQ=255)
                                                                     cGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTATCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGa                                                                                 >  2:743793/1‑150 (MQ=255)
                                                                     cGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGa                                                                                 >  1:269234/1‑150 (MQ=255)
                                                                       agaACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCcaa                                                                                                                     >  1:1558286/1‑112 (MQ=255)
                                                                       agaACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCcaa                                                                                                                     <  2:1558286/112‑1 (MQ=255)
                                                                          aCATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGACATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATgg                                                                            <  1:331687/150‑1 (MQ=255)
                                                                          aCAGTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATgg                                                                            <  2:693303/150‑1 (MQ=255)
                                                                           cATTCACCAGTACACCGCCGACCTCCCCGACGCTGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATggg                                                                           <  2:1904108/150‑1 (MQ=255)
                                                                              tCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTc                                                                                                                                               <  1:870746/79‑1 (MQ=255)
                                                                              tCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTc                                                                                                                                               >  2:870746/1‑79 (MQ=255)
                                                                                 ccAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGggcgaagg                                                                     >  1:1482387/1‑150 (MQ=255)
                                                                                  cAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGggcgaaggc                                                                    >  1:1939669/1‑150 (MQ=255)
                                                                                   aGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGggcgaaggcc                                                                   >  2:1559518/1‑149 (MQ=255)
                                                                                     tACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGggcgaaggcgaa                                                                 >  1:229092/1‑150 (MQ=255)
                                                                                      acacCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAAc                                                                <  1:860999/150‑1 (MQ=255)
                                                                                      acacCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAAc                                                                <  1:1316041/150‑1 (MQ=255)
                                                                                        acCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCa                                                              >  2:1034120/1‑150 (MQ=255)
                                                                                        acCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCa                                                              <  2:1939669/150‑1 (MQ=255)
                                                                                          cgccgACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTg                                                                                                     <  2:1440516/109‑1 (MQ=255)
                                                                                          cgccgACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTg                                                                                                     >  1:1440516/1‑109 (MQ=255)
                                                                                           gccgACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTg                                                           >  2:1768586/1‑150 (MQ=255)
                                                                                                         cGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGc                                             <  2:1482387/150‑1 (MQ=255)
                                                                                                            cGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTa                                          >  2:1068401/1‑150 (MQ=255)
                                                                                                               aaTGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACaa                                       >  2:974186/1‑150 (MQ=255)
                                                                                                                  gATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGATCTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGa                                    <  2:1259725/150‑1 (MQ=255)
                                                                                                                         ccTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCacga                              >  1:955582/1‑150 (MQ=255)
                                                                                                                                  gAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAAACACTAAAACGATACGCGGCCTACAACGAGCACGACGAAgcgcg                     >  2:645918/1‑150 (MQ=255)
                                                                                                                                       aCCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATg                <  1:1242735/150‑1 (MQ=255)
                                                                                                                                        ccGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGt               >  2:771777/1‑150 (MQ=255)
                                                                                                                                         cGCTCCACCGGCTGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGTc              <  1:1130662/150‑1 (MQ=255)
                                                                                                                                          gCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGggcgaagg                                                                     <  1:1142378/94‑1 (MQ=255)
                                                                                                                                          gCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGggcgaagg                                                                     >  2:1142378/1‑94 (MQ=255)
                                                                                                                                           cTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGTCGt            <  1:974186/150‑1 (MQ=255)
                                                                                                                                           cTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGTCGt            <  1:771777/150‑1 (MQ=255)
                                                                                                                                            tCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGCGCGGTGGATCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGTCGTg           <  1:821722/150‑1 (MQ=255)
                                                                                                                                              cACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGTCGTGGa         >  2:78919/1‑150 (MQ=255)
                                                                                                                                               aCCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGTCGTGGaa        <  1:743793/150‑1 (MQ=255)
                                                                                                                                               aCCGGCGGCATCCTCAAAGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGTCGTGGaa        <  2:931867/150‑1 (MQ=255)
                                                                                                                                                   gcggcATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGggcgaaggcgaa                                                                 >  1:857381/1‑89 (MQ=255)
                                                                                                                                                   gcggcATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGggcgaaggcgaa                                                                 <  2:857381/89‑1 (MQ=255)
                                                                                                                                                   gcggcATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGTCGTGGAAACCa    <  1:1615205/150‑1 (MQ=255)
                                                                                                                                                     ggcATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGTCGTGGAAACCATc  <  2:229092/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGGTGGCGACAGCCTGATGGCGGTAGGCGATGACGACCAGTCGATCTACGGCTGGCGCGGCGCCAAGATCGAGAACATTCACCAGTACACCGCCGACTTCCCCGACGCCGAAATGATCCGCCTGGAGCAGAACTACCGCTCCACCGGCGGCATCCTCAAGGCGGCCAACGCGCTGATCGCCAACAACAGCGGGCGCCTGGGCAAGGAACTGTGGACCGACATGGGCGAAGGCGAACCACTGACCCTGTACGCGGCCTACAACGAGCACGACGAAGCGCGCTATGTCGTGGAAACCATC  >  NC_002947/6101627‑6101924

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: