Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 694,941 G→A L5L (CTC→CTT miaB ← isopentenyl‑adenosine A37 tRNA methylthiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913694,9410GA96.2% 149.4 / ‑5.6 54L5L (CTC→CTTmiaBisopentenyl‑adenosine A37 tRNA methylthiolase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (26/24);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATT  >  NC_000913/694700‑695186
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
aTGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCCGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGCAGTTTTTTg                                                                                                                                                                                                                                               >  2:384437/1‑251 (MQ=255)
aTGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTg                                                                                                                                                                                                                                               <  1:16567/251‑1 (MQ=255)
    cACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCCGGTTTTAATATGAAGTTTTTTGGta                                                                                                                                                                                                                                           >  2:263526/1‑249 (MQ=255)
              ccGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGct                                                                                                                                                                                                                                >  1:124025/1‑251 (MQ=255)
                   gCCGAAAACCAGGTCGGGATTTTTCTCTTTTAACAGTTTCCAGCGACCCAACCGATGGAAGACTTTTTCCTTAGCCTTCTCGCGGATGGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGc                                                                                                                                                                                                                           <  2:301787/251‑1 (MQ=255)
                              ggTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGc                                                                                                                                                                                                                >  2:45584/1‑251 (MQ=255)
                               gTCTGGATTCGTCGCTTTTAACAGTCGCCAGCTACCCAACGGATGGCAGACTTTTTCCTGAGCCTGCTCGCGGATTAGGCAGGTGGTCAGCAGCAGCACACCCGCTTCTGCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGcc                                                                                                                                                                                                               <  2:7459/251‑1 (MQ=255)
                                   ggtttcttcTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCATTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGAGGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAg                                                                                                                                                                                                             <  2:394734/246‑1 (MQ=255)
                                      ttcttcTCTTTTAACAGTTTCCAGCGACCAAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGACGCACATCCGCTTATTCCGCAACGTCGGTCAGTTGTTAGCCGTGGGTGACCTCCAGCAGATCGGCCAACTTCGATGAATCGTACTCGTTCATTTGCCAGACACAGGTTTTATTATGAAGTTTTTTGATCATCGCATTGCTATTGCGCAAtagtatataggaagt                                                                                                                                                                                                        >  2:264917/1‑241 (MQ=255)
                                             ctTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGgcg                                                                                                                                                                                                 >  1:208960/1‑251 (MQ=255)
                                                ttAACAGTTTCCAGCGACCCAACTAATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCa                                                                                                                                                                                               <  1:274273/250‑1 (MQ=255)
                                                  atctgTTTCCAGCGCCCCAACTGATGGAAGACTTTTTCCTGGGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACAGCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATa                                                                                                                                                                                             <  2:449974/246‑1 (MQ=255)
                                                     aGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCCATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTTCGAAATAGTAGCCAGGAATGCAGGGCGCAGAgtg                                                                                                                                                                                          >  2:255202/1‑250 (MQ=255)
                                                      gTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATa                                                                                                                                                                                             >  1:449974/1‑246 (MQ=255)
                                                       tttCCAGCGACCCAACTGATGGGAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGCATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGCAGGTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCCGGGCGCCTAGTGTa                                                                                                                                                                                        >  2:97256/1‑250 (MQ=255)
                                                                 cccAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTtgctg                                                                                                                                                                             >  2:219901/1‑251 (MQ=255)
                                                                   cAACTGATGGAAGTCTTTGTCCTGAGCTTTCTCGCGGATTGAGCAGTTGTTCAGCAGCAGCACACCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGGGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGATCTTGCGAAATAGTAGCCAGGAAGGCAGGGCGCATAGTGTAATGCTTtgctgct                                                                                                                                                                           <  2:88208/251‑2 (MQ=255)
                                                                         aTGGAAGATTTTTTCCTAAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCACCAGGTTTTAATATGGCGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGCATGCAGGGCGCTTAGTGTAATGCTTTGCTGCCGTtgtg                                                                                                                                                                     >  1:75794/1‑251 (MQ=255)
                                                                           ggAAGACTTTGTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGAc                                                                                                                                                                   >  1:384275/1‑251 (MQ=255)
                                                                             aaGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGAcc                                                                                                                                                                  >  1:264192/1‑250 (MQ=255)
                                                                                               ttCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGtt                                                                                                                                                                        >  2:129146/1‑226 (MQ=255)
                                                                                               ttCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGtt                                                                                                                                                                        <  1:129146/226‑1 (MQ=255)
                                                                                                  cgacctaTTGGGCAGGTTTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCtt                                                                                                                                            <  2:271382/245‑1 (MQ=255)
                                                                                                    gcgaaTTGAGCAGGTGCTAGACACCAGCCCATCCGCTGCTTCCGCCCTGTCGGTCAGTTGATAGCCGTGGGTCGCAGCCAGCCGATCGGCCATCGTCGATGAATCGTACTCGTTCATCGGACAGCCCCAGGTTTTAATAGGAAGTGTTTTGGTCATCGACTTGCTCTTGCGAAAGAGAAGCCATCAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCTTTATCAGTCCTTAg                                                                                                                                          <  2:462892/247‑1 (MQ=255)
                                                                                                    gcgGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAg                                                                                                                                          <  2:264192/251‑1 (MQ=255)
                                                                                                        aTTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCtt                                                                                                                                            >  1:271382/1‑245 (MQ=255)
                                                                                                            aGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAAtagt                                                                                                                                                                                                                   <  1:130457/170‑1 (MQ=255)
                                                                                                            aGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGAACTCGTTCATCTGACAGCCCCAGGTTTTAATATGACGTTTTTTGGTCATCGACTTGCTCTTGCGAAAtagt                                                                                                                                                                                                                   >  2:130457/1‑170 (MQ=255)
                                                                                                                      cagcagcagCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTaaa                                                                                                                        <  1:327736/251‑1 (MQ=255)
                                                                                                                           agcagcaCATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTa                                                                                                                    <  1:3897/251‑1 (MQ=255)
                                                                                                                             cagctcacccGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTaaa                                                                                                                  <  1:374566/243‑1 (MQ=255)
                                                                                                                                  caTCCGCTTCTTCCGCCACGTCGGACAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCTTCGACTTGCTCTTGCGAAATAGTAGCCAGGCATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCat                                                                                                              >  2:76960/1‑250 (MQ=255)
                                                                                                                                             tcCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAgg                                                                                                  >  1:191227/1‑251 (MQ=255)
                                                                                                                                              ccGCCACGGCGGTCAGGTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTCCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGga                                                                                                 <  1:381031/251‑1 (MQ=255)
                                                                                                                                                     gTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAt                                                                                                                                                                                              <  1:267625/151‑1 (MQ=255)
                                                                                                                                                     gTCGGACAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGCAATAGTAGCCAGGAATGCAGGGCGCAt                                                                                                                                                                                              >  2:267625/1‑151 (MQ=255)
                                                                                                                                                        ggTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGAcc                                                                                        >  1:406642/1‑250 (MQ=255)
                                                                                                                                                                          ggtggAATCCAGCCGCTCGGCCTTCTTCGATGAATCGTACTCGTTCATCTGACAGCCACAGGTTTTCATATGAGGTTTTTTGGTCATCGACTTGCTCTTGCGAATTAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCCGCCCTTAGGGGTGAAAATCCTGTAAACTTACAGCATATTGCCAACAGGTTGTATGACCATGACAAATCCACCAACgg                                                                     >  2:203482/1‑251 (MQ=255)
                                                                                                                                                                             ggCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGaa                                                                   >  1:102143/1‑250 (MQ=255)
                                                                                                                                                                               cATCCAGCAGATCGTCCTGTGTCGATGAATTGTCCGCGTTCATCTGACAGCCCATGGTTTTAATATGGATTTTTTGTGTGATCGGCTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGAAATGCTTTTCTTCCGTTGTGACCAGTATGAGCGTTATAAGCCCTTAGGTTTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGAGGATTTACCATGACAACTCAACCATCGTAAAtt                                                                <  2:13643/251‑1 (MQ=255)
                                                                                                                                                                                aTCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTGCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTg                                                               <  1:89607/251‑1 (MQ=255)
                                                                                                                                                                                        gATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTGAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTc                                                       <  1:255202/251‑1 (MQ=255)
                                                                                                                                                                                                     tCGATTAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGGTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGaa                                             <  2:449320/249‑1 (MQ=255)
                                                                                                                                                                                                       gATGAATCTTTCTCGTTCATCTGACATCCCCTGGTTTTAATATGAAGTTTTTTGTTCATCGACTTGCTCTTGCGAAATAGTAGCCATGAATGCAGGGCGCATAGTTTAATGCTTTGCTTCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACCGGATGATTGACCATGACAAATCAACCAACGTAAATTGCCATTGTCGGCTGAGGAATgg                                          <  2:333884/250‑1 (MQ=255)
                                                                                                                                                                                                        aTGAATCGTAACAGTTAATCTGACAGACCCAGGTTTTAATAAGAAGTTTGTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATGGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTc                                        <  2:522/251‑1 (MQ=255)
                                                                                                                                                                                                                        cATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAgg                                                                                                                                         <  2:382871/138‑1 (MQ=255)
                                                                                                                                                                                                                        cATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAgg                                                                                                                                         >  1:382871/1‑138 (MQ=255)
                                                                                                                                                                                                                           cTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAACTAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTgg                     >  2:451528/1‑251 (MQ=255)
                                                                                                                                                                                                                             gACAGCCGCAGGATTTAATATGAAGTTTTTTGGTCATCGGCTTGCTCGTGCGAAATAGTAGCCAGGAATGCAGGGCGCATATTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAATCATATTGCTAACAGGATGAGTGACCATGACAAATCAACCAACGGAAAGTGCCATTGTCGGCGTATGAATGGTCGGCGTCGCTCTTGCGCTgt                     <  2:44888/249‑2 (MQ=255)
                                                                                                                                                                                                                              aCAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCt                 >  1:469364/1‑251 (MQ=255)
                                                                                                                                                                                                                              aCAGCCACAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGAGCTggagc                  >  2:193108/1‑248 (MQ=255)
                                                                                                                                                                                                                                    ccAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGcac            >  1:148280/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGcac            >  2:184826/1‑249 (MQ=255)
                                                                                                                                                                                                                                              atatGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGAtt  >  2:72323/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
ATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGATT  >  NC_000913/694700‑695186

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: