Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 961,037 A→G 6.2% I115T (ATT→ACT)  selU ← tRNA 2‑selenouridine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947961,0370AG 6.2% 184.3 / 2.2 64I115T (ATT→ACT) selUtRNA 2‑selenouridine synthase
Reads supporting (aligned to +/- strand):  ref base A (34/26);  new base G (2/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

TCCAGCTGGTGCAGCACATCGGTCTTGCCGGTACCGGTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGGCTTGCTTGGTGGCGCCACTGACCAGTTG  >  NC_002947/960893‑961182
                                                                                                                                                |                                                                                                                                                 
gccAGCTCCTGCAGCACATCGGCCTTGCCGGTACCGGTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATAcc                                                                                                                                               <  8:210650/148‑1 (MQ=255)
  cAGCTGGTGCAGCACATCGGTCTTGCCGGTACCGGTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATAccc                                                                                                                                              <  4:310160/148‑1 (MQ=255)
  cAGCTGGTGCAGCACATCGGTCTTGCCGGTACCGGTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAAAACCCg                                                                                                                                             >  3:269652/1‑149 (MQ=255)
   aGCTGGTGCAGCACATCGGTCTTGCCGGTACCGGTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGc                                                                                                                                            <  4:129515/149‑1 (MQ=255)
                            cGGTACCGGTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAAGACCCGCCTCATCGAGCAGCCAGCCCTGCAc                                                                                                                    >  6:18293/1‑148 (MQ=255)
                              gTACCGGTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGa                                                                                                                  >  6:209290/1‑148 (MQ=255)
                               tACTGGTCGGTGGACTGACCGGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGCATACACGCCTCATCGCGCAGCCAGCCCTGCACGATc                                                                                                                >  6:37453/1‑149 (MQ=255)
                                 ccGGTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTg                                                                                                              <  6:132719/149‑1 (MQ=255)
                                 ccGGTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTg                                                                                                              >  3:10943/1‑149 (MQ=255)
                                   ggTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTATTCAGTACCCGcgt                                                                                                                                          >  4:167132/1‑117 (MQ=255)
                                   ggTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTATTCAGTACCCGcgt                                                                                                                                          <  3:167132/119‑3 (MQ=255)
                                        gTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGc                                                                                                       <  6:87378/149‑1 (MQ=255)
                                           cACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAgg                                                                                                    >  7:292896/1‑149 (MQ=255)
                                           cACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAgg                                                                                                    >  7:288759/1‑149 (MQ=255)
                                           cACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAgg                                                                                                    >  6:158984/1‑149 (MQ=255)
                                           cACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAgg                                                                                                    >  5:276402/1‑149 (MQ=255)
                                           cACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAg                                                                                                     <  7:282667/148‑1 (MQ=255)
                                           cACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCATTGTACCCCCCTTTGACACGCGGGTAGTGAATACCCGCCTCATCGCCCAGCCACCCCTGAACGAACTGCGATCGCAgg                                                                                                    >  2:56909/1‑149 (MQ=255)
                                                         aaaTCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAg                                                                                        >  7:113125/1‑147 (MQ=255)
                                                          aaTCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGTCCGCCACGGAAGCAg                                                                                     >  8:285060/1‑149 (MQ=255)
                                                              gCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTAc                                                                                  <  4:10943/148‑1 (MQ=255)
                                                                aTTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAgg                                                                               >  6:38383/1‑149 (MQ=255)
                                                                    cTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAAc                                                                           >  7:202752/1‑149 (MQ=255)
                                                                    cTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAAc                                                                           >  6:230937/1‑149 (MQ=255)
                                                                     tCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAAcc                                                                          >  4:170376/1‑149 (MQ=255)
                                                                     tCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGATACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAAcc                                                                          >  7:299381/1‑149 (MQ=255)
                                                                      ccACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAAcag                                                                         >  6:177699/1‑147 (MQ=255)
                                                                        aCCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTc                                                                       <  8:113125/149‑1 (MQ=255)
                                                                         ccGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCACAGGCCGCCACGGAAGCAGTACAGGTAACCGt                                                                        >  5:24388/1‑147 (MQ=255)
                                                                            cTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAggg                                                                   >  1:261570/1‑149 (MQ=255)
                                                                                cTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTgc                                                                <  8:288759/148‑1 (MQ=255)
                                                                                cTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTgc                                                                <  8:292896/148‑1 (MQ=255)
                                                                                   ggTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTgc                                                                >  8:305816/1‑145 (MQ=255)
                                                                                   ggTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTgc                                                                <  7:305816/145‑1 (MQ=255)
                                                                                     gactctcCAGCATGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTg                                                          <  1:56909/146‑1 (MQ=255)
                                                                                          tcCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCaaa                                                     >  5:313233/1‑149 (MQ=255)
                                                                                          tcCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCaaa                                                     >  3:176254/1‑149 (MQ=255)
                                                                                          tcCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCaaa                                                     >  3:296035/1‑149 (MQ=255)
                                                                                          tcCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCaaa                                                     >  6:248718/1‑149 (MQ=255)
                                                                                          tcCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTAGGCaaa                                                     >  6:30745/1‑149 (MQ=255)
                                                                                              gcagGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCt                                                 <  5:248718/149‑1 (MQ=255)
                                                                                                 ggAAGGTACGCATGGCCTTGTAGCCGCCTTTGACAGGCGGGTACTGAATACCCGCCTCATCGCTCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGTAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGcc                                              <  4:269652/149‑1 (MQ=255)
                                                                                                  gAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCCGATACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGccc                                             <  8:299381/149‑1 (MQ=255)
                                                                                                    aGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAg                                           <  5:37453/149‑1 (MQ=255)
                                                                                                      gTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGc                                         <  5:209290/149‑1 (MQ=255)
                                                                                                          gCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTc                                      >  8:286338/1‑148 (MQ=255)
                                                                                                               gCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGTGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGAc                                 <  5:172975/148‑1 (MQ=255)
                                                                                                                ccTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCTCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACgg                               <  6:276402/149‑1 (MQ=255)
                                                                                                                ccTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACgg                               <  7:285060/149‑1 (MQ=255)
                                                                                                                  ttGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGgc                             >  8:122245/1‑149 (MQ=255)
                                                                                                                    gTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCa                                            <  1:322879/132‑1 (MQ=255)
                                                                                                                    gTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCa                                            >  2:322879/1‑132 (MQ=255)
                                                                                                                    gTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGgc                             <  5:158984/147‑1 (MQ=255)
                                                                                                                     tAGGCGCCTTTGACACGCGGGTACTGAATACCCGCCTAATCGCGCAGCCAGAGCTGCACGATCTGCGTTCGCAGGCCGCCTGGGATGGAGTACAGGTTACCGTCAGGGTGCGCGTTGGCATACGCTGAGCAGGCCTCGAGACggccagg                          >  1:283219/1‑144 (MQ=255)
                                                                                                                     tAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGgctt                           <  5:38383/148‑1 (MQ=255)
                                                                                                                         cgccTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGGcttgcttg                      <  5:230937/149‑1 (MQ=255)
                                                                                                                                gACACGCTGGTACTGAATACCCGCCTCATCTAGCAGCCAGACCTGCACGATCTGCGATCGCAGGACGCCACGGAAGCAGTACAGGTAACCGTCAGGCTGCGCCTTGGCAAACGCTGACCAGGCCTCCAGCCGGGCTTGCTTGGTTgcgc               >  6:242831/1‑149 (MQ=255)
                                                                                                                                 acacGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGGCTTGCTTGGTGgcgc               >  1:119709/1‑148 (MQ=255)
                                                                                                                                   aggcgGGTATTCAGTACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGGCTTGCTTGGTGGCGCCa             <  1:78764/146‑1 (MQ=255)
                                                                                                                                   aggcgGGTATTCAGTACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGGCTTGCTTGGTGGCGCCa             >  2:78764/3‑148 (MQ=255)
                                                                                                                                        ggTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGGCTTGCTTGGTGGCGCCACTGa         <  6:304450/147‑1 (MQ=255)
                                                                                                                                        ggTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGGCTTGCTTGGTGGCGCCACTGa         >  5:304450/1‑147 (MQ=255)
                                                                                                                                            cTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGGCTTGCTTGGTGGCGCCACTGACCAGtt   <  2:119709/149‑1 (MQ=255)
                                                                                                                                             tGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGGCTTGCTTGGTGGCGCCACTGACCAGTtg  >  4:208975/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TCCAGCTGGTGCAGCACATCGGTCTTGCCGGTACCGGTCAGTCCACCGACCAGTACGAAATCGCATTGCTCCACCGCTTGCTGGGTAGTCTCCAGCAGGAAGGTACGCATGGCCTTGTAGCCGCCTTTGACACGCGGGTACTGAATACCCGCCTCATCGCGCAGCCAGCCCTGCACGATCTGCGATCGCAGGCCGCCACGGAAGCAGTACAGGTAACCGTCAGGGTGCGCCTTGGCAAACGCTGCCCAGGCCTCCAGACGGGCTTGCTTGGTGGCGCCACTGACCAGTTG  >  NC_002947/960893‑961182

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