Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,007,723 T→G 5.6% D472A (GAT→GCT)  PP_2628 ← ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,007,7230TG 5.6% 271.2 / 5.3 89D472A (GAT→GCT) PP_2628ABC transporter ATP‑binding 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 = 8.51e-01

ATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTCGCCCGGACGAATGTGCAGC  >  NC_002947/3007576‑3007872
                                                                                                                                                   |                                                                                                                                                     
atgaAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATcc                                                                                                                                                     >  1:2220068/1‑150 (MQ=255)
    aaCCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCg                                                                                                                                                 >  1:2003882/1‑150 (MQ=255)
     aCCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTCGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGt                                                                                                                                                <  1:758582/150‑1 (MQ=255)
       cATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGa                                                                                                                                              >  2:1221038/1‑150 (MQ=255)
       cATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTACGTGa                                                                                                                                              >  2:1794551/1‑150 (MQ=255)
        aTCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGGAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATACCGATGGAGCTGCGCACTCTATGCttgctgac                                                                                                                                             >  2:1123850/1‑143 (MQ=255)
        aTCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGATAGCCTGCgcagc                                                                                                                                             <  1:850155/150‑5 (MQ=255)
            gCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCttc                                                                                                                                         >  2:1210254/1‑149 (MQ=255)
            gCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTg                                                                                                                                         <  2:1247561/150‑1 (MQ=255)
             ccTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGATCTTgg                                                                                                                                        <  2:94055/150‑1 (MQ=255)
              cTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGc                                                                                                                                       >  2:727485/1‑150 (MQ=255)
               tGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCg                                                                                                                                      >  2:2235286/1‑150 (MQ=255)
               tGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCg                                                                                                                                      >  2:865289/1‑150 (MQ=255)
                    aCTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGc                                                                                                                                       <  2:2250456/144‑1 (MQ=255)
                    aCTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGc                                                                                                                                       >  1:2250456/1‑144 (MQ=255)
                     cTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTc                                                                                                                                <  2:1963236/150‑1 (MQ=255)
                      tGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGAc                                                                                                                                             <  2:569416/136‑1 (MQ=255)
                      tGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGAc                                                                                                                                             >  1:569416/1‑136 (MQ=255)
                         ccGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGcc                                                                                                                            <  1:1296193/150‑1 (MQ=255)
                          cGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCg                                                                                                                                                 >  1:907531/1‑128 (MQ=255)
                          cGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCg                                                                                                                                                 <  2:907531/128‑1 (MQ=255)
                           gCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCg                                                                                                                          >  2:294597/1‑150 (MQ=255)
                                   aCGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCg                                                                                                                  >  2:1388234/1‑150 (MQ=255)
                                    cGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTgg                                                                                                                                        >  2:243450/1‑127 (MQ=255)
                                    cGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTgg                                                                                                                                        <  1:243450/127‑1 (MQ=255)
                                     ggATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTa                                                                                                                <  1:1388234/150‑1 (MQ=255)
                                       aTCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTAtt                                                                                                              >  1:808377/1‑150 (MQ=255)
                                         cTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGACGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTc                                                                                                            <  2:1094597/150‑1 (MQ=255)
                                          tgtgCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGATGCGCAGGCTATTCTGCCTGACCGGGGTGGCGATCTGTGCGGGGATGCGTATTTcc                                                                                                           <  2:25278/150‑1 (MQ=37)
                                           gtgCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTccc                                                                                                          >  1:1425438/1‑150 (MQ=255)
                                            tgCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTcccc                                                                                                         <  2:1165809/150‑1 (MQ=255)
                                            tgCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTcccc                                                                                                         >  1:992916/1‑150 (MQ=255)
                                                tCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTg                                                                                                     <  2:1553788/150‑1 (MQ=255)
                                                  ggTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTc                                                                                                   <  2:1786744/150‑1 (MQ=255)
                                                       cATCGGGGCGGCCATATGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCAcgt                                                                                              <  1:1243827/150‑1 (MQ=255)
                                                              gcggcCATAGGCAATATTGTCGCGTAAGCATCGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGaa                                                                                       <  1:611768/150‑1 (MQ=255)
                                                              gcggcCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGaa                                                                                       <  1:517766/150‑1 (MQ=255)
                                                                   cATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGt                                                                                                                         <  2:1012874/111‑1 (MQ=255)
                                                                   cATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGt                                                                                                                         >  1:1012874/1‑111 (MQ=255)
                                                                   cATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCgca                                                                                  >  2:1968178/1‑150 (MQ=255)
                                                                   cATAGGCAATATTGTCGCGAATAGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGt                                                                                                                         <  2:202965/111‑1 (MQ=255)
                                                                   cATAGGCAATATTGTCGCGAATAGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGt                                                                                                                         >  1:202965/1‑111 (MQ=255)
                                                                     tAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCgcagc                                                                                >  1:152961/1‑150 (MQ=255)
                                                                         cAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGAACATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCgcagcagca                                                                            >  1:970987/1‑150 (MQ=255)
                                                                          aaTATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGcagcagcag                                                                           >  1:2262959/1‑150 (MQ=255)
                                                                             atTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGtt                                                                        <  2:152961/150‑1 (MQ=255)
                                                                                 tCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGAGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCACCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATc                                                                    >  1:20080/1‑150 (MQ=255)
                                                                                   gcgAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGt                                                                                                                         >  1:528824/1‑95 (MQ=255)
                                                                                   gcgAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGt                                                                                                                         <  2:528824/95‑1 (MQ=255)
                                                                                    cgAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTgg                                                                                                                                        <  1:127497/79‑1 (MQ=255)
                                                                                    cgAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTgg                                                                                                                                        >  2:127497/1‑79 (MQ=255)
                                                                                    cgAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAg                                                                 >  2:697200/1‑150 (MQ=255)
                                                                                     gAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAgg                                                                <  1:727485/150‑1 (MQ=255)
                                                                                      aaTCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGt                                                               >  2:1570468/1‑150 (MQ=255)
                                                                                          gAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCTTCGTTCTGCCTGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGAc                                                           <  1:294597/150‑1 (MQ=255)
                                                                                           aaCGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACt                                                          <  1:697200/150‑1 (MQ=255)
                                                                                                 gCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCgg                                                    <  2:2003882/150‑1 (MQ=255)
                                                                                                   aaCAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGgcg                                                  >  1:89178/1‑150 (MQ=255)
                                                                                                     cAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGaa                                                                                       <  2:1405583/111‑1 (MQ=255)
                                                                                                     cAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGaa                                                                                       >  1:1405583/1‑111 (MQ=255)
                                                                                                      aGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCg                                               >  2:156861/1‑150 (MQ=255)
                                                                                                      aGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCg                                               >  2:668207/1‑150 (MQ=255)
                                                                                                             tgtCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGc                                        <  2:808377/150‑1 (MQ=255)
                                                                                                               tCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCg                                      >  1:1969769/1‑150 (MQ=255)
                                                                                                               tCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCg                                      >  1:720484/1‑150 (MQ=255)
                                                                                                                 cTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTcccc                                                                                                         >  2:741450/1‑81 (MQ=255)
                                                                                                                 cTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTcccc                                                                                                         <  1:741450/81‑1 (MQ=255)
                                                                                                                   gAGTGACCATCCCCATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGTCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCa                                  <  1:1123850/150‑1 (MQ=255)
                                                                                                                    aGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGTCCGACCAg                                 >  2:1869385/1‑150 (MQ=255)
                                                                                                                     gTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAgg                                <  2:453893/150‑1 (MQ=255)
                                                                                                                      tGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCaa                                                                  >  2:161988/1‑115 (MQ=255)
                                                                                                                      tGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCaa                                                                  <  1:161988/115‑1 (MQ=255)
                                                                                                                       gACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGcc                              <  1:1030302/150‑1 (MQ=255)
                                                                                                                         ccATGCCGATGGCGCTGTGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGa                            >  2:912635/1‑150 (MQ=255)
                                                                                                                          cATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGAc                           <  1:650836/150‑1 (MQ=255)
                                                                                                                          cATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAAGCCGAc                           >  2:242995/1‑150 (MQ=255)
                                                                                                                           aTGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACt                          >  1:393376/1‑150 (MQ=255)
                                                                                                                              ccGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGTGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCg                                          >  2:2377232/1‑131 (MQ=255)
                                                                                                                              ccGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGTGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCg                                          <  1:2377232/131‑1 (MQ=255)
                                                                                                                                gATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTTGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGATCGGCCGACCAGGCCGACTTTTTc                     <  2:89178/150‑1 (MQ=255)
                                                                                                                                gATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTc                     <  1:2109065/150‑1 (MQ=255)
                                                                                                                                    gcgcTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTCGccc                 <  2:2262959/150‑1 (MQ=255)
                                                                                                                                      gcTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAgg                                                                <  2:1785581/101‑1 (MQ=255)
                                                                                                                                      gcTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAgg                                                                >  1:1785581/1‑101 (MQ=255)
                                                                                                                                        cacgcAGGATAGCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCa                                  >  2:1054219/3‑129 (MQ=255)
                                                                                                                                        cacgcAGGATAGCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCa                                  <  1:1054219/127‑1 (MQ=255)
                                                                                                                                         acgcAGGATAGCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCCGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCa                                  >  1:1571309/2‑128 (MQ=255)
                                                                                                                                         acgcAGGATAGCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCa                                  <  2:1571309/127‑1 (MQ=255)
                                                                                                                                               gCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTCGCCCGGACGAAtgtg      >  2:1857278/1‑150 (MQ=255)
                                                                                                                                                   tCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTCGCCCGGACGAAt         <  2:1344620/143‑1 (MQ=255)
                                                                                                                                                   tCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTCGCCCGGACGAAt         >  1:1344620/1‑143 (MQ=255)
                                                                                                                                                   tCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTCGCCCGGACGAATGTGCAGc  <  1:668207/150‑1 (MQ=255)
                                                                                                                                                   tCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTCGCCCGGACGAATGTGCAGc  <  2:970987/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTCGCCCGGACGAATGTGCAGC  >  NC_002947/3007576‑3007872

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