Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 159,811 C→T 5.6% W101* (TGG→TAG)  gluQ ← glutamyl‑Q tRNA(Asp) synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913159,8110CT 5.6% 236.1 / 2.8 71W101* (TGG→TAG) gluQglutamyl‑Q tRNA(Asp) synthetase
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base T (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

ACCGGATGCTGCTGGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGTCGATATCTTCTATGCG  >  NC_000913/159667‑159956
                                                                                                                                                |                                                                                                                                                 
aCCGGATGCTGCTGGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATCTAACCAgg                                                                                                                                                >  3:159036/1‑148 (MQ=255)
   ggATGCTGCTGGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGc                                                                                                                                               >  1:269897/1‑146 (MQ=255)
    gATGCTGCTGGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCg                                                                                                                                              <  2:269897/146‑1 (MQ=255)
      tgctgcTGGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTg                                                                                                                                          <  4:135142/148‑1 (MQ=255)
               cGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGa                                                                                                                                  >  3:183097/1‑147 (MQ=255)
                     cGCACTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATa                                                                                                                                >  7:107770/1‑143 (MQ=255)
                     cGCACTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATa                                                                                                                                <  8:107770/143‑1 (MQ=255)
                       cACTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGcgtcg                                                                                                                         <  3:161953/148‑1 (MQ=255)
                       cACTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGcgtcg                                                                                                                         >  4:257749/1‑148 (MQ=255)
                         cTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTg                                                                                                                       >  4:224935/1‑148 (MQ=255)
                           gcggcgTTGTCTGGTCCATGATGCAACACCCTGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCa                                                                                                                                                  <  4:133742/119‑1 (MQ=255)
                           gcggcgTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCa                                                                                                                                                  >  3:133742/1‑119 (MQ=255)
                            cggcgTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACg                                                                                                                    >  7:141019/1‑148 (MQ=255)
                              gcgTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGtt                                                                                                                  >  2:23457/1‑148 (MQ=255)
                              gcgTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGtt                                                                                                                  >  5:234217/1‑148 (MQ=255)
                                  tGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGa                                                                                                              >  3:225035/1‑148 (MQ=255)
                                    tCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGAtt                                                                                                            <  6:142507/148‑1 (MQ=255)
                                       ggTCCATGATGTAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGc                                                                                                          >  1:195843/1‑147 (MQ=255)
                                        gTCCATGATGTAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGcc                                                                                                         <  2:195843/147‑1 (MQ=255)
                                        gTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTg                                                                                                           <  4:248346/145‑1 (MQ=255)
                                        gTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGcc                                                                                                         <  8:60231/147‑1 (MQ=255)
                                            atgatgCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGcgt                                                                                                                           >  5:119966/1‑125 (MQ=255)
                                            atgatgCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGcgt                                                                                                                           <  6:119966/125‑1 (MQ=255)
                                            atgatgCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCagag                                                                                                     >  1:74658/1‑147 (MQ=255)
                                                        ccGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTcc                                                                                        <  8:141019/148‑1 (MQ=255)
                                                         cGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTccc                                                                                       <  5:314872/148‑1 (MQ=255)
                                                          ggCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTcc                                                                                        >  6:144988/1‑146 (MQ=255)
                                                          ggCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCa                                                                                      <  4:225035/148‑1 (MQ=255)
                                                                 aCCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGcc                                                                                                         <  7:32741/122‑1 (MQ=255)
                                                                 aCCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGcc                                                                                                         >  8:32741/1‑122 (MQ=255)
                                                                 aCCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGa                                                                               <  2:74658/148‑1 (MQ=255)
                                                                     tcgtAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGt                                                                           <  7:21511/148‑1 (MQ=255)
                                                                       gtAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTa                                                                          <  5:144988/147‑1 (MQ=255)
                                                                          aaTACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCa                                                                                      <  8:219722/132‑1 (MQ=255)
                                                                          aaTACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCa                                                                                      >  7:219722/1‑132 (MQ=255)
                                                                          aaTACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATg                                                                       <  1:23457/147‑1 (MQ=255)
                                                                                     tGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGgcg                                                           >  7:39427/1‑148 (MQ=255)
                                                                                      gCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTGCTTGTAACCAGGCGAGTGCTTCCCGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGgcgc                                                          >  1:296141/1‑148 (MQ=255)
                                                                                        tttGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGt                                                                           >  1:83727/1‑129 (MQ=255)
                                                                                        tttGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGt                                                                           <  2:83727/129‑1 (MQ=255)
                                                                                        tttGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAg                                                        >  6:128301/1‑148 (MQ=255)
                                                                                          tGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAgg                                                       <  8:39427/147‑1 (MQ=255)
                                                                                            aaTACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGAt                                                     <  4:183097/147‑1 (MQ=255)
                                                                                             atccgcgCACGCGTACAGGCTCAGTAATAACTTAGTCCTTGTTCAGGTAACCAGGCGAGTTCTTCACGATAGGCGGCGTGACGTGGCGATTGCCAGAGAACATCGCCGTCCCAGTGAAGACCGTAATGTTCCAGCTGGCGCAGTATAg                                                   <  4:159036/145‑1 (MQ=255)
                                                                                             aTACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAg                                                   >  3:205855/1‑148 (MQ=255)
                                                                                             aTACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAg                                                   <  5:247481/148‑1 (MQ=255)
                                                                                                cgcgCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGttt                                                <  5:128301/148‑1 (MQ=255)
                                                                                                          tACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCAcc                                      <  7:253753/148‑1 (MQ=255)
                                                                                                           aCAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTc                                                                    <  1:188013/117‑1 (MQ=255)
                                                                                                           aCAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTc                                                                    >  2:188013/1‑117 (MQ=255)
                                                                                                               gTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGaa                                 >  2:229102/1‑148 (MQ=255)
                                                                                                                      aataaCTTAGTCCTTGTTCATGTAACTAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGgcgc                                                          >  1:240957/1‑116 (MQ=255)
                                                                                                                      aataaCTTAGTCCTTGTTCATGTAACTAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGgcgc                                                          >  1:300375/1‑116 (MQ=255)
                                                                                                                      aataaCTTAGTCCTTGTTCATGTAACTAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGgcgc                                                          <  2:240957/116‑1 (MQ=255)
                                                                                                                      aataaCTTAGTCCTTGTTCATGTAACTAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGgcgc                                                          <  2:300375/116‑1 (MQ=255)
                                                                                                                      aataaCTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCAc                           >  5:127639/1‑147 (MQ=255)
                                                                                                                      aataaCTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATACGCGTCGTGACGTAGCGATTGCCAGAGAACATCGCCCACCCAGTGCAGACCCTAAAGTTCCAGCTGGCGCAGGATACTTTCTGCGGCACCGGGAACTTAAc                           >  5:36753/1‑147 (MQ=255)
                                                                                                                        taaCTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAg                        >  3:210281/1‑148 (MQ=255)
                                                                                                                         aaCTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAgg                       >  4:51401/1‑148 (MQ=255)
                                                                                                                           cTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCg                     <  7:144205/148‑1 (MQ=255)
                                                                                                                               gTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGt                  <  4:205855/147‑1 (MQ=255)
                                                                                                                                  cTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGTCGa               <  4:267857/147‑1 (MQ=255)
                                                                                                                                     gTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGTCGatat            <  4:210281/147‑1 (MQ=255)
                                                                                                                                     gTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGTCGATAtc           >  3:128871/1‑148 (MQ=255)
                                                                                                                                         aTGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGTCGATAtcttc        <  6:234217/147‑1 (MQ=255)
                                                                                                                                            tAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGTCGATAtctt         <  2:185054/143‑1 (MQ=255)
                                                                                                                                            tAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGTCGATAtctt         >  1:185054/1‑143 (MQ=255)
                                                                                                                                            tAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGTCGATATCTTCTATg    <  3:51401/148‑1 (MQ=255)
                                                                                                                                             aaCCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTc                             <  7:174250/122‑1 (MQ=255)
                                                                                                                                             aaCCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTc                             >  8:174250/1‑122 (MQ=255)
                                                                                                                                              aCCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGTCGATATCTTCTATGCg  >  3:105301/1‑148 (MQ=255)
                                                                                                                                               ccAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGTCGATATCTTCTATGCg  >  2:44108/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ACCGGATGCTGCTGGCGGATACGCACTGCGGCGTTGTCTGGTCCATGATGCAACACCCGGCAATGACCGTCGTAAATACCGCCAATGCTTTGAATACGCGCACGCGTACAGGTGCAGTAATAACTTAGTCCTTGTTCATGTAACCAGGCGAGTGCTTCACGATAGGCGTCGTGACGTTGCGATTGCCAGAGAACATCGCCGTCCCAGTGCAGACCGTAATGTTCCAGCTGGCGCAGGATAGTTTCTGCGGCACCGGGAACTTCACGAGGCGGGTCGATATCTTCTATGCG  >  NC_000913/159667‑159956

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