Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,203,029 C→A 7.1% L430M (CTG→ATG) ‡ PP_4583 → putative Peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,203,0290CA 7.1% 309.1 / 9.1 99L430M (CTG→ATG) ‡PP_4583putative Peptidase
Reads supporting (aligned to +/- strand):  ref base C (40/52);  new base A (4/3);  total (44/55)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

AAGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCTTTGC  >  NC_002947/5202883‑5203177
                                                                                                                                                  |                                                                                                                                                    
aaGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTgg                                                                                                                                                   >  1:53905/1‑150 (MQ=255)
aaGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTgg                                                                                                                                                   <  1:1620385/150‑1 (MQ=255)
    gcTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGcc                                                                                                                                               >  2:1824522/1‑150 (MQ=255)
      tCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGa                                                                                                                                             >  1:699147/1‑150 (MQ=255)
      tCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGACCACGAACGCGCACATGGTGGTCAAGGAAACCCCGGTGCTCGGGGCGTTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGa                                                                                                                                             >  2:1655985/1‑150 (MQ=255)
       cAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCAGGGTGCCgtc                                                                                                                                            <  2:1643520/150‑3 (MQ=255)
       cAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCAGGGTGCCgtc                                                                                                                                            <  2:460966/150‑3 (MQ=255)
        aaCCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATc                                                                                                                                           >  1:942776/1‑150 (MQ=255)
        aaCCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATc                                                                                                                                           <  2:834527/150‑1 (MQ=255)
         aCCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCAAGTGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCa                                                                                                                                          <  1:1041157/150‑1 (MQ=255)
            gCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGcc                                                                                                                                       <  2:132418/150‑1 (MQ=255)
            gCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGcc                                                                                                                                       <  2:130017/150‑1 (MQ=255)
             cccGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCt                                                                                                                                      >  2:1448746/1‑150 (MQ=255)
                   ccAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAcg                                                                                                                                >  1:1345539/1‑150 (MQ=255)
                   ccAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCAGGGTGCCgtc                                                                                                                                            >  2:204276/1‑136 (MQ=255)
                   ccAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCAGGGTGCCgtc                                                                                                                                            <  1:204276/138‑3 (MQ=255)
                    cAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAcgc                                                                                                                               >  1:1557661/1‑150 (MQ=255)
                           gcCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGACACCCTGGTGCCGATCAGCCTGGTACgccgccgcc                                                                                                                        >  2:482677/1‑150 (MQ=255)
                            ccAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCa                                                                                                                       <  2:24711/150‑1 (MQ=255)
                            ccAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCa                                                                                                                       <  2:1845013/150‑1 (MQ=255)
                             cAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAg                                                                                                                      >  2:1576185/1‑150 (MQ=255)
                                cTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGAc                                                                                                                   <  2:1465869/150‑1 (MQ=255)
                                  gACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCt                                                                                                                 <  2:1364030/150‑1 (MQ=255)
                                    cGTTGCCCACGGGCGCCCAGGTGGTGCTCCAGCACAAGCCGGTGCTCGGTGCATTCGAGGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTgg                                                                                                               <  2:1893392/150‑1 (MQ=255)
                                    cGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAACCCCCGGTGCTCCTTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTgg                                                                                                               <  1:206090/150‑1 (MQ=255)
                                    cGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTgg                                                                                                               <  2:62783/150‑1 (MQ=255)
                                    cGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTgg                                                                                                               <  2:1124864/150‑1 (MQ=255)
                                    cGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTgg                                                                                                               <  2:1376745/150‑1 (MQ=255)
                                      ttGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGc                                                                                                             >  2:31075/1‑150 (MQ=255)
                                        gCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCTGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCaa                                                                                                           <  2:1674743/150‑1 (MQ=255)
                                          ccACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGa                                                                                                         >  1:1813388/1‑150 (MQ=255)
                                          ccACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGa                                                                                                         >  2:1856712/1‑150 (MQ=255)
                                          ccACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGa                                                                                                         >  1:786254/1‑150 (MQ=255)
                                          ccACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGCCCAACCGCAGCGGCCGGCACCCTCGTCCCGATCAGCCTGGTACGCCGCCGCCAGCACCTGGGCAAGa                                                                                                         >  2:64179/1‑150 (MQ=255)
                                            aCGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGAcc                                                                                                       >  1:100552/1‑150 (MQ=255)
                                              gggCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGt                                                                                                     >  1:1036045/1‑150 (MQ=255)
                                              gggCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGt                                                                                                     <  1:677870/150‑1 (MQ=255)
                                              gggCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGt                                                                                                     <  1:1448746/150‑1 (MQ=255)
                                               ggCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTg                                                                                                    >  1:626282/1‑150 (MQ=255)
                                                gcgcCCAGGTGGTGCTCAAGCAAACCCAGGTGAACGGAACCCATGAGGCTGATGACTATGTCAGTGAGCGCATTGGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGc                                                                                                   <  1:1867298/150‑1 (MQ=255)
                                                 cgcCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGcc                                                                                                  >  1:642294/1‑150 (MQ=255)
                                                  gcCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCTATGACTCTGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCa                                                                                                 <  1:586437/150‑1 (MQ=255)
                                                    ccAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTggg                                                                                                              >  1:1101178/1‑135 (MQ=255)
                                                    ccAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTggg                                                                                                              <  2:1101178/135‑1 (MQ=255)
                                                    ccAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCAct                                                                                               >  2:1111950/1‑150 (MQ=255)
                                                        ggggtgCTCAAGCAAGCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtac                                                                                           <  2:950111/148‑1 (MQ=255)
                                                            tgCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtacctgt                                                                                       <  2:1801692/150‑1 (MQ=255)
                                                             gctacaGCAACCCCCTGCGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtacctgta                                                                                      <  2:987189/145‑1 (MQ=255)
                                                              cTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtacctgtac                                                                                     >  2:321944/1‑150 (MQ=255)
                                                                   gCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCt                                                                                                                                                     <  1:889585/81‑1 (MQ=255)
                                                                   gCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCt                                                                                                                                                     >  2:889585/1‑81 (MQ=255)
                                                                   gCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGGGGCACTGTACCTGTACGGCtg                                                                                >  1:1509930/1‑149 (MQ=255)
                                                                       aCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAg                                                                                                                      >  2:1925884/1‑108 (MQ=255)
                                                                       aCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAg                                                                                                                      <  1:1925884/108‑1 (MQ=255)
                                                                       aCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGt                                                                            >  1:850545/1‑150 (MQ=255)
                                                                        ccccGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTg                                                                           <  2:1345539/150‑1 (MQ=255)
                                                                         cccGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGc                                                                          >  2:145494/1‑150 (MQ=255)
                                                                         cccGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGc                                                                          <  2:1557661/150‑1 (MQ=255)
                                                                          ccGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCCACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGcc                                                                         <  2:1036045/150‑1 (MQ=255)
                                                                             gTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTAt                                                                      <  1:1111950/150‑1 (MQ=255)
                                                                             gTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTAt                                                                      <  2:626282/150‑1 (MQ=255)
                                                                             gTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTAt                                                                      <  2:1009092/150‑1 (MQ=255)
                                                                             gTGCTCGGTGCATTCGATGCCGATGACCATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTAt                                                                      <  2:1009034/150‑1 (MQ=255)
                                                                              tGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATg                                                                     <  1:188880/150‑1 (MQ=255)
                                                                                  cGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGa                                                                 >  2:1732430/1‑150 (MQ=255)
                                                                                    gTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGaaa                                                               <  1:1576185/150‑1 (MQ=255)
                                                                                    gTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGaaa                                                               <  1:198825/150‑1 (MQ=255)
                                                                                       cATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATg                                                                              <  2:1134639/132‑1 (MQ=255)
                                                                                       cATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATg                                                                              >  1:1134639/1‑132 (MQ=255)
                                                                                         ttCNATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTc                                                          >  2:1186516/1‑150 (MQ=255)
                                                                                         ttCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTc                                                          <  2:1813388/150‑1 (MQ=255)
                                                                                         ttCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTc                                                          <  2:100552/150‑1 (MQ=255)
                                                                                                ccGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACAGCGCCGCCAGGCCCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGTCTCGACCCCt                                                   >  1:891369/1‑150 (MQ=255)
                                                                                                        tATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTca                                           <  1:1186516/150‑1 (MQ=255)
                                                                                                        tATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTca                                           <  1:64179/150‑1 (MQ=255)
                                                                                                         aTGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTcac                                          <  2:541997/150‑1 (MQ=255)
                                                                                                          tGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTg                                                  <  1:375013/141‑1 (MQ=255)
                                                                                                          tGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTg                                                  >  2:375013/1‑141 (MQ=255)
                                                                                                            tCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATg                                       <  1:1824522/150‑1 (MQ=255)
                                                                                                            tCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATg                                       >  2:1480021/1‑150 (MQ=255)
                                                                                                              aGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCGTTTCCTTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATgcg                                     >  1:1368691/1‑150 (MQ=255)
                                                                                                               gTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATgcgc                                    <  2:1717623/150‑1 (MQ=255)
                                                                                                                 gAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATgcgcgc                                  <  2:786254/150‑1 (MQ=255)
                                                                                                                 gAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGAGACCAAGGAGATAGCCGCGACCCCTGGTTCTCACATgcccgc                                  >  2:152293/1‑150 (MQ=255)
                                                                                                                  aGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGcc                                 >  2:463127/1‑150 (MQ=255)
                                                                                                                    cgcCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGt                                                                                                     >  2:616281/1‑80 (MQ=255)
                                                                                                                    cgcCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGt                                                                                                     <  1:616281/80‑1 (MQ=255)
                                                                                                                     gcCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGa                              <  2:699147/150‑1 (MQ=255)
                                                                                                                       cTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGc                            >  2:242554/1‑150 (MQ=255)
                                                                                                                               aaCCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTCTGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGa                    <  2:62566/150‑1 (MQ=255)
                                                                                                                                       tgatcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCgcggcg            >  2:520310/5‑150 (MQ=255)
                                                                                                                                       tgatcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCgcggcg            >  2:525798/5‑150 (MQ=255)
                                                                                                                                         gACGGCACCCTGGTGCCGGTCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTg          <  2:254270/150‑1 (MQ=255)
                                                                                                                                         gACGGCACCCTGGTGCCGATCAGCCTGGTCCGCCGCCGCCCGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTg          <  2:1490262/150‑1 (MQ=255)
                                                                                                                                          tcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTgg         >  1:63970/2‑150 (MQ=255)
                                                                                                                                          actgCACCCTGATGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTgg         <  1:505670/147‑1 (MQ=255)
                                                                                                                                            ggCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGcc       >  1:1042743/1‑150 (MQ=255)
                                                                                                                                              cACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCtt     <  1:242554/150‑1 (MQ=255)
                                                                                                                                               aCCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCttt    >  1:1607332/1‑150 (MQ=255)
                                                                                                                                                 ccTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCTTTGc  <  2:642294/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCTTTGC  >  NC_002947/5202883‑5203177

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