Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,141,212 T→A 7.0% L44Q (CTG→CAG)  PP_3645 → L‑proline amide hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,141,2120TA 7.0% 309.3 / 10.8 101L44Q (CTG→CAG) PP_3645L‑proline amide hydrolase
Reads supporting (aligned to +/- strand):  ref base T (48/43);  new base A (5/2);  total (54/47)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 116 total aligned reads displayed.

CTTCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCGGC  >  NC_002947/4141064‑4141360
                                                                                                                                                    |                                                                                                                                                       
cTTCAAGGAGCGCCCTCCACGTTTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGAGGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGGAATGTGGGGCGGTCAGTTCGTtggcctg                                                                                                                                                        <  2:1582203/150‑1 (MQ=255)
 ttCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctgg                                                                                                                                                       >  2:1092918/1‑150 (MQ=255)
 ttCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctgg                                                                                                                                                       <  1:564657/150‑1 (MQ=255)
      ggAGCGCCCTCCATGATTTTGCCTTTCTCTAAACGCCCGCCTGCCGCGACCTGCTCTCTGGTTGTCGGCCAGTGCCCACCCGTGGCGCTGTTCCACGCTGTGGGCCTCTACACGGAAAGGGGGGGCCGTTGGTTCGTTGGCAGGGCCaaa                                                                                                                                                  <  2:1703624/150‑2 (MQ=255)
      ggAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCCGGCCAAc                                                                                                                                                  >  1:654062/1‑150 (MQ=255)
         gcgcCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑Ac                                                                                                                                              <  2:772596/150‑1 (MQ=255)
               tCCATGATTGAGCCTGTCGCTGAACGCACACCTGCCAGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTCATCCACGGCGTGGGCCTGAACAAGGAAACCTATCGGGGTCAATTCGGTGGCCTAGCCAACG‑ACT‑aagca                                                                                                                                       >  1:738183/1‑146 (MQ=255)
               tCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACAAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGGGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACcgc                                                                                                                                       <  1:489161/150‑1 (MQ=255)
                ccATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACcgcg                                                                                                                                      <  2:1427917/150‑1 (MQ=255)
                 cATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCAGGCCAACG‑aa                                                                                                                                              >  2:1775555/1‑141 (MQ=255)
                 cATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCAGGCCAACG‑aa                                                                                                                                              <  1:1775555/142‑2 (MQ=255)
                    gATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCAt                                                                                                                                  <  1:1347116/150‑1 (MQ=255)
                         aGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCt                                                                                                                             <  2:2236605/150‑1 (MQ=255)
                         aGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCt                                                                                                                             <  2:2033404/150‑1 (MQ=255)
                           ccTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑Ac                                                                                                                                          <  1:102001/135‑1 (MQ=255)
                           ccTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTAc                                                                                                                           >  1:1187344/1‑150 (MQ=255)
                           ccTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGACCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑Ac                                                                                                                                          >  2:102001/1‑135 (MQ=255)
                                cGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACAt                                                                                                                      <  2:1187344/150‑1 (MQ=255)
                                 gCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATg                                                                                                                     >  1:1085686/1‑150 (MQ=255)
                                  cTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACG‑ACCGCGTCATCGCCTACGACATGc                                                                                                                    >  1:1811436/1‑150 (MQ=255)
                                       cGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGc                                                                                                              <  2:1237892/150‑1 (MQ=255)
                                            acCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑AGCGCGTGATCGCCTACGACATGCT‑CGGCCACgg                                                                                                         >  1:1964547/1‑150 (MQ=255)
                                              cTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctggc                                                                                                                                                      >  1:853556/1‑106 (MQ=255)
                                              cTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTtggcctggc                                                                                                                                                      <  2:853556/106‑1 (MQ=255)
                                              cTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑Cg                                                                                                                >  2:1876497/1‑141 (MQ=255)
                                              cTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑Cg                                                                                                                <  1:1876497/141‑1 (MQ=255)
                                                      aCCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCcgc                                                                                               <  1:1092918/150‑1 (MQ=255)
                                                      aCCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCcgc                                                                                               <  1:1421216/150‑1 (MQ=255)
                                                        cAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGt                                                                                             >  2:350334/1‑150 (MQ=255)
                                                                  gATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGgccgcc                                                                                   <  2:1914984/150‑1 (MQ=255)
                                                                    tGTCGGCCAGGGCCAACCCGGGGTGCTGATCCCCGGCGTGGGCCTCAACAAGGAAATGTGGGCCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGa                                                                                 <  2:210379/150‑1 (MQ=255)
                                                                        ggCCAGGTCCAACACGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGacaca                                                                             <  1:498786/150‑1 (MQ=255)
                                                                              ggCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGa                                                                                                                         <  2:1511851/101‑1 (MQ=255)
                                                                              ggCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGa                                                                                                                         >  1:1511851/1‑101 (MQ=255)
                                                                                 cAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCcgc                                                                                               <  2:140802/123‑1 (MQ=255)
                                                                                 cAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCcgc                                                                                               >  1:140802/1‑123 (MQ=255)
                                                                                 cAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGaa                                                                    <  1:350334/150‑1 (MQ=255)
                                                                                 cAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGaa                                                                    >  2:846591/1‑150 (MQ=255)
                                                                                  aaCCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAg                                                                   >  2:1354655/1‑150 (MQ=255)
                                                                                   aCCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAgg                                                                  >  2:909625/1‑150 (MQ=255)
                                                                                   aCCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCCGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAgg                                                                  <  2:654062/150‑1 (MQ=255)
                                                                                    cccGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGc                                                                 <  2:1811436/150‑1 (MQ=255)
                                                                                     ccGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCt                                                                        <  1:88412/142‑1 (MQ=255)
                                                                                     ccGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCt                                                                        >  2:88412/1‑142 (MQ=255)
                                                                                     ccGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCt                                                                >  2:888210/1‑150 (MQ=255)
                                                                                     ccGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCt                                                                >  1:965384/1‑150 (MQ=255)
                                                                                     ccGTGGTGCTGATCCACGGCGTGGGCCTCATCAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCt                                                                <  1:1811947/150‑1 (MQ=255)
                                                                                     ccGTGGTGCTGATCCACGGCGTGGGCCCGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCt                                                                >  2:118319/1‑150 (MQ=255)
                                                                                      cGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTa                                                               >  1:329595/1‑150 (MQ=255)
                                                                                        tggtgCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGCCATGGT‑CGCGTGTTGTCAGAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATg                                                             >  2:1269131/1‑150 (MQ=255)
                                                                                            gCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCcg                                                                                                >  2:1653985/1‑111 (MQ=255)
                                                                                            gCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCcg                                                                                                <  1:1653985/111‑1 (MQ=255)
                                                                                             cTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGAc                                                        <  2:1174432/150‑1 (MQ=255)
                                                                                               gATCCACGGCGTGGGCCTGAACAAGGAAAAGTGGGGCGGTCAGTCCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCTAAGAGCGTCGTGTAGGGC‑CGGCGCGCCCGCCGACACACCGCTGGCAGGCTATGCCGACCa                                                      >  2:674499/1‑150 (MQ=255)
                                                                                                aTCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGcc                                                           >  1:1393988/1‑144 (MQ=255)
                                                                                                aTCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGcc                                                           <  2:1393988/144‑1 (MQ=255)
                                                                                                 tCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGc                                                    <  1:846591/150‑1 (MQ=255)
                                                                                                   cACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGc                                                                                                                    <  1:1613919/85‑1 (MQ=255)
                                                                                                   cACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGc                                                                                                                    >  2:1613919/1‑85 (MQ=255)
                                                                                                   cACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTa                                                               >  2:2077025/1‑137 (MQ=255)
                                                                                                   cACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTa                                                               <  1:2077025/137‑1 (MQ=255)
                                                                                                    aCGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTAt                                                              <  1:1958046/137‑1 (MQ=255)
                                                                                                    aCGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTAt                                                              >  2:1958046/1‑137 (MQ=255)
                                                                                                    aCGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCg                                                 >  1:2314936/1‑150 (MQ=255)
                                                                                                          tGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAAc                                           <  1:1988477/150‑1 (MQ=255)
                                                                                                              ccTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCt                                       >  1:1002005/1‑150 (MQ=255)
                                                                                                               cTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCtt                                      <  2:161081/150‑1 (MQ=255)
                                                                                                                tGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTg                                     >  2:2141537/1‑150 (MQ=255)
                                                                                                                   aCAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGacc                                  >  1:752189/1‑150 (MQ=255)
                                                                                                                     aaGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGaccac                                >  1:930709/1‑150 (MQ=255)
                                                                                                                     aaGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGaccac                                >  1:946698/1‑150 (MQ=255)
                                                                                                                     aaGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGaccac                                >  2:2271154/1‑150 (MQ=255)
                                                                                                                     aaGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGaccac                                >  1:256328/1‑150 (MQ=255)
                                                                                                                     aaGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCCCCGAACTGCTTGaccac                                >  2:69970/1‑150 (MQ=255)
                                                                                                                       ggAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCt                              <  1:888210/150‑1 (MQ=255)
                                                                                                                       ggAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCt                              <  1:909625/150‑1 (MQ=255)
                                                                                                                        gAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCATAGGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTg                             >  2:1791621/1‑150 (MQ=255)
                                                                                                                            tgtgGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCaaa                         <  1:557122/150‑1 (MQ=255)
                                                                                                                                 ggCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCCCGGTCAAAGCCGCGTACCGGCGGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGca                    >  2:2099819/1‑150 (MQ=255)
                                                                                                                                  gCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGc                                                                 <  2:1692951/104‑1 (MQ=255)
                                                                                                                                  gCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGc                                                                 >  1:1692951/1‑104 (MQ=255)
                                                                                                                                  gCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGcac                   >  2:1623627/1‑150 (MQ=255)
                                                                                                                                        aGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTgt                                        >  2:140510/1‑122 (MQ=255)
                                                                                                                                        aGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGc                                        <  1:140510/123‑1 (MQ=255)
                                                                                                                                          ttCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCAcc           >  2:158513/1‑150 (MQ=255)
                                                                                                                                          ttCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCAcc           >  2:1908777/1‑150 (MQ=255)
                                                                                                                                           tCGTTGGCCAGGCCAACGAACTGACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCAc            >  1:558488/1‑150 (MQ=255)
                                                                                                                                           tCGTTGGCCAGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCg          >  2:1618880/1‑150 (MQ=255)
                                                                                                                                           tCGTTGGCCAGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGGCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCg          >  2:1974734/1‑150 (MQ=255)
                                                                                                                                             gTTGGCCAGGCCAACGAACTGACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGaa                                                                    >  1:2092398/1‑92 (MQ=38)
                                                                                                                                             gTTGGCCAGGCCAACGAACTGACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGaa                                                                    <  2:2092398/92‑1 (MQ=38)
                                                                                                                                                 gcctggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CTGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGcac                   <  2:400900/135‑1 (MQ=255)
                                                                                                                                                 gcctggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCg    <  1:164808/150‑1 (MQ=255)
                                                                                                                                                 gcctggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCg    >  1:2338520/1‑150 (MQ=255)
                                                                                                                                                 gccgggccAACG‑ACT‑ACCACGTCAGCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCCCCGACACACCGCTGGAAGGCTATGCCGACAAGCTCGCCGAACTGTTTGCACACCTGCAAATTGCACGAGCCACCGTGATCg    <  2:1002005/150‑1 (MQ=255)
                                                                                                                                                   ctggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAAt                        <  2:1315522/128‑1 (MQ=255)
                                                                                                                                                   ctggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAAt                        >  1:1315522/1‑128 (MQ=255)
                                                                                                                                                   ctggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACaa                 >  1:1787000/1‑135 (MQ=255)
                                                                                                                                                   ctggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACaa                 <  2:1787000/135‑1 (MQ=255)
                                                                                                                                                   ctggccAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCGGc  <  1:1184181/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
CTTCAAGGAGCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACG‑ACT‑ACCGCGTCATCGCCTACGACATGCT‑CGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACCGTGATCGGC  >  NC_002947/4141064‑4141360

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