Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,018,187 A→T 27.6% 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,1870AT27.6% 706.6 / inf 432intergenic (+63/+68)prfC/potF‑Ipeptide chain release factor 3/putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (140/158);  new base T (43/75);  total (190/242)
Fisher's exact test for biased strand distribution p-value = 6.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-02

Only 100 of 449 total aligned reads displayed.

TACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGC  >  NC_002947/1018038‑1018336
                                                                                                                                                     |                                                                                                                                                      
tacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCa                                                                                                                                                        <  1:331652/150‑1 (MQ=255)
tacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCa                                                                                                                                                        <  2:2892660/150‑1 (MQ=255)
tacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCa                                                                                                                                                        <  2:731489/150‑1 (MQ=255)
tacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACACGCATCACCa                                                                                                                                                        >  1:6101378/1‑150 (MQ=255)
    tGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGAcc                                                                                                                                                    >  2:65868/1‑150 (MQ=255)
    tGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCCCAGCCAACAACAGAcc                                                                                                                                                    >  1:1821939/1‑150 (MQ=39)
     ggCGCCGACCCGTGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCt                                                                                                                                                   <  2:3741760/150‑1 (MQ=255)
       cgcCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCtcc                                                                                                                                                 <  1:7061370/150‑3 (MQ=255)
         ccGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAg                                                                                                                                               <  1:5896311/150‑1 (MQ=255)
            aCCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATGGGGCCGCCTAGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTct                                                                                                                                           <  1:1641896/150‑2 (MQ=21)
              ccGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATAACCAGACCTGGAGGGCAt                                                                                                                                          >  2:24960/1‑150 (MQ=255)
               cGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                                         <  1:6688818/150‑1 (MQ=255)
               cGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                                         <  1:4584468/150‑1 (MQ=255)
               cGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                                         <  2:7887077/150‑1 (MQ=255)
                gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTctggt                                                                                                                                        <  2:576024/150‑5 (MQ=37)
                gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                        <  1:5104184/150‑1 (MQ=255)
                gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                        <  1:1099059/150‑1 (MQ=255)
                gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                        <  1:4260515/150‑1 (MQ=255)
                gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                                        <  1:3934675/150‑1 (MQ=255)
                 ggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACATGACCTCCAGGTctggtg                                                                                                                                       <  1:125751/150‑6 (MQ=37)
                  gTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGa                                                                                                                                      <  1:6585965/150‑1 (MQ=255)
                  gTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGa                                                                                                                                      <  1:2669038/150‑1 (MQ=255)
                   tCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGAt                                                                                                                                     <  2:4860746/150‑1 (MQ=255)
                   tCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGAt                                                                                                                                     <  2:7980466/150‑1 (MQ=255)
                   tCAACCTGGCGCTGATGTAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATa                                                                                                                                    <  1:4067417/150‑1 (MQ=255)
                     aaCCTGGCGCTGATGGAAGAGCGTTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                                                                                   <  1:2995626/150‑1 (MQ=255)
                     aaCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGAttc                                                                                                                                   <  2:3916818/150‑3 (MQ=25)
                     aaCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                                                                                   <  2:1504349/150‑1 (MQ=255)
                       ccTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCt                                                                                                                                 <  1:3709989/150‑1 (MQ=25)
                       ccTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTCGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCt                                                                                                                                 <  1:5506259/150‑1 (MQ=255)
                        cTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATGGGGCCGCCTAGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgt                                                                                                                               <  2:6407916/150‑1 (MQ=17)
                         tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgt                                                                                                                               <  1:7906652/150‑1 (MQ=25)
                         tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgt                                                                                                                               <  2:8470959/150‑1 (MQ=25)
                         tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                               <  2:7766700/150‑1 (MQ=255)
                         tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                               <  2:6743035/150‑1 (MQ=255)
                         tGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                               <  1:2073458/150‑1 (MQ=255)
                         tGGCGCTGATGGAAGAGCGCTGGCCCGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgt                                                                                                                               <  1:1436770/150‑1 (MQ=25)
                          ggCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgtg                                                                                                                              <  2:716954/150‑1 (MQ=25)
                           gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTgg                                                                                                                             <  1:1973933/150‑1 (MQ=255)
                           gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTgg                                                                                                                             <  1:2218501/150‑1 (MQ=255)
                           gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATGGGGCCGCCTAGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTggg                                                                                                                            <  2:5079623/150‑1 (MQ=17)
                              cTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAg                                                                                                                          <  1:3566214/150‑1 (MQ=255)
                                 aTGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCgg                                                                                                                       <  1:8294575/150‑1 (MQ=21)
                                   ggAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCAGAAATCGGGGCTGCTATGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGCAAGCACCAGACCAGGAGGGCCTGTCATACCGCTGGGAGCGGGc                                                                                                                     >  2:1891026/1‑150 (MQ=37)
                                       gagCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCTTGTGATACCTGTGGGAGCGGGCGTGc                                                                                                                 <  2:2929270/150‑1 (MQ=255)
                                            cTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGGCCGCCCGCTCCCACAGGAATCACCACAACAGGAGGTCAAGTGATAACTGTGGGGGCGGGGGTGCCcgcc                                                                                                            >  1:2310400/1‑149 (MQ=38)
                                              ggCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAAGGGGGGCGGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACGTGTGGGCGCGGGCGTGCCCGCGaa                                                                                                          <  2:6147376/150‑1 (MQ=37)
                                                   gACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAACTCGGGGCTGCTTTGCAGCCCTCCGCGGGCACGCCCGCTCCCACCGGTATCACATGCCCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGct                                                                                                     <  2:8655487/150‑2 (MQ=17)
                                                        cAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGCAACACCCGACCGGGGGGTCAGGTGGTACCGGGGGGGGCGGGCGGGGCCGCGAAGGGCCGCGAg                                                                                                >  2:1138258/1‑150 (MQ=25)
                                                                    tACCCGCGAGCACCATTAAGGCCTTAAATCGGGGCTGCTTTGAAGCCCTTCGCGGGCACGCCCACTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCAttt                                                                                   <  2:3008943/150‑1 (MQ=255)
                                                                      cccGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCt                                                                                                                                 >  1:8558147/1‑103 (MQ=14)
                                                                       ccGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCCTGTGATACCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCaa                                                                                >  2:5530191/1‑150 (MQ=255)
                                                                        cgcgAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAGGCGGCCCCATGTCaag                                                                               >  2:4930536/1‑149 (MQ=255)
                                                                              cACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGAttc                                                                                                                                   >  1:7410169/1‑91 (MQ=14)
                                                                                  aTTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCCGTTCAATTTCCGGCttt                                                                     >  2:7514232/1‑150 (MQ=255)
                                                                                         ccTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACCGGAATCACCAGACCTGGTGGTCTTGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTg                                                              <  2:8194306/150‑1 (MQ=255)
                                                                                           tGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGt                                                            <  2:3584742/150‑1 (MQ=255)
                                                                                           tGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCGCCCCCACGAATCACCAGACCGGGCGGTCATGGGGAACCGGGGGGGGGGGGGGGGCCCGCGCGGGGGCGCGAGGCGGGCCCCGTTCAAGTTCCGGGTGTGTTACTGGt                                                            >  2:809573/1‑150 (MQ=16)
                                                                                              gcgaagggCTGCAAAGCAGCCCTTCGCGGGCACGCCCCCTCCCACACGAATCACCAGACCGGGAGGGCATGTGTTACCGGGGGGGGCGGGGGTGCCCGCGCAGGGCCGCGAGGCGGCCCCATTTCAATTTCCGGCTGTGTTACTGGTCCa                                                         >  2:1778167/6‑150 (MQ=255)
                                                                                                  agggCTGCAAAGCAGCCCTTCGCGGGCACGCCCGCGCCCACAGGAATCACCAGACCGGGGGGTCAAGTGATACCGGGGGGGGCGGGCGTGCCCGCGCAGGGCCGCGAGGCGGCCCCAATTCAATTTCCGGCTTTGGTACTGGTCCAGAtt                                                     >  1:8032325/2‑150 (MQ=25)
                                                                                                     gctgctTTGCCGCCCTGCGCGGGCACGCCCGCTCCCACAGCAATCACCCGACCTGGAGGTCATGTTATACCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTcgc                                                  <  2:308668/150‑1 (MQ=37)
                                                                                                     gctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCCCAGGAATCACCAGACCTGGAGGTCTTGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTcgc                                                  <  2:5286266/150‑1 (MQ=255)
                                                                                                     gctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTcgc                                                  <  2:2051716/150‑1 (MQ=255)
                                                                                                             gCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATAACCAGACCTGGAGGTCATGTGATACCGGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATGTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAAt                                          >  2:6049453/1‑150 (MQ=255)
                                                                                                              cAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGGGGTCATGTGATACCGGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATa                                         >  2:4368040/1‑150 (MQ=255)
                                                                                                                   cTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTc                                    <  2:5216022/150‑1 (MQ=255)
                                                                                                                          ggCACGCCCGCTCCCACCGGAATCACCAGACCTGGAGGTCATGTGTTCCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTc                             <  2:4667008/150‑1 (MQ=255)
                                                                                                                          ggCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTc                             <  1:6502232/150‑1 (MQ=255)
                                                                                                                             aCGCCCGCTCCCACCGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGcc                          <  2:5035913/150‑1 (MQ=255)
                                                                                                                               gCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCgg                        >  2:5984892/1‑150 (MQ=21)
                                                                                                                                cccGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGc                       >  1:7218052/1‑150 (MQ=255)
                                                                                                                                   gCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATc                    >  1:8772819/1‑150 (MQ=255)
                                                                                                                                     tCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGc                  >  2:2632419/1‑150 (MQ=255)
                                                                                                                                      cccACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCGATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGc                  >  2:2690906/1‑150 (MQ=14)
                                                                                                                                       ccACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt                >  2:4150087/1‑150 (MQ=25)
                                                                                                                                       ccACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt                >  2:2286358/1‑150 (MQ=25)
                                                                                                                                       ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt                >  1:8425950/1‑150 (MQ=255)
                                                                                                                                        cacaGGTATCACATGACCTCCAGGTCTGGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc               >  1:6274798/1‑150 (MQ=25)
                                                                                                                                        cacaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc               >  2:1653238/1‑150 (MQ=255)
                                                                                                                                         acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc              >  2:5464373/1‑150 (MQ=25)
                                                                                                                                         acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc              >  2:1606968/1‑150 (MQ=25)
                                                                                                                                         acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc              >  1:152996/1‑150 (MQ=255)
                                                                                                                                         acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc              >  2:1336140/1‑150 (MQ=255)
                                                                                                                                         acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTcc              >  1:3242171/1‑150 (MQ=255)
                                                                                                                                          caGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCa             <  2:1393612/150‑1 (MQ=255)
                                                                                                                                          caGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCa             >  1:2817849/1‑150 (MQ=255)
                                                                                                                                          caGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCa             >  2:8262998/1‑150 (MQ=255)
                                                                                                                                           aGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAg            >  1:1032436/1‑150 (MQ=25)
                                                                                                                                           aGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAg            >  2:7076081/1‑150 (MQ=25)
                                                                                                                                           aGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCCTCGCCTCCAg            >  1:468260/1‑150 (MQ=255)
                                                                                                                                           aGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAg            >  1:4975669/1‑150 (MQ=255)
                                                                                                                                            ggtaTCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgc           >  1:4842747/4‑150 (MQ=25)
                                                                                                                                            ggAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgc           <  2:8737679/150‑1 (MQ=255)
                                                                                                                                             gtaTCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgcg          >  1:3156398/3‑150 (MQ=25)
                                                                                                                                             gAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCCATGTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATCCGGTCGATCTTCGCCGGCATCGCCTCCcgcg          >  2:4637092/1‑150 (MQ=255)
                                                                                                                                              taTCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgcgc         >  2:6339778/2‑150 (MQ=25)
                                                                                                                                              aaTCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgcgc         >  1:6905245/1‑150 (MQ=255)
                                                                                                                                                 cACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa      <  2:3349362/150‑1 (MQ=255)
                                                                                                                                                   ccAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACa    >  1:525670/1‑150 (MQ=255)
                                                                                                                                                   ccAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACa    >  2:7272286/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: