Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,870,870 T→A 5.7% E145D (GAA→GAT) ‡ truD ← tRNA(Glu) pseudouridine(13) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,870,8700TA 5.7% 369.6 / 9.5 105E145D (GAA→GAT) ‡truDtRNA(Glu) pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base T (51/48);  new base A (4/2);  total (55/50)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 114 total aligned reads displayed.

ACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCT  >  NC_000913/2870722‑2871018
                                                                                                                                                    |                                                                                                                                                    
aCCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACgttgt                                                                                                                                                      <  2:198366/149‑1 (MQ=255)
  cAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                    >  3:151075/1‑149 (MQ=255)
     aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                 <  2:746429/149‑1 (MQ=255)
     aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                 <  1:238545/149‑1 (MQ=255)
     aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                 <  1:855597/149‑1 (MQ=255)
     aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                 <  2:365713/149‑1 (MQ=255)
     aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                 <  3:208748/149‑1 (MQ=255)
     aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                 >  4:237127/1‑149 (MQ=255)
           cGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATc                                                                                                                                           >  4:689078/1‑149 (MQ=255)
            gTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGAACAACGTctgat                                                                                                                                          <  4:590997/149‑5 (MQ=255)
             tttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACgttgt                                                                                                                                                      <  1:418239/136‑1 (MQ=255)
             tttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACgttgt                                                                                                                                                      >  2:418239/1‑136 (MQ=255)
                  tGCGATGGCGCACCGGAGTATTGGTTTCCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTGGTTCAACGTCATCGCGATTg                                                                                                                                    >  4:305659/1‑149 (MQ=255)
                     gATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTc                                                                                                                                 >  3:504181/1‑149 (MQ=255)
                           gcACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcg                                                                                                                           >  1:835377/1‑149 (MQ=255)
                            cACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcgc                                                                                                                          <  1:269337/149‑1 (MQ=255)
                            cACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcgc                                                                                                                          <  2:315579/149‑1 (MQ=255)
                               cGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcgc                                                                                                                          <  1:760560/146‑1 (MQ=255)
                               aggAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTAGCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcgc                                                                                                                          >  2:760560/2‑146 (MQ=255)
                                ggAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAg                                                                                                                      >  2:375209/1‑149 (MQ=255)
                                      ttGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTCGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAg                                                                                                                >  2:693760/1‑149 (MQ=255)
                                       tGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAgg                                                                                                               <  1:554628/149‑1 (MQ=255)
                                       tGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAgg                                                                                                               <  3:277573/149‑1 (MQ=255)
                                         gTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGt                                                                                                             <  3:662885/149‑1 (MQ=255)
                                               gcCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAACCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGc                                                                                                       <  4:493456/149‑1 (MQ=255)
                                               gcCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGc                                                                                                       <  3:460618/149‑1 (MQ=255)
                                                 ccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCCAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGt                                                                                                     >  3:709544/1‑149 (MQ=255)
                                                 ccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGt                                                                                                     <  4:190663/149‑1 (MQ=255)
                                                 ccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGt                                                                                                     >  1:745797/1‑149 (MQ=255)
                                                 ccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGt                                                                                                     >  1:515527/1‑149 (MQ=255)
                                                 ccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGt                                                                                                     <  1:459208/149‑1 (MQ=255)
                                                    gcgccgcgccCCCTGCAAGTTGCTACCGCCAATCCCAAACCGGTGGGAACCGAAGTAGGTCGGTACACCTTTTACGCAAATATCTATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAGCCAGTGTAAAGGctgtac                                                                                                  <  4:279411/149‑6 (MQ=255)
                                                    gcgcTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTAc                                                                                                  >  3:837067/1‑149 (MQ=255)
                                                    gcgcTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTAc                                                                                                  <  2:650527/149‑1 (MQ=255)
                                                     cgcTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTAcc                                                                                                 <  3:648167/149‑1 (MQ=255)
                                                            cccccTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAgc                                                                                          <  3:689078/149‑1 (MQ=255)
                                                              cccTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAgcgc                                                                                        >  1:461474/1‑149 (MQ=255)
                                                                cTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAgcgcgc                                                                                      >  4:588232/1‑149 (MQ=255)
                                                                  gCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCt                                                                                    <  2:515527/149‑1 (MQ=255)
                                                                      gTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAAc                                                                                >  1:719884/1‑149 (MQ=255)
                                                                      gTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAAc                                                                                <  3:237127/149‑1 (MQ=255)
                                                                      gTTGCTAACGCCGATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAGCGTCATCGCGATTACTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAAc                                                                                >  4:746752/1‑149 (MQ=255)
                                                                         gCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCa                                                                             <  4:116375/149‑1 (MQ=255)
                                                                               gCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTgc                                                                       >  4:800374/1‑149 (MQ=255)
                                                                                 cAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTgcgc                                                                     <  1:130021/149‑1 (MQ=255)
                                                                                  aaTCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCt                                                                    >  4:137115/1‑149 (MQ=255)
                                                                                   aTCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCtt                                                                   <  2:854547/149‑1 (MQ=255)
                                                                                   aTCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCtt                                                                   <  4:686289/149‑1 (MQ=255)
                                                                                    tCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTtg                                                                  <  2:719884/149‑1 (MQ=255)
                                                                                     cccAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTtgt                                                                 <  4:709544/149‑1 (MQ=255)
                                                                                      ccAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTtgtg                                                                >  4:125692/1‑149 (MQ=255)
                                                                                       cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGc                                                               >  1:245316/1‑149 (MQ=255)
                                                                                       cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGc                                                               >  2:465378/1‑149 (MQ=255)
                                                                                       cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGc                                                               >  3:336308/1‑149 (MQ=255)
                                                                                        aaaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCt                                                                                                >  4:434962/1‑115 (MQ=255)
                                                                                        aaaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCt                                                                                                <  3:434962/115‑1 (MQ=255)
                                                                                        aaaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGcc                                                              >  2:762521/1‑149 (MQ=255)
                                                                                          aaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCcgc                                                            <  3:305659/149‑1 (MQ=255)
                                                                                               tGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATa                                                       <  4:504181/149‑1 (MQ=255)
                                                                                                gggCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATAc                                                      <  1:385310/149‑1 (MQ=255)
                                                                                                gggCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATAc                                                      >  3:778787/1‑149 (MQ=255)
                                                                                                 ggCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACt                                                     >  2:277031/1‑149 (MQ=255)
                                                                                                   cACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTcc                                                   <  3:137115/149‑1 (MQ=255)
                                                                                                   cACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTCACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTcc                                                   >  3:507626/1‑149 (MQ=255)
                                                                                                    aCCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCa                                                  <  1:29202/149‑1 (MQ=255)
                                                                                                    aCCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCa                                                  >  1:105584/1‑149 (MQ=255)
                                                                                                      cGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTtgtg                                                                >  4:841276/1‑133 (MQ=255)
                                                                                                      cGAAGTAGCTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGc                                                >  4:256631/1‑149 (MQ=255)
                                                                                                       gAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCa                                               >  3:863052/1‑149 (MQ=255)
                                                                                                         agtagtTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCAcc                                             <  2:105584/149‑1 (MQ=255)
                                                                                                          gtagtTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCt                                            <  2:228386/149‑1 (MQ=255)
                                                                                                           tagtTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTg                                           >  2:159603/1‑149 (MQ=255)
                                                                                                                  gTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCt                                    >  1:368666/1‑149 (MQ=255)
                                                                                                                   tACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCtt                                   >  1:551119/1‑149 (MQ=255)
                                                                                                                      acCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCa                                >  2:313730/1‑149 (MQ=255)
                                                                                                                      acCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCa                                <  1:886037/149‑1 (MQ=255)
                                                                                                                      aatttttACGCAAATATCGATCAGACGTTGTGCGACGTCATCGCGATTGCTCACTTCGCCCAGGACCAGTGTAAAGGCCTAACCTTTCAGAGCGCCTAAATGCAGCTTGCGCTTGTGCAGTGAATACTCCAGCACATGGCAGCCTTCCa                                <  2:782909/146‑1 (MQ=255)
                                                                                                                        cTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTc                                  >  4:326652/1‑145 (MQ=255)
                                                                                                                         ttttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGtt                             >  4:406554/1‑149 (MQ=255)
                                                                                                                         ttttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGtt                             >  1:659456/1‑149 (MQ=255)
                                                                                                                          tttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCACTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTg                            <  2:648105/149‑1 (MQ=255)
                                                                                                                             aCGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGaaa                         <  4:170559/149‑1 (MQ=255)
                                                                                                                              cGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAg                        <  1:313730/149‑1 (MQ=255)
                                                                                                                              cGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAg                        <  3:454340/149‑1 (MQ=255)
                                                                                                                               gCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAgg                       >  3:614525/1‑149 (MQ=255)
                                                                                                                               gCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAgg                       >  3:309142/1‑149 (MQ=255)
                                                                                                                                         aTCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCTTTCCAGTTGAAAGGCGCTCAGATc             <  2:689787/149‑1 (MQ=255)
                                                                                                                                         aTCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATc             >  4:503724/1‑149 (MQ=255)
                                                                                                                                          tCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCg            >  2:434176/1‑149 (MQ=255)
                                                                                                                                           atgACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCgg           >  4:656398/3‑149 (MQ=255)
                                                                                                                                           atgACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCgg           >  4:220797/3‑149 (MQ=255)
                                                                                                                                             gACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGaaa                         >  1:186117/1‑133 (MQ=255)
                                                                                                                                             gACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGaaa                         <  2:186117/133‑1 (MQ=255)
                                                                                                                                              aCGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCa        >  4:531779/1‑149 (MQ=255)
                                                                                                                                                 ttgttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCAttt     >  3:525814/1‑149 (MQ=255)
                                                                                                                                                   gttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGTGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTcc   <  2:551119/149‑1 (MQ=255)
                                                                                                                                                   gttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTcc   <  3:854431/149‑1 (MQ=255)
                                                                                                                                                   gttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTcc   <  1:693760/149‑1 (MQ=255)
                                                                                                                                                   gttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTcc   <  2:245316/149‑1 (MQ=255)
                                                                                                                                                    ttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCt  >  1:424859/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCT  >  NC_000913/2870722‑2871018

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: