Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,693,853 T→C 12.0% V34A (GTC→GCC)  PP_2362 → usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,693,8530TC12.0% 187.1 / 14.6 75V34A (GTC→GCC) PP_2362usher protein
Reads supporting (aligned to +/- strand):  ref base T (31/35);  new base C (4/5);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

TCAACACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGCGGGCCGGGCGGTTGT  >  NC_002947/2693704‑2694002
                                                                                                                                                     |                                                                                                                                                     
tCAACACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGTCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGt                                                                                                                                                       >  2:2230920/1‑150 (MQ=255)
tCAACACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGt                                                                                                                                                       >  1:657880/1‑150 (MQ=255)
  aaCACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCa                                                                                                                                                     <  2:890484/150‑1 (MQ=255)
    cacTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGc                                                                                                                                                   <  1:1219163/150‑1 (MQ=255)
     acTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGcc                                                                                                                                                  >  2:60790/1‑150 (MQ=255)
       tCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGcccc                                                                                                                                                >  1:1417411/1‑150 (MQ=255)
       tCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGcccc                                                                                                                                                >  2:1129196/1‑150 (MQ=255)
         ggcggACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGACCATGGCCAGccccag                                                                                                                                              <  2:1491974/150‑3 (MQ=255)
             gACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTg                                                                                                                                          <  2:817210/150‑1 (MQ=255)
               cGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGt                                                                                                                                        >  1:331395/1‑150 (MQ=255)
                   aTCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTg                                                                                                                                    >  1:708108/1‑150 (MQ=255)
                   aTCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTg                                                                                                                                    >  2:1279719/1‑150 (MQ=255)
                       gCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCCTGTGTGGAATGGGTCCCGGCGACGACGTGCCGCTGGTGCCCGGGGCTGGCCGTAGTCGGGCCCGCGCGGGTACTGGCcg                                                                                                                                >  2:113847/1‑150 (MQ=255)
                         cAAGGCCAATGAGCCAACGGTCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCcgac                                                                                                                              <  1:2230920/150‑1 (MQ=255)
                           aGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCcgacga                                                                                                                            <  1:1780494/150‑1 (MQ=255)
                               cAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTg                                                                                                                        >  1:559361/1‑150 (MQ=255)
                                      ccAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGGCCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCAc                                                                                                                 >  1:583847/1‑150 (MQ=255)
                                      ccAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCAc                                                                                                                 >  2:1933733/1‑150 (MQ=255)
                                       cAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCAcc                                                                                                                >  1:1135638/1‑150 (MQ=255)
                                               gAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACgg                                                                                                        >  2:160237/1‑150 (MQ=255)
                                               gAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACgg                                                                                                        >  1:1130262/1‑150 (MQ=255)
                                                aaTGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGga                                                                                                       <  2:1130262/150‑1 (MQ=255)
                                                aaTGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGga                                                                                                       <  2:708108/150‑1 (MQ=255)
                                                  tGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGgaga                                                                                                     >  1:844124/1‑150 (MQ=255)
                                                       ggCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTcc                                                                                                <  2:1620646/149‑1 (MQ=255)
                                                         cAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGc                                                                                              <  1:60790/150‑1 (MQ=255)
                                                             ggCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATc                                                                                          >  1:2225549/1‑150 (MQ=255)
                                                               ccAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGACCATGGCCAGccccag                                                                                                                                              >  1:52235/1‑94 (MQ=255)
                                                               ccAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGACCATGGCCAGccccag                                                                                                                                              <  2:52235/96‑3 (MQ=255)
                                                                       gaCCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGcc                                                                                                                      >  2:1288149/1‑112 (MQ=255)
                                                                       gaCCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGcc                                                                                                                      <  1:1288149/112‑1 (MQ=255)
                                                                       gaCCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCa                                                                                >  1:956921/1‑150 (MQ=255)
                                                                         cccGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCAcc                                                                              <  1:1847211/150‑1 (MQ=255)
                                                                          ccGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCAccc                                                                             >  1:1032766/1‑150 (MQ=255)
                                                                          ccGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCAccc                                                                             <  2:891145/150‑1 (MQ=255)
                                                                                aTGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGATATGTCCGCCATCGCCGAGGCCACCCTTTAcc                                                                       <  2:527462/150‑1 (MQ=255)
                                                                                 tGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCt                                                                      <  1:1529921/150‑1 (MQ=255)
                                                                                  ggTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTc                                                                     <  2:1417411/150‑1 (MQ=255)
                                                                                    ttGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGa                                                                   <  2:2010221/150‑1 (MQ=255)
                                                                                      ggAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGAcc                                                                 >  1:637126/1‑150 (MQ=255)
                                                                                        aaaCCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCtg                                                               <  1:541726/150‑1 (MQ=255)
                                                                                         aaCCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCtgg                                                              <  1:380234/150‑1 (MQ=255)
                                                                                            cGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCtggtgg                                                           <  2:956921/150‑1 (MQ=255)
                                                                                                tgtgAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGaa                                                       >  1:1837318/1‑150 (MQ=255)
                                                                                                tgtgAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGaa                                                       >  2:2139291/1‑150 (MQ=255)
                                                                                                  tgAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAAcc                                                     <  2:657880/150‑1 (MQ=255)
                                                                                                    aaaTGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGAc                                                                  <  2:1247949/135‑1 (MQ=255)
                                                                                                    aaaTGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGAc                                                                  >  1:1247949/1‑135 (MQ=255)
                                                                                                    aaaTGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAg                                                   >  2:332558/1‑150 (MQ=255)
                                                                                                             ccGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAgg                                          <  2:637126/150‑1 (MQ=255)
                                                                                                             ccGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAgg                                          <  1:1618303/150‑1 (MQ=255)
                                                                                                              cGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAggg                                         <  2:559361/150‑1 (MQ=255)
                                                                                                              cGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAggg                                         <  1:1187106/150‑1 (MQ=255)
                                                                                                                   acgacgTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGAc                                    <  2:1837318/150‑1 (MQ=255)
                                                                                                                   acgacgTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGAc                                    <  2:844124/150‑1 (MQ=255)
                                                                                                                        gTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGAt                               >  1:1449857/1‑150 (MQ=255)
                                                                                                                        gTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGAt                               >  1:965026/1‑150 (MQ=255)
                                                                                                                              cTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGt                         <  1:1279719/150‑1 (MQ=255)
                                                                                                                               tGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTg                        <  2:1449857/150‑1 (MQ=255)
                                                                                                                                 gTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTgca                      <  1:1129196/150‑1 (MQ=255)
                                                                                                                                 cagcgcggggCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGGCCAGGGTTGACCTGATACTGATgca                      >  2:218401/7‑150 (MQ=255)
                                                                                                                                       ggggCTGACCATGGCCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGcg                >  2:389682/1‑150 (MQ=255)
                                                                                                                                        gggCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGcgg               <  1:218401/150‑1 (MQ=255)
                                                                                                                                        gggCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGcgg               <  2:1135638/150‑1 (MQ=255)
                                                                                                                                        gggCTGACCATGGCCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAgg                                          <  1:986518/123‑1 (MQ=255)
                                                                                                                                        gggCTGACCATGGCCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAgg                                          >  2:986518/1‑123 (MQ=255)
                                                                                                                                         ggCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGcggg              >  2:1985752/1‑150 (MQ=255)
                                                                                                                                         ggCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGcggg              <  2:331395/150‑1 (MQ=255)
                                                                                                                                            tGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGAACACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGcgggccg           >  2:1796439/1‑150 (MQ=255)
                                                                                                                                             gaccATGGCCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGcgggccgg          <  1:389682/148‑1 (MQ=255)
                                                                                                                                               ccATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGcgggccgggc        >  1:954964/1‑150 (MQ=255)
                                                                                                                                                cATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATAc                             >  2:2255964/1‑128 (MQ=255)
                                                                                                                                                cATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCGTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATAc                             <  1:2255964/128‑1 (MQ=255)
                                                                                                                                                cATGGCCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATAc                             >  2:1630587/1‑128 (MQ=255)
                                                                                                                                                cATGGCCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATAc                             <  1:1630587/128‑1 (MQ=255)
                                                                                                                                                   ggTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGTAGCGGGCCGGGCGGtt    <  2:965026/150‑1 (MQ=255)
                                                                                                                                                    gTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGCGGGCCGGCCGGTTg   >  2:1065683/1‑150 (MQ=255)
                                                                                                                                                     tCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGCGGGCCGGGCGGTTGt  >  1:599450/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAACACTCGGCGGACGCCATCCGCCAAGGCCAATGAGCCAACGGCCGAATGAAGTGCAGGGGCCAGGGAGGACCCGGCAATGGTTGGAAACCGTGTGTGAAATGGCTCCCGGCGACGACGTGCCGCTGGTGCCTGGGGCTGGCCATGGTCAGCCCCGCGCTGGTACTGGCCGACGACTTGCCCCCACCGCCCACGGAGATGTCCGCCATCGCCGAGGCCACCCTTTACCTCGACCTGGTGGTGAACCAGATGCCCAGGGTTGACCTGATACCGGTGCAGCAGCGGGCCGGGCGGTTGT  >  NC_002947/2693704‑2694002

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