Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,018,187 A→T 28.5% intergenic (+63/+68) prfC → / ← potF‑I peptide chain release factor 3/putrescine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,018,1870AT28.5% 267.3 / 134.5 173intergenic (+63/+68)prfC/potF‑Ipeptide chain release factor 3/putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (51/68);  new base T (21/28);  total (73/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.79e-01

Only 100 of 177 total aligned reads displayed.

TACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGC  >  NC_002947/1018038‑1018336
                                                                                                                                                     |                                                                                                                                                      
tacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCa                                                                                                                                                        <  1:1486952/150‑1 (MQ=255)
tacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCa                                                                                                                                                        >  1:1842819/1‑150 (MQ=255)
  cctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGa                                                                                                                                                      <  1:346913/150‑1 (MQ=255)
       cgcCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCATTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTgg                                                                                                                                                 <  1:1708150/150‑1 (MQ=255)
        gcCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGa                                                                                                                                                <  2:92997/150‑1 (MQ=255)
        gcCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGa                                                                                                                                                <  1:4993502/150‑1 (MQ=255)
        gcCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGa                                                                                                                                                <  1:4993469/150‑1 (MQ=255)
        gcCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGa                                                                                                                                                >  1:2834348/1‑150 (MQ=255)
         ccGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAg                                                                                                                                               >  2:2899859/1‑150 (MQ=255)
           gACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGt                                                                                                                                             <  1:15469/150‑1 (MQ=255)
            aCCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTc                                                                                                                                            <  2:1033542/150‑1 (MQ=255)
              ccGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAt                                                                                                                                          <  1:2183765/150‑1 (MQ=255)
              ccGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAt                                                                                                                                          <  2:2333179/150‑1 (MQ=255)
               cGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                                         <  2:5549394/150‑1 (MQ=255)
               cGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACACGAATCACCAGACCTGGAGGTCAtg                                                                                                                                         >  2:2659555/1‑150 (MQ=255)
                gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTctggt                                                                                                                                        <  2:2258225/150‑5 (MQ=37)
                gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTctggt                                                                                                                                        <  2:3533143/150‑5 (MQ=37)
                gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                        <  1:1011864/150‑1 (MQ=255)
                gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                        <  2:2834348/150‑1 (MQ=255)
                  gTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGa                                                                                                                                      <  1:2888036/150‑1 (MQ=255)
                   tCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGACTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGAt                                                                                                                                     <  2:1785386/150‑1 (MQ=255)
                    cAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATa                                                                                                                                    <  1:1647809/150‑1 (MQ=255)
                     aaCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGAttc                                                                                                                                   <  2:1918028/150‑3 (MQ=25)
                     aaCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                                                                                   <  1:5325027/150‑1 (MQ=255)
                     aaCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCCGCCTAGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGAttc                                                                                                                                   <  1:4436470/150‑3 (MQ=18)
                      aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTcc                                                                                                                                  <  2:2462869/150‑1 (MQ=25)
                      aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAAGCACCACACCGGGAGGTCCTGTGAGAcc                                                                                                                                  >  2:2051232/1‑150 (MQ=38)
                       ccTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCt                                                                                                                                 <  2:5537064/150‑1 (MQ=255)
                       ccTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCt                                                                                                                                 <  2:3160485/150‑1 (MQ=255)
                       ccTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCCGCCTAGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCt                                                                                                                                 <  1:2646142/150‑1 (MQ=18)
                        cTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtg                                                                                                                                <  2:272324/150‑1 (MQ=25)
                        cTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtg                                                                                                                                >  2:418298/1‑150 (MQ=255)
                         tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGTGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                               <  1:5051229/150‑1 (MQ=255)
                         tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgt                                                                                                                               <  1:2293271/150‑1 (MQ=25)
                         tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                               <  2:4632953/150‑1 (MQ=255)
                         tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                               <  2:3735726/150‑1 (MQ=255)
                         tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                               <  1:518671/150‑1 (MQ=255)
                          ggCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgtg                                                                                                                              <  2:2546380/150‑1 (MQ=25)
                          ggCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgtg                                                                                                                              <  2:2546685/150‑1 (MQ=25)
                           gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTgg                                                                                                                             <  2:1057506/150‑1 (MQ=25)
                           gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTgg                                                                                                                             <  1:1980190/150‑1 (MQ=25)
                           gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTgg                                                                                                                             <  2:4353137/150‑1 (MQ=255)
                           gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATGCAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTgg                                                                                                                             <  2:3336599/150‑1 (MQ=255)
                             gcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGa                                                                                                                           <  1:3780249/150‑1 (MQ=21)
                             gcTGATGGAAGAGCGCTGGCCGGACATCAAGGTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCCGCCTAGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGa                                                                                                                           <  1:660286/150‑1 (MQ=17)
                                 aTGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCgg                                                                                                                       <  2:5297278/150‑1 (MQ=255)
                                  tGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCggg                                                                                                                      <  2:1297509/150‑1 (MQ=255)
                                  tGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCggg                                                                                                                      <  2:1297531/150‑1 (MQ=255)
                                       gagCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGc                                                                                                                 <  2:3774318/150‑1 (MQ=255)
                                         gcgcTGGCCGGACATCCAGTTCCGCGCGACCCGCGAGCACCATTCAGGCCTGAACTCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGccc                                                                                                               <  2:3466669/150‑1 (MQ=17)
                                                  ggACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTgg                                                                                                                             <  1:438370/127‑1 (MQ=18)
                                                  ggACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTgg                                                                                                                             >  2:438370/1‑127 (MQ=18)
                                                       tCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCAGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCtg                                                                                                >  2:601086/1‑149 (MQ=255)
                                                        cAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGa                                                                                                                                                      <  2:4670180/96‑1 (MQ=255)
                                                             tCCGCGCTCCCCGCGAGCACAATTACGGCCTTACTTGGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCCCCCACCGGAATCACCAGACCTGGCGGTCATGTGATACCGGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAggcggc                                                                                          <  1:4690175/150‑1 (MQ=25)
                                                                  gcTACCCGCGAGCACCATTACGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAAGCACCAGACCTGGAGGACATGTGTTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCCTAGCGGcccttc                                                                                     <  2:1706672/150‑4 (MQ=25)
                                                                    tACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAACAGACCTGGAGGTGATGTGATAGCTTTGGGTGCGGGCGTGCCCGCGAATGGCCGCGTTGGGGGCCCAAtt                                                                                    >  1:913165/1‑150 (MQ=25)
                                                                              cACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCgg                                                                         <  2:4563197/150‑1 (MQ=255)
                                                                                 cATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGAAACACCACACCAGGAGGTCAAGTGATACCGGTGGGAGCGGGCGTGCCCGCCAAGGGCCGCTCGGCGGCCCCAATTCAAATTCCGGCtt                                                                      >  2:5342492/1‑150 (MQ=37)
                                                                                 cATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCCCAGGAATCACCACACCTGGAGGTCATGTGATACCTGTGGGCGCGGGCGTGCCCCCGAAGGGCCGCGAGGCGGCCCCCTTTAAATTTCCGGCtg                                                                      >  2:1502310/1‑149 (MQ=37)
                                                                                        gCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTACCACAGGAAAAACCAGACCTGGAGGTCAAGTGATAACGGTGGGGGCGGGCGTGGCCGCGGAGGGCCGCTAGGCGGCCCCCATTCAAATTACGGCTTTGATaat                                                               >  2:5194834/1‑148 (MQ=25)
                                                                                           tGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                        <  1:3908229/75‑1 (MQ=255)
                                                                                           tGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                        >  2:3908229/1‑75 (MQ=255)
                                                                                                    ggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGCATAACCAGACCAGGAGGTCATGTGATACCGGTGGGAGCGGGGGGGCCCGCGAAAGGCCGCGAGGCGGGCCCCTGTCCAATACCCGCGTTGAGAAGGGGCCAGAttgg                                                   >  2:683491/1‑148 (MQ=16)
                                                                                                        gctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGTAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTa                                               <  1:3675642/150‑1 (MQ=255)
                                                                                                        gctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGTTGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTa                                               <  1:2299507/150‑1 (MQ=255)
                                                                                                        gctGTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGTGTTCAGGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTa                                               <  1:4551301/150‑1 (MQ=255)
                                                                                                         ctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTTATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTAc                                              <  1:1833947/150‑1 (MQ=255)
                                                                                                         ctTGGCCGCCCTCCGCGGGCCCGCCCGCTCCCACCGCAACCCCCCGCCGTGGCGGTCCTGGCAGCCCTGTGGGCGGGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTAc                                              <  1:1226168/150‑1 (MQ=255)
                                                                                                              cAGCCCTGCGCGGGCACGCCCGCTCCCACAGTAACCACCAGACGTGGTGGTCATTTTATCCGTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATa                                         <  1:2803191/150‑1 (MQ=38)
                                                                                                                   cTTCGCGGGCCCCCCCGCCCCCCCCGTAAACCCCCGCCCTGGAGGTCCGGTGATACCTGTGGGCGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTc                                    <  1:4316833/150‑1 (MQ=255)
                                                                                                                     ccgcgGGCCCGCCCGCCCCCCCCGGAACCCCCCGACGTGGAGGTCATTTGTTACCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGa                                  <  1:2659555/149‑1 (MQ=255)
                                                                                                                             aCGCCCGCTCCCACCGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGcc                          <  2:3871169/150‑1 (MQ=255)
                                                                                                                                 ccGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCa                      >  1:5406653/1‑150 (MQ=255)
                                                                                                                                  cGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCAt                     <  2:5220368/150‑1 (MQ=255)
                                                                                                                                       ccACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt                >  1:4600888/1‑150 (MQ=25)
                                                                                                                                       ccACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt                >  1:387667/1‑150 (MQ=25)
                                                                                                                                       ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt                >  2:2131784/1‑150 (MQ=255)
                                                                                                                                       ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt                >  2:3829232/1‑150 (MQ=255)
                                                                                                                                       ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTGCAAGTTCCGGCGTTGATCCGGGTCCAGATTCGCGTACGAATACGGTCGAGCGTCGCCGGCATCGccg                >  1:2167451/1‑149 (MQ=38)
                                                                                                                                        cacaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc               >  1:3716393/1‑150 (MQ=25)
                                                                                                                                        cacaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCAGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc               >  2:2817109/1‑150 (MQ=21)
                                                                                                                                         acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc              >  1:1748147/1‑150 (MQ=25)
                                                                                                                                         acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc              >  1:3555954/1‑150 (MQ=25)
                                                                                                                                         acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc              >  1:5012454/1‑150 (MQ=255)
                                                                                                                                         acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc              >  1:4417526/1‑150 (MQ=255)
                                                                                                                                         acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc              >  1:443137/1‑150 (MQ=255)
                                                                                                                                         acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc              >  1:4575894/1‑150 (MQ=255)
                                                                                                                                          caGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCa             >  1:4806568/1‑150 (MQ=255)
                                                                                                                                          caGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCa             >  2:5030268/1‑150 (MQ=255)
                                                                                                                                           aGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAg            >  2:5203976/1‑150 (MQ=25)
                                                                                                                                           aGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTACAg            >  1:510072/1‑150 (MQ=25)
                                                                                                                                           aGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAg            >  2:1377345/1‑150 (MQ=255)
                                                                                                                                             gtaTCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgcg          >  1:4114231/3‑150 (MQ=25)
                                                                                                                                              aaTCCCCTGCCCTGGTGGTCTGGTGAGACGTGGGGGCGGGGGGGGGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgcgc         <  1:2899859/150‑1 (MQ=255)
                                                                                                                                               aTCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCa        >  2:434474/1‑150 (MQ=255)
                                                                                                                                                tCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCAt       >  1:730733/1‑150 (MQ=255)
                                                                                                                                                 cACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa      <  2:4430013/150‑1 (MQ=255)
                                                                                                                                                 cACCAGACCTGGAGGTCATGCGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa      <  2:180807/150‑1 (MQ=255)
                                                                                                                                                  aCCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATAc     >  1:1293087/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGC  >  NC_002947/1018038‑1018336

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