Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,306,272 T→C 6.9% S110G (AGC→GGC)  PP_2028 ← TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,306,2720TC 6.9% 163.2 / 2.7 59S110G (AGC→GGC) PP_2028TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (33/21);  new base C (2/2);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

CAAATTGCGGGTGGTGGCCAGGTCGTCGGACAGCGGCACCAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGCCCAGCGCCACCCACGGCA  >  NC_002947/2306123‑2306421
                                                                                                                                                     |                                                                                                                                                     
cAAATTGCGGGTGGTGGCCAGGTCGTCGGACAGCGGCACCAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCt                                                                                                                                                       <  1:537437/150‑1 (MQ=255)
 aaaTTGCGGGTGGTGGCCAGGTCGTCGGACAGCGGCACCAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTg                                                                                                                                                      >  2:147276/1‑150 (MQ=255)
   aTTGCGGGTGGTGGCCAGGTCGTCGGACAGCGGCACCAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTc                                                                                                                                                    >  2:1012012/1‑150 (MQ=255)
        gggTGGTGGCCAGGTCGTCGGACAGCGGCACCAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCCCgg                                                                                                                                               <  2:705896/150‑1 (MQ=255)
                  cAGGTCGTCGGACAGCGGCACCAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTGGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCTATGCAGTGGGGCGGGCTGTCCTCGGCGAGCATTTg                                                                                                                                     >  2:641224/1‑150 (MQ=255)
                      tcgtcgGACAGCGGCACCAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCggg                                                                                                                                 <  2:537215/150‑1 (MQ=255)
                         tcgGACAGCGGCACCAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTc                                                                                                                              >  2:625133/1‑150 (MQ=255)
                                     aCCAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAgg                                                                                                                  >  2:549246/1‑150 (MQ=255)
                                     aCCAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAgg                                                                                                                  >  2:117544/1‑150 (MQ=255)
                                       cAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCgg                                                                                                                                  >  1:566558/1‑132 (MQ=255)
                                       cAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCgg                                                                                                                                  <  2:566558/132‑1 (MQ=255)
                                          ggTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCTAGCATTTGCGGGGTCAGTTCGAGCAGGATCAc                                                                                                             <  1:641224/150‑1 (MQ=255)
                                            tgtgCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCACCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCa                                                                                                           >  1:741686/1‑150 (MQ=255)
                                                   ggAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTc                                                                                                    <  1:794292/150‑1 (MQ=255)
                                                    gAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCg                                                                                                   >  2:438566/1‑150 (MQ=255)
                                                     aGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCgg                                                                                                  >  1:475816/1‑150 (MQ=255)
                                                            aTACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGa                                                                                                                       >  2:774303/1‑122 (MQ=255)
                                                            aTACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGa                                                                                                                       <  1:774303/122‑1 (MQ=255)
                                                                accaccAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTGGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGa                                                                                       >  1:258314/1‑150 (MQ=255)
                                                                accaccAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGa                                                                                       >  1:867474/1‑150 (MQ=255)
                                                                accaccAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGAGATGTGTGCTGTTATGATCACCAGCGGGTCGGCCGGGCGCTGa                                                                                       >  2:533463/1‑150 (MQ=255)
                                                                   accaGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTCACgg                                                                                    >  2:233819/1‑150 (MQ=255)
                                                                       gcgcCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGa                                                                                                                                            >  1:43196/1‑90 (MQ=255)
                                                                       gcgcCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGa                                                                                                                                            <  2:43196/90‑1 (MQ=255)
                                                                         gcCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTttct                                                                              <  1:552518/150‑1 (MQ=255)
                                                                            gTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGACTTTCTTcc                                                                           >  2:7832/1‑150 (MQ=255)
                                                                            gTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCCGTCCTCGGCGAGCATTTg                                                                                                                                     >  1:453844/1‑92 (MQ=255)
                                                                            gTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCCGTCCTCGGCGAGCATTTg                                                                                                                                     >  1:453766/1‑92 (MQ=255)
                                                                            gTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCCGTCCTCGGCGAGCATTTg                                                                                                                                     <  2:453844/92‑1 (MQ=255)
                                                                            gTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCCGTCCTCGGCGAGCATTTg                                                                                                                                     <  2:453766/92‑1 (MQ=255)
                                                                                 gCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCg                                                                      >  1:939520/1‑150 (MQ=255)
                                                                                       gcggcgGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCa                                                                <  2:867474/150‑1 (MQ=255)
                                                                                           cggTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTg                                                            <  1:438566/150‑1 (MQ=255)
                                                                                            ggTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGCCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTgg                                                           <  2:564618/150‑1 (MQ=255)
                                                                                               gtTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGcc                                                        <  1:117544/150‑1 (MQ=255)
                                                                                                ttCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGccc                                                       >  2:980166/1‑150 (MQ=255)
                                                                                                     gCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCGCCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGagca                                                  >  1:606311/1‑150 (MQ=255)
                                                                                                           cGAGGATCTTGCGCCGTGCCTTTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCa                                            <  2:758110/150‑1 (MQ=255)
                                                                                                            gAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAg                                           <  1:549246/150‑1 (MQ=255)
                                                                                                             aGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAgg                                          >  2:1014588/1‑150 (MQ=255)
                                                                                                              ggATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGa                                         >  2:165139/1‑150 (MQ=255)
                                                                                                                   ttGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGGTGATCGTCGGGCAGACCCGCGCCTGACTGGAACAGTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCAcc                                    >  1:520285/1‑150 (MQ=255)
                                                                                                                    tGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCGCCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCa                                   <  2:606311/150‑1 (MQ=255)
                                                                                                                         cGGGACTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAAc                              >  2:68557/1‑150 (MQ=255)
                                                                                                                           ggCCTGTTCCAGGCGGTTGGGCGGGCGGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAACGCCTGGATCACCAGTAACCa                            >  2:564998/1‑150 (MQ=255)
                                                                                                                            gCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAg                           <  1:1014588/150‑1 (MQ=255)
                                                                                                                             ccTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAgg                          >  2:760555/1‑150 (MQ=255)
                                                                                                                                gTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGAt                                        >  1:941179/1‑133 (MQ=255)
                                                                                                                                gTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGAt                                        <  2:941179/133‑1 (MQ=255)
                                                                                                                                  tCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAgg                     <  1:399903/150‑1 (MQ=255)
                                                                                                                                    cAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGcc                   >  1:726561/1‑150 (MQ=255)
                                                                                                                                    cAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGACGCTTCCATCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGcc                   >  1:520682/1‑150 (MQ=255)
                                                                                                                                     aGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGccc                  >  2:1071307/1‑150 (MQ=255)
                                                                                                                                     aGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGccc                  >  1:97942/1‑150 (MQ=255)
                                                                                                                                     aGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGccc                  <  1:1012012/150‑1 (MQ=255)
                                                                                                                                      ggCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGCCCa                 <  2:741686/150‑1 (MQ=255)
                                                                                                                                        cGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGCCCAgc               <  1:165139/150‑1 (MQ=255)
                                                                                                                                               gcgggcTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGCCCAGCGccaccc        >  1:736820/1‑150 (MQ=255)
                                                                                                                                                 gggcTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGCCCAGCGccacccac      >  2:1067761/1‑150 (MQ=255)
                                                                                                                                                   gcTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGCCCAGCGCCACCCACgg    >  2:843776/1‑150 (MQ=255)
                                                                                                                                                     tGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCAGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGCCCAGCGCCACCCACGGCa  <  2:475816/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAATTGCGGGTGGTGGCCAGGTCGTCGGACAGCGGCACCAGGGTGTGCGCGGAGCCGGCATACACCACCAGCGCCGTCTGGCTGTCGCGGCGGTGTTCAAGCAGGTCGAGGATCTTGCGCCGGGCCTGTTCCAGGCGGTTGGGCGGGCTGTCCTCGGCGAGCATTTGCGGGGTCAGTTCGAGCAGGATCACCAGCGGGTCGGCCGGGCGCTGACGGCTTTCTTCCAGCCGCTGCCAGCTGGGCCCGAGCAGCGCCAGGATCACCAGTAACCAGGCAAGGCCCAGCGCCACCCACGGCA  >  NC_002947/2306123‑2306421

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