Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,008,805 A→C 6.2% I137S (ATC→AGC)  PP_3536 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,008,8050AC 6.2% 240.6 / 6.2 82I137S (ATC→AGC) PP_3536hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/46);  new base C (2/3);  total (33/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

TCGCGGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGATCTGCTTTGGCGTGTCTG  >  NC_002947/4008661‑4008954
                                                                                                                                                |                                                                                                                                                     
tcGCGGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGGTCTTg                                                                                                                                                  >  1:1390878/1‑150 (MQ=255)
 cgcgGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGa                                                                                                                                                 <  2:1545378/150‑1 (MQ=255)
  gcgGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGAt                                                                                                                                                <  1:1809821/150‑1 (MQ=255)
   cgGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGATCGACGCTCTTGagc                                                                                                                                               <  1:1421711/150‑3 (MQ=255)
   cgGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCATATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGAtt                                                                                                                                               <  2:2276665/150‑1 (MQ=255)
    gggggAAAGTATATATATTGTCATTCTCAATCACATTTTTTTCCTCGCTGAGGTCCTTGACTGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGAttt                                                                                                                                              <  2:265897/150‑1 (MQ=255)
       tgAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGaa                                                                                                                                           <  1:1307530/149‑1 (MQ=255)
       ggAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGaa                                                                                                                                           <  1:1288123/150‑1 (MQ=255)
        gAAAGGATATTTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTCGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTTTAGATCAAGAGTTTTGATCATGATTTGAAg                                                                                                                                          >  2:2117924/1‑150 (MQ=255)
            gTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGcc                                                                                                                                      <  1:1147534/150‑1 (MQ=255)
             tataTCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCt                                                                                                                                     >  2:1590349/1‑150 (MQ=255)
              atatCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTc                                                                                                                                    <  2:1009249/150‑1 (MQ=255)
                 tCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCAc                                                                                                                                 >  1:1876647/1‑150 (MQ=255)
                   tAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCt                                                                                                                               >  1:121299/1‑150 (MQ=255)
                      tGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAg                                                                                                                            <  2:1405741/150‑1 (MQ=255)
                      tGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAg                                                                                                                            >  2:1193871/1‑150 (MQ=255)
                        tATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAgc                                                                                                                                         <  1:1618572/135‑1 (MQ=255)
                        tATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAgc                                                                                                                                         >  2:1618572/1‑135 (MQ=255)
                             cTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGgaga                                                                                                                     <  1:1136000/150‑1 (MQ=255)
                                aGCATCATCTTTATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGActc                                                                                                                  <  1:463820/150‑1 (MQ=255)
                                  catcatCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTg                                                                                                                <  2:2043784/150‑1 (MQ=255)
                                  catcatCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTg                                                                                                                >  2:1812949/1‑150 (MQ=255)
                                   atcatcTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTgg                                                                                                               <  1:356333/150‑1 (MQ=255)
                                    tcatcTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGgt                                                                                                              <  1:2009135/150‑1 (MQ=255)
                                         ttCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGATCGACGCTCTTGa                                                                                                                                                 >  2:1039534/1‑110 (MQ=255)
                                         ttCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGATCGACGCTCTTGa                                                                                                                                                 <  1:1039534/110‑1 (MQ=255)
                                         ttCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGtctc                                                                                                         <  1:1635640/150‑1 (MQ=255)
                                                  gcTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCaa                                                                                                >  1:506448/1‑150 (MQ=255)
                                                  gcTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCaa                                                                                                <  2:1876647/150‑1 (MQ=255)
                                                       gggCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAg                                                                                                                                          <  1:471842/103‑1 (MQ=255)
                                                       gggCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAg                                                                                                                                          >  2:471842/1‑103 (MQ=255)
                                                       gggCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCAt                                                                                           <  2:58076/150‑1 (MQ=255)
                                                          ccTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAg                                                                                        <  1:1590349/150‑1 (MQ=255)
                                                            tGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGc                                                                                      <  2:1827871/150‑1 (MQ=255)
                                                                 ggCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACt                                                                                 >  2:1009406/1‑150 (MQ=255)
                                                                    ttGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGATCGACGCTCTTGagc                                                                                                                                               >  2:901871/1‑83 (MQ=255)
                                                                    ttGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGATCGACGCTCTTGagc                                                                                                                                               <  1:901871/85‑3 (MQ=255)
                                                                      gAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGga                                                                                                                       <  2:2171225/107‑1 (MQ=255)
                                                                      gAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGga                                                                                                                       >  1:2171225/1‑107 (MQ=255)
                                                                      gAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGcaa                                                                            <  1:2357553/150‑1 (MQ=255)
                                                                      gAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGaaa                                                                            >  1:925350/1‑150 (MQ=255)
                                                                      gAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGaaa                                                                            <  1:2155315/150‑1 (MQ=255)
                                                                      gAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTAGTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGaaa                                                                            <  1:538133/150‑1 (MQ=255)
                                                                          gACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGa                                                                        <  1:1259360/150‑1 (MQ=255)
                                                                          gACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGa                                                                        <  2:121299/150‑1 (MQ=255)
                                                                          gACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGa                                                                        <  1:1009406/150‑1 (MQ=255)
                                                                             cTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTTACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATAc                                                                     <  2:1646264/150‑1 (MQ=255)
                                                                                 tCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACtt                                                                                <  2:833100/135‑1 (MQ=255)
                                                                                 tCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACtt                                                                                >  1:833100/1‑135 (MQ=255)
                                                                                       cGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAAc                                                           >  2:630528/1‑150 (MQ=255)
                                                                                        gagaTTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGa                                                                                                                                                 <  1:1596461/63‑1 (MQ=255)
                                                                                        gagaTTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGa                                                                                                                                                 >  2:1596461/1‑63 (MQ=255)
                                                                                               tATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCa                                                   >  1:1199357/1‑150 (MQ=255)
                                                                                                 tCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGca                                                                                     <  2:2212695/114‑1 (MQ=255)
                                                                                                 tCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGca                                                                                     >  1:2212695/1‑114 (MQ=255)
                                                                                                 tCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACg                                                 >  2:929569/1‑150 (MQ=255)
                                                                                                   gCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACggg                                               >  2:995381/1‑150 (MQ=255)
                                                                                                    cTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGc                                              <  2:506448/150‑1 (MQ=255)
                                                                                                          cTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGa                                                                        >  1:1758256/1‑118 (MQ=255)
                                                                                                          cTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGa                                                                        <  2:1758256/118‑1 (MQ=255)
                                                                                                             aCGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGtc                                                                                                           >  2:1605862/1‑80 (MQ=255)
                                                                                                             aCGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGtc                                                                                                           <  1:1605862/80‑1 (MQ=255)
                                                                                                                  tATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTcaca                                >  2:2184721/1‑150 (MQ=255)
                                                                                                                    ttGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGa                                                                        <  2:533815/108‑1 (MQ=255)
                                                                                                                    ttGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGa                                                                        >  1:533815/1‑108 (MQ=255)
                                                                                                                      gACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAg                               >  1:804459/1‑147 (MQ=255)
                                                                                                                      gACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAg                               <  2:804459/147‑1 (MQ=255)
                                                                                                                        cGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCt                          <  1:1193871/150‑1 (MQ=255)
                                                                                                                         gggCCTGTCGCTCAAGAGCGTCGATCTTTATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTg                         <  1:1812949/150‑1 (MQ=255)
                                                                                                                            ccTGTCGCTCAAGATGATCGATCTTGATTTGACGCGCCTCCACCTCAGCTAGAGACTCTGGTGTCTCAGGCTCCAAATCCTGCGGCACACTTCAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACACGGTCTGTaa                      <  2:1199357/150‑1 (MQ=255)
                                                                                                                              tGTCGCTCAAGAGCGTCGATCTTGATTTGCAGCGCATCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGa                    >  1:2268872/1‑150 (MQ=255)
                                                                                                                                  gCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGaaa                                                                            <  1:2252226/90‑1 (MQ=255)
                                                                                                                                  gCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGaaa                                                                            >  2:2252226/1‑90 (MQ=255)
                                                                                                                                      aaGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGATCTGCttt            >  1:2260984/1‑150 (MQ=255)
                                                                                                                                      aaGAGCGGCGATCTTGATTTGAATCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGTTCCAAGTCATGAGGCACACTTGACTTAGCTCCAGTCCTCAACCCTCGCGACGGGCCACACTCACTCACAGGGCCTGTAAGATCTGCctt            <  2:695396/150‑1 (MQ=255)
                                                                                                                                       agagCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGATCTGCTTTg           <  1:630528/150‑1 (MQ=255)
                                                                                                                                       agagCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTAGAGACTCTCGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGATCTGCTTTg           <  1:2069766/150‑1 (MQ=255)
                                                                                                                                        gagCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGAt                                                                       <  2:180800/89‑1 (MQ=255)
                                                                                                                                        gagCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGAt                                                                       >  1:180800/1‑89 (MQ=255)
                                                                                                                                            gTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGATCTGCTTTGGCgtg      >  1:1605661/1‑150 (MQ=255)
                                                                                                                                              cGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGATCTGCTTTGGCGTGTc    <  2:1390878/150‑1 (MQ=255)
                                                                                                                                               gATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGATCTGCTTTGGCGTGTCt   >  1:849878/1‑150 (MQ=255)
                                                                                                                                                aTCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGt                            <  2:38952/124‑1 (MQ=255)
                                                                                                                                                aTCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGt                            >  1:1996035/1‑124 (MQ=255)
                                                                                                                                                aTCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGt                            >  1:38952/1‑124 (MQ=255)
                                                                                                                                                aTCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGt                            <  2:1996035/124‑1 (MQ=255)
                                                                                                                                                aTCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGATCTGCTTTGGCGTGTCTg  <  1:2184721/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TCGCGGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGATCTGCTTTGGCGTGTCTG  >  NC_002947/4008661‑4008954

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