Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 672,245 C→A 6.6% V847L (GTG→TTG) ‡ leuS ← leucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913672,2450CA 6.6% 461.8 / 15.3 136V847L (GTG→TTG) ‡leuSleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (67/60);  new base A (4/5);  total (71/65)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

Only 100 of 136 total aligned reads displayed.

GAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCT  >  NC_000913/672099‑672393
                                                                                                                                                  |                                                                                                                                                    
gaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACttt                                                                                                                                                   <  2:794809/150‑1 (MQ=255)
 catcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa                                                                                                                                                  >  2:197421/4‑150 (MQ=255)
 aaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTa                                                                                                                                                  >  2:152989/1‑150 (MQ=255)
   ccTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACg                                                                                                                                                >  1:692504/1‑150 (MQ=255)
    cTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGt                                                                                                                                               <  1:403421/150‑1 (MQ=255)
    cTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACACCAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGt                                                                                                                                               >  2:299586/1‑150 (MQ=255)
      gCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTAc                                                                                                                                             >  1:1039032/1‑150 (MQ=255)
      gCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTAc                                                                                                                                             <  1:1133713/150‑1 (MQ=255)
      gCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAGTGATTTACGTAc                                                                                                                                             <  1:1345907/150‑1 (MQ=255)
         tcgTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAGTGATTTACGTAccag                                                                                                                                          <  1:102969/148‑4 (MQ=255)
         tGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAGTGATTTACGTAccag                                                                                                                                          <  2:342261/150‑4 (MQ=255)
           gTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAGTGATTTACGTAccaggt                                                                                                                                        <  1:563470/150‑6 (MQ=255)
             aTCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACg                                                                                                                                      >  2:515940/1‑150 (MQ=255)
               cACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt                                                                                                                                                    <  1:215193/134‑1 (MQ=255)
               cACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt                                                                                                                                                    >  2:215193/1‑134 (MQ=255)
                aCGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCa                                                                                                                                   >  2:113344/1‑150 (MQ=255)
                aCGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCa                                                                                                                                   <  1:41867/150‑1 (MQ=255)
                    aGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCaa                                                                                                                               >  1:1119260/1‑150 (MQ=255)
                       tGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGat                                                                                                                            >  2:1213174/1‑150 (MQ=255)
                             ccaCACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTg                                                                                                                      <  2:898088/150‑1 (MQ=255)
                              cacaCCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGc                                                                                                                     >  2:1457000/1‑150 (MQ=255)
                                 ctctggcggTGTTAATCACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTAc                                                                                                                  >  1:420081/5‑150 (MQ=255)
                                  cccGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTAcc                                                                                                                 >  1:1268617/1‑150 (MQ=255)
                                   ccGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCa                                                                                                                >  2:134334/1‑150 (MQ=255)
                                      gcggTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGAt                                                                                                             <  1:380222/150‑1 (MQ=255)
                                          tGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTc                                                                                                         >  1:386144/1‑150 (MQ=255)
                                           gATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTcc                                                                                                        <  2:1039032/150‑1 (MQ=255)
                                            aTTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCt                                                                                                       >  1:475891/1‑150 (MQ=255)
                                             ttAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTg                                                                                                      <  1:152989/150‑1 (MQ=255)
                                              tAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTgg                                                                                                     >  1:680872/1‑150 (MQ=255)
                                                acacCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGcc                                                                                                   <  2:895453/150‑1 (MQ=255)
                                                   ccGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGc                                                                                                <  1:1430005/150‑1 (MQ=255)
                                                     gCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCAc                                                                                              >  1:323571/1‑150 (MQ=255)
                                                          agaTAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTcg                                                                                         <  1:323370/150‑1 (MQ=255)
                                                           gaTAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTcgc                                                                                        <  2:51840/150‑1 (MQ=255)
                                                           gaTAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTcgc                                                                                        >  1:1244834/1‑150 (MQ=255)
                                                           gaTAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTcgc                                                                                        <  2:546217/150‑1 (MQ=255)
                                                             tAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGa                                                                                      >  1:1084758/1‑150 (MQ=255)
                                                                caacaaTGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAAcc                                                                                   <  2:1368357/150‑1 (MQ=255)
                                                                caacaaTGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAAcc                                                                                   >  1:591051/1‑150 (MQ=255)
                                                                caacaaTGTTGCCAGATATCGCACGCTTACTCCCGCGCTTAGACAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAAcc                                                                                   >  2:503683/1‑150 (MQ=255)
                                                                     aTGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGttc                                                                              <  2:386144/150‑1 (MQ=255)
                                                                      tGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGttct                                                                             <  2:549239/150‑1 (MQ=255)
                                                                             aGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTg                                                                      <  2:1309154/150‑1 (MQ=255)
                                                                             aGAAATCGCACGCTTCCTCCCGCGCTTAGCCAACGACAAGATTGAGGAGTTTACCTGCTACGTAAATCACTTTACGTACACTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTg                                                                      >  1:1361164/1‑150 (MQ=255)
                                                                              gATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGc                                                                     <  2:420081/150‑1 (MQ=255)
                                                                                  tCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTcc                                                                 <  1:515940/150‑1 (MQ=255)
                                                                                  tCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTcc                                                                 <  1:875311/150‑1 (MQ=255)
                                                                                   cGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCa                                                                >  1:1219566/1‑150 (MQ=255)
                                                                                    gCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCAc                                                               >  2:1425333/1‑150 (MQ=255)
                                                                                    gCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCAc                                                               >  2:1122079/1‑150 (MQ=255)
                                                                                       cGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCgg                                                            <  2:440415/150‑1 (MQ=255)
                                                                                          ttCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTCCGTCCACCGGAAc                                                         <  1:1139117/150‑1 (MQ=255)
                                                                                            cctccCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACgg                                                       <  2:627170/150‑1 (MQ=255)
                                                                                            cctccCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACgg                                                       <  2:1071815/150‑1 (MQ=255)
                                                                                            cctccCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACgg                                                       <  2:1119260/150‑1 (MQ=255)
                                                                                              tccCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTg                                                     <  1:822069/150‑1 (MQ=255)
                                                                                               cccGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGa                                                    >  2:1221918/1‑150 (MQ=255)
                                                                                                     cTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTgg                                              <  1:1425333/150‑1 (MQ=255)
                                                                                                     cTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTgg                                              <  1:1430493/150‑1 (MQ=255)
                                                                                                     cTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTgg                                              <  1:6212/150‑1 (MQ=255)
                                                                                                     cTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTgg                                              <  2:692504/150‑1 (MQ=255)
                                                                                                        aGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCAc                                           >  1:679592/1‑150 (MQ=255)
                                                                                                         gCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACg                                          >  1:1357282/1‑150 (MQ=255)
                                                                                                         gCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACg                                          >  2:1251003/1‑150 (MQ=255)
                                                                                                         gCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACg                                          <  2:344850/150‑1 (MQ=255)
                                                                                                          ccAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACgg                                         <  2:25724/150‑1 (MQ=255)
                                                                                                          ccAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACgg                                         >  1:684330/1‑150 (MQ=255)
                                                                                                           cAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGa                                        >  2:970779/1‑150 (MQ=255)
                                                                                                           cAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGa                                        >  1:763559/1‑150 (MQ=255)
                                                                                                              cGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACtt                                     >  1:538382/1‑150 (MQ=255)
                                                                                                               gACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTgg                                              <  1:720916/140‑1 (MQ=255)
                                                                                                               gACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACttt                                    <  2:636436/150‑1 (MQ=255)
                                                                                                               gACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCAAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTgg                                              >  2:720916/1‑140 (MQ=255)
                                                                                                                aCCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTa                                   >  1:689742/1‑150 (MQ=255)
                                                                                                                 ccAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTAc                                  >  1:1273066/1‑150 (MQ=255)
                                                                                                                 ccAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTAc                                  >  2:228071/1‑150 (MQ=255)
                                                                                                                   aGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCg                                <  1:197421/150‑1 (MQ=255)
                                                                                                                    gATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGt                               >  2:628408/1‑150 (MQ=255)
                                                                                                                         aggagTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAAcc                          <  2:689742/150‑1 (MQ=255)
                                                                                                                         aggagTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAAcc                          <  1:500896/150‑1 (MQ=255)
                                                                                                                         aggagTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAAcc                          <  2:1424187/150‑1 (MQ=255)
                                                                                                                             gTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGca                      >  1:1200108/1‑150 (MQ=255)
                                                                                                                               ttACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGcacc                    <  1:203529/150‑1 (MQ=255)
                                                                                                                                tACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGcacca                   >  1:1318063/1‑150 (MQ=255)
                                                                                                                                tACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGcacca                   >  2:1428926/1‑150 (MQ=255)
                                                                                                                                    tGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGAc               <  1:487826/150‑1 (MQ=255)
                                                                                                                                    tGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGAc               >  1:771363/1‑150 (MQ=255)
                                                                                                                                     tgTACGTAAAGTGATTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGAcc              >  1:245078/2‑150 (MQ=255)
                                                                                                                                     tgTACGTAAAGTGATTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGAcc              >  2:1453161/2‑150 (MQ=255)
                                                                                                                                     ggTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGAcc              >  1:492152/1‑150 (MQ=255)
                                                                                                                                      gTACGTAAAGTGATTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGTACCACGACCa             >  1:312957/1‑150 (MQ=255)
                                                                                                                                       tACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAg            >  1:1115797/1‑150 (MQ=255)
                                                                                                                                        aCGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGc           >  2:1011123/1‑150 (MQ=255)
                                                                                                                                         cGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCg          >  1:1041490/1‑150 (MQ=255)
                                                                                                                                           tAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTg        <  2:1244834/150‑1 (MQ=255)
                                                                                                                                            aaaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTAGCGTTAACCTGCACCACGTCCTGCGTgg       <  2:1361164/150‑1 (MQ=255)
                                                                                                                                             aaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCt                         >  2:940585/1‑131 (MQ=255)
                                                                                                                                             aaTCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCt                         <  1:940585/131‑1 (MQ=255)
                                                                                                                                                 aCTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCt  >  2:1437146/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACGCCATCAAGATATTTTGCTACCAGATGTTCCTGGCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCT  >  NC_000913/672099‑672393

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