Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,007,720 C→A 5.5% S473I (AGC→ATC)  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,7200CA 5.5% 291.5 / 2.8 89S473I (AGC→ATC) PP_2628ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (44/40);  new base A (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.95e-01

ATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTCGCCCGGACGAATGTG  >  NC_002947/3007576‑3007868
                                                                                                                                                |                                                                                                                                                    
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                                                                                                                                   >  1:2250456/1‑144 (MQ=255)
                    aCTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGc                                                                                                                                   <  2:2250456/144‑1 (MQ=255)
                     cTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTc                                                                                                                            <  2:1963236/150‑1 (MQ=255)
                      tGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGAc                                                                                                                                         >  1:569416/1‑136 (MQ=255)
                      tGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGAc                                                                                                                                         <  2:569416/136‑1 (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                                                                                                                     <  2:528824/95‑1 (MQ=255)
                                                                                   gcgAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGt                                                                                                                     >  1:528824/1‑95 (MQ=255)
                                                                                    cgAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTgg                                                                                                                                    >  2:127497/1‑79 (MQ=255)
                                                                                    cgAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTgg                                                                                                                                    <  1:127497/79‑1 (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:720484/1‑150 (MQ=255)
                                                                                                               tCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCg                                  >  1:1969769/1‑150 (MQ=255)
                                                                                                                 cTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTcccc                                                                                                     <  1:741450/81‑1 (MQ=255)
                                                                                                                 cTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTcccc                                                                                                     >  2:741450/1‑81 (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                                      <  1:2377232/131‑1 (MQ=255)
                                                                                                                              ccGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGTGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCg                                      >  2:2377232/1‑131 (MQ=255)
                                                                                                                                gATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTTGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGATCGGCCGACCAGGCCGACTTTTTc                 <  2:89178/150‑1 (MQ=255)
                                                                                                                                gATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTc                 <  1:2109065/150‑1 (MQ=255)
                                                                                                                                    gcgcTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTCGccc             <  2:2262959/150‑1 (MQ=255)
                                                                                                                                      gcTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAgg                                                            >  1:1785581/1‑101 (MQ=255)
                                                                                                                                      gcTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAgg                                                            <  2:1785581/101‑1 (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)
                                                                                                                                                |                                                                                                                                                    
ATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCCTTGTCCACGTCGTAGAAGCGCAGCAGCAGGTTGATCAAGGTGGACTTGCCGGCGCCGGAGCGGCCGACCAGGCCGACTTTTTCGCCCGGACGAATGTG  >  NC_002947/3007576‑3007868

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