Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,870,871 T→A 5.8% E145V (GAA→GTA) ‡ 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,8710TA 5.8% 360.6 / 9.1 105E145V (GAA→GTA) ‡truDtRNA(Glu) pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base T (51/47);  new base A (4/2);  total (55/50)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.34e-01

Only 100 of 113 total aligned reads displayed.

CAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTT  >  NC_000913/2870724‑2871019
                                                                                                                                                   |                                                                                                                                                    
cAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                     >  3:151075/1‑149 (MQ=255)
cAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                     >  2:385310/1‑149 (MQ=255)
cAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                                                                                                                                     <  1:819767/149‑1 (MQ=255)
   aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                  <  2:746429/149‑1 (MQ=255)
   aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                  <  1:238545/149‑1 (MQ=255)
   aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                  <  3:208748/149‑1 (MQ=255)
   aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                  <  2:365713/149‑1 (MQ=255)
   aaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAAc                                                                                                                                                  <  1:855597/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)
             tATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGa                                                                                                                                        <  2:78426/149‑1 (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                                                                                                                <  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:459208/149‑1 (MQ=255)
                                               ccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGt                                                                                                      >  1:515527/1‑149 (MQ=255)
                                               ccAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGt                                                                                                      >  1:745797/1‑149 (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)
                                                                             gCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCTGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTgc                                                                        >  3:534700/1‑149 (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                                                                >  2:465378/1‑149 (MQ=255)
                                                                                     cAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGc                                                                >  1:245316/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)
                                                                                      aaaaCGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGcc                                                               >  2:652413/1‑149 (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:105584/1‑149 (MQ=255)
                                                                                                  aCCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCa                                                   <  1:29202/149‑1 (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                                 <  1:524438/149‑1 (MQ=255)
                                                                                                                      cTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTc                                   >  4:326652/1‑145 (MQ=255)
                                                                                                                      cTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTc                                   <  3:326652/145‑1 (MQ=255)
                                                                                                                       ttttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGtt                              >  4:406554/1‑149 (MQ=255)
                                                                                                                       ttttACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGtt                              >  1:659456/1‑149 (MQ=255)
                                                                                                                           aCGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGaaa                          <  4:170559/149‑1 (MQ=255)
                                                                                                                            cGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAg                         <  3:454340/149‑1 (MQ=255)
                                                                                                                            cGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAg                         <  1:313730/149‑1 (MQ=255)
                                                                                                                             gCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAgg                        >  3:614525/1‑149 (MQ=255)
                                                                                                                             gCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAgg                        >  3:309142/1‑149 (MQ=255)
                                                                                                                                   atCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCa                  >  2:213313/1‑149 (MQ=255)
                                                                                                                                     cgcgatgACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGa                >  4:655791/7‑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:220797/3‑149 (MQ=255)
                                                                                                                                           gACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGaaa                          <  2:186117/133‑1 (MQ=255)
                                                                                                                                           gACGTTGAACAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGaaa                          >  1:186117/1‑133 (MQ=255)
                                                                                                                                            aCGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCa         >  4:531779/1‑149 (MQ=255)
                                                                                                                                                 gttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGTGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTcc    <  2:551119/149‑1 (MQ=255)
                                                                                                                                                 gttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTcc    <  2:245316/149‑1 (MQ=255)
                                                                                                                                                 gttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTcc    <  3:854431/149‑1 (MQ=255)
                                                                                                                                                 gttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTcc    <  1:693760/149‑1 (MQ=255)
                                                                                                                                                  ttCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCt   >  1:424859/1‑149 (MQ=255)
                                                                                                                                                   tCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCtt  <  4:833617/149‑1 (MQ=255)
                                                                                                                                                   tCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCtt  >  4:203619/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTT  >  NC_000913/2870724‑2871019

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