Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,183,216 T→G 5.6% H189P (CAC→CCC)  PP_4563 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,183,2160TG 5.6% 265.0 / 6.2 89H189P (CAC→CCC) PP_4563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (44/40);  new base G (2/3);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

TCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCTTGTAGTCGGCCGGGGTGACA  >  NC_002947/5183067‑5183365
                                                                                                                                                     |                                                                                                                                                      
tctTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGt                                                                                                                                                        >  1:592510/1‑150 (MQ=255)
tctTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGt                                                                                                                                                        <  2:901842/150‑1 (MQ=255)
tctTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGt                                                                                                                                                        <  2:1248168/150‑1 (MQ=255)
tctTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGt                                                                                                                                                        <  2:1527496/150‑1 (MQ=255)
tctTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAAACTTTCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGt                                                                                                                                                        >  2:571744/1‑150 (MQ=255)
    gTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGCAGCGCGTCTCTTGCAt                                                                                                                                                    <  2:599820/150‑1 (MQ=255)
     tCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGCGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATg                                                                                                                                                   >  1:1329261/1‑150 (MQ=255)
               ccAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGGGCATGCACGAgaagg                                                                                                                                         <  2:420848/150‑5 (MQ=255)
               ccAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGGGCATGCACGAgaag                                                                                                                                          <  1:1342657/149‑4 (MQ=255)
               ccAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGGGCATGCACGAgaag                                                                                                                                          >  2:1342657/1‑146 (MQ=255)
                      cccGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCCGGGTGTCAGCCCGGGAGGCGGCTACCGGCTGCAGTCTCAACGAAATCGTTCCGTCCGGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACtt                                                                                                                                  >  2:1478728/1‑150 (MQ=255)
                      cccGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACtt                                                                                                                                  >  2:447653/1‑150 (MQ=255)
                      cccGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACtt                                                                                                                                  >  2:1409688/1‑150 (MQ=255)
                      cccGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACtt                                                                                                                                  >  2:1754399/1‑150 (MQ=255)
                       ccGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTg                                                                                                                                 >  1:835660/1‑150 (MQ=255)
                        cGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGa                                                                                                                                >  1:1148387/1‑150 (MQ=255)
                          ccTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGCGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAac                                                                                                                              <  2:1329261/150‑1 (MQ=255)
                             tCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACAcc                                                                                                                           <  1:1754399/150‑1 (MQ=255)
                               ggCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTc                                                                                                                         <  1:1973126/150‑1 (MQ=255)
                                gCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCg                                                                                                                        >  2:1864457/1‑150 (MQ=255)
                                  ttCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGtt                                                                                                                      <  1:1827955/150‑1 (MQ=255)
                                  ttCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGtt                                                                                                                      >  1:253317/1‑150 (MQ=255)
                                   cccAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCTGTCTCAGTGAAATCTTGCGGTCCGCGCGCACTCCCTTGCTGACCCCCTTCTCGTTCAAGCGCGAGCGCGAGAACTTGACCATCTCGtgc                                                                                                                     <  2:1483570/149‑3 (MQ=255)
                                     cAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCtt                                                                                                                   >  2:1862310/1‑150 (MQ=255)
                                      aGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTtg                                                                                                                  <  1:960303/150‑1 (MQ=255)
                                              cATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAg                                                                                                          >  1:666579/1‑150 (MQ=255)
                                              cATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAg                                                                                                          >  2:1737427/1‑150 (MQ=255)
                                               aTGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAgg                                                                                                         >  2:1243648/1‑150 (MQ=255)
                                                tGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGaa                                                                                                                                     >  2:933376/1‑121 (MQ=255)
                                                tGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGc                                                                                                        <  2:401366/150‑1 (MQ=255)
                                                tGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCACTTCTCGTGCATGCCCGAGCGCAAGaa                                                                                                                                     <  1:933376/121‑1 (MQ=255)
                                                  aTCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGaa                                                                                                                               >  1:690112/1‑125 (MQ=255)
                                                  aTCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGaa                                                                                                                               <  2:690112/125‑1 (MQ=255)
                                                  aTCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGccc                                                                                                      <  2:666579/150‑1 (MQ=255)
                                                           tGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCgg                                                                                             >  2:1915406/1‑150 (MQ=255)
                                                            gCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGa                                                                                                                                              >  1:1106400/1‑100 (MQ=255)
                                                            gCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGa                                                                                                                                              <  2:1106400/100‑1 (MQ=255)
                                                              aGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGgttg                                                                                          >  1:919094/1‑150 (MQ=255)
                                                                ggTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGgttgtt                                                                                        <  1:447653/150‑1 (MQ=255)
                                                                    tCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCtt                                                                                                                   <  1:1154895/119‑1 (MQ=255)
                                                                    tCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCtt                                                                                                                   >  2:1154895/1‑119 (MQ=255)
                                                                    tCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTtgatg                                                                                   <  2:253317/150‑1 (MQ=255)
                                                                     ccGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCGTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTtgatga                                                                                  <  1:1737427/150‑1 (MQ=255)
                                                                      cGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCATGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTGGTGCATTCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGGGGGTGTTgatgat                                                                                 >  1:2298416/1‑150 (MQ=255)
                                                                       gCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGatgatg                                                                                >  2:376917/1‑150 (MQ=255)
                                                                          aGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACGTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGcc                                                                             <  2:592510/150‑1 (MQ=255)
                                                                             gcggcgCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACTCCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTCATCCTGACCTTg                                                                          <  2:835660/150‑1 (MQ=255)
                                                                                    cTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCa                                                                   >  2:888423/1‑150 (MQ=255)
                                                                                     tACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTGTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCAt                                                                  >  1:1366337/1‑150 (MQ=255)
                                                                                              cAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTc                                                         <  1:2337276/150‑1 (MQ=255)
                                                                                               aGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTACTTGATAAAGCCCTTCTCGTGCATGCCCGAGCGTTGGAACTTGAAAACCCCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCg                                                        >  1:1375911/1‑150 (MQ=255)
                                                                                                ggCTCAGTGAAATCTTGCCGTCCTGGCGCACTTGCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGAGGATGGCCTTGAAGCCCATCGGCGTTTCGc                                                       <  2:2252549/150‑1 (MQ=255)
                                                                                                 gCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTtgt                                                                                                                 <  2:1518493/92‑1 (MQ=255)
                                                                                                 gCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTtgt                                                                                                                 >  1:1518493/1‑92 (MQ=255)
                                                                                                 gCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGcc                                                      >  1:1624753/1‑150 (MQ=255)
                                                                                                     aGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGt                                                                                                                                                        >  2:1876854/1‑49 (MQ=255)
                                                                                                     aGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGt                                                                                                                                                        <  1:1876854/49‑1 (MQ=255)
                                                                                                       tGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGccac                                                >  2:1312020/1‑150 (MQ=255)
                                                                                                            tCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTtga                                                                                     >  2:1655319/1‑108 (MQ=255)
                                                                                                            tCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTtga                                                                                     <  1:1655319/108‑1 (MQ=255)
                                                                                                            tCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACcagca                                           <  1:749148/150‑1 (MQ=255)
                                                                                                                   gTCCTGGCGCTCTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGcc                                                      >  2:1299316/1‑132 (MQ=255)
                                                                                                                   gTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGa                                                                         >  1:591561/1‑113 (MQ=255)
                                                                                                                   gTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGa                                                                         <  2:591561/113‑1 (MQ=255)
                                                                                                                   gTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGcc                                                      <  1:1299316/132‑1 (MQ=255)
                                                                                                                   gTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGAc                                    >  2:1748770/1‑150 (MQ=255)
                                                                                                                   gTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGAc                                    <  1:571744/150‑1 (MQ=255)
                                                                                                                       tGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGc                                >  2:1500933/1‑150 (MQ=255)
                                                                                                                        ggCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCt                               <  1:1862310/150‑1 (MQ=255)
                                                                                                                         gcgcACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTtgt                                                                                                                 <  1:1850023/68‑1 (MQ=255)
                                                                                                                         gcgcACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTtgt                                                                                                                 >  2:1850023/1‑68 (MQ=255)
                                                                                                                           gcACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGc                            >  1:671660/1‑150 (MQ=255)
                                                                                                                           gcACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGc                            <  1:1864457/150‑1 (MQ=255)
                                                                                                                               ttCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTtgtg                                                                                                                >  2:1786601/1‑63 (MQ=255)
                                                                                                                               ttCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTtgtg                                                                                                                <  1:1786601/63‑1 (MQ=255)
                                                                                                                                  cTTGATGAAGCCCTTCTCGTGCATGCCTGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGCTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGCCCACCAGCAGCTCGACCGGCTGGCCAACCtt                     <  2:2298416/150‑1 (MQ=255)
                                                                                                                                  cTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGt                                                                                                  >  1:446857/1‑74 (MQ=255)
                                                                                                                                  cTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGt                                                                                                  <  2:446857/74‑1 (MQ=255)
                                                                                                                                  cTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCtt                     <  2:919094/150‑1 (MQ=255)
                                                                                                                                  cTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCtt                     >  1:248348/1‑150 (MQ=255)
                                                                                                                                   ttGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCTTg                    <  1:1312020/150‑1 (MQ=255)
                                                                                                                                   ttGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCTTg                    >  2:1113804/1‑150 (MQ=255)
                                                                                                                                   ttGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACAAGCAGCTCGACCGGCTGGCCAACCTTg                    <  1:1915406/150‑1 (MQ=255)
                                                                                                                                       tgaAGCCCTGCTCGTGCAGGCCCGTGGGCAAGTGCTTGAACACCTCGTTCTTGTGGATCAGGCCGCCGTGGCGCTTGTTGATGATGGCCTTGAACCGCATCGGCGTCTCGCGGGCCACCAGCCTGTCGACCGGCTGGCCTACCTTGTAGt                >  2:681226/1‑150 (MQ=255)
                                                                                                                                         aaGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCTTGTAGTCg              <  1:376917/150‑1 (MQ=255)
                                                                                                                                          aGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGAc                                    >  1:1019359/1‑127 (MQ=255)
                                                                                                                                          aGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGAc                                    <  2:1019359/127‑1 (MQ=255)
                                                                                                                                          aGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTGCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCTTGTAGTCgg             >  2:1502045/1‑150 (MQ=255)
                                                                                                                                           gCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCTTGTAGTCGGc            >  2:380021/1‑150 (MQ=255)
                                                                                                                                               ttCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTTGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCTTGTAGTCGGCCggg        <  2:542740/150‑1 (MQ=255)
                                                                                                                                               ttCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCTTGTAGTCGGCCggg        <  1:1113804/150‑1 (MQ=255)
                                                                                                                                                gctcGGGCATGCACGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTtgat                                                                                    <  1:344052/72‑1 (MQ=255)
                                                                                                                                                gctcGGGCATGCACGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTtgat                                                                                    >  2:344052/2‑73 (MQ=255)
                                                                                                                                                  tcGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGCAAGAAGCCTTGAAGCCCATCGGCGTTTCGCCGGCCCCCAGCAGCTCGACCGGCTGGCCAACCTTGTAGTCGGCCgggggg     >  1:1740137/1‑148 (MQ=255)
                                                                                                                                                    gTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTc                                                         <  1:322466/96‑1 (MQ=255)
                                                                                                                                                    gTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTc                                                         >  2:322466/1‑96 (MQ=255)
                                                                                                                                                     tGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCTTGTAGTCGGCCGGGGTGACa  >  1:1287697/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTTGATGAAGCCCTTCTCGTGCATGCCCGAGCGCAAGAACTTGAACACCTCGTTCTTGTGGATCAGGCCCCAGTGGCGGTTGTTGATGATGGCCTTGAAGCCCATCGGCGTTTCGCCGGCCACCAGCAGCTCGACCGGCTGGCCAACCTTGTAGTCGGCCGGGGTGACA  >  NC_002947/5183067‑5183365

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