Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,008,802 T→A 6.0% D138V (GAC→GTC)  PP_3536 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,008,8020TA 6.0% 259.7 / 5.8 84D138V (GAC→GTC) PP_3536hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/46);  new base A (2/3);  total (34/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CTTCATCATCGCGGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGATCTGCTTTGGCGTG  >  NC_002947/4008653‑4008950
                                                                                                                                                     |                                                                                                                                                    
cTTCATCATCGCGGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGt                                                                                                                                                      >  2:1136000/1‑150 (MQ=255)
 ttCATCATCGCGGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTc                                                                                                                                                     >  2:1259360/1‑150 (MQ=255)
 ttCATCATCGCGGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTc                                                                                                                                                     <  1:1460935/150‑1 (MQ=255)
        tcGCGGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTc                                                                                                                                                     <  1:974857/143‑1 (MQ=255)
        tcGCGGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTc                                                                                                                                                     >  2:974857/1‑143 (MQ=255)
        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                                                                                                                                     >  2:1618572/1‑135 (MQ=255)
                                tATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAgc                                                                                                                                     <  1:1618572/135‑1 (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                                                                                                                                             <  1:1039534/110‑1 (MQ=255)
                                                 ttCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGATCGACGCTCTTGa                                                                                                                                             >  2:1039534/1‑110 (MQ=255)
                                                 ttCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGtctc                                                                                                     <  1:1635640/150‑1 (MQ=255)
                                                          gcTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCaa                                                                                            <  2:1876647/150‑1 (MQ=255)
                                                          gcTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCaa                                                                                            >  1:506448/1‑150 (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                                                                    <  2:121299/150‑1 (MQ=255)
                                                                                  gACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGa                                                                    <  1:1259360/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                                                                                 >  1:2212695/1‑114 (MQ=255)
                                                                                                         tCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGca                                                                                 <  2:2212695/114‑1 (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                                                                    >  1:533815/1‑108 (MQ=255)
                                                                                                                            ttGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGa                                                                    <  2:533815/108‑1 (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)
                                                                                                                                                     |                                                                                                                                                    
CTTCATCATCGCGGGGGAAAGTATATCTAATGTATTTCTGAGCATCATCTTCATCGGCGCTGAGGGCCTGGACGGCTTGAACGACCTCGTCCAAACGAGATTCTATCGCTGGATCTAACGCGTATTGACGGGCCTGTCGCTCAAGAGCGTCGATCTTGATTTGAAGCGCCTCCACCTCAGCTGGAGACTCTGGTGTCTCAGGCTCCAAGTCATGAGGCACACTTGAAATAGATACATTCCTCAACTCGCCCCACGGGCCCCACTCACTCACAGGGTCTGTAAGATCTGCTTTGGCGTG  >  NC_002947/4008653‑4008950

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: