Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,375,092 G→C 7.1% D44H (GAC→CAC)  ycjQ → putative zinc‑binding dehydrogenase YcjQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,375,0920GC 7.1% 382.8 / 14.4 113D44H (GAC→CAC) ycjQputative zinc‑binding dehydrogenase YcjQ
Reads supporting (aligned to +/- strand):  ref base G (53/52);  new base C (3/5);  total (56/57)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 117 total aligned reads displayed.

TTCTTTACAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCGACATTATCGAGTGCG  >  NC_000913/1374943‑1375239
                                                                                                                                                     |                                                                                                                                                   
ttCTTTACAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCg                                                                                                                                                     <  1:281701/150‑1 (MQ=255)
ttCTTTACAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCg                                                                                                                                                     <  2:766535/150‑1 (MQ=255)
 tCTTTACAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGa                                                                                                                                                    >  2:419439/1‑150 (MQ=255)
     tACAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTc                                                                                                                                                >  2:78343/1‑150 (MQ=255)
       cAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCcg                                                                                                                                              >  1:233182/1‑150 (MQ=255)
       cAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAACTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCcg                                                                                                                                              >  1:1368930/1‑150 (MQ=255)
         aGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCcgcg                                                                                                                                            <  2:1529323/150‑1 (MQ=255)
         aGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCcgcg                                                                                                                                            >  2:758569/1‑150 (MQ=255)
         aGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTCGACCACTTCcgtt                                                                                                                                            <  1:1095121/150‑3 (MQ=255)
          ggAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTCGACCACTTCcgttc                                                                                                                                           <  1:1192677/150‑4 (MQ=255)
            aGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTCGACCACTTCcgttccg                                                                                                                                         <  1:509004/150‑6 (MQ=255)
            aGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTCGACCACTTCcgttccg                                                                                                                                         <  1:1071303/150‑6 (MQ=255)
               tttAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTCGACCACTTCcgttc                                                                                                                                           <  1:559272/145‑4 (MQ=255)
               tttAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTCGACCACTTCcgttc                                                                                                                                           >  2:559272/1‑142 (MQ=255)
                   aaTGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGccc                                                                                                                                  >  1:1568802/1‑150 (MQ=255)
                        aaaaGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTa                                                                                                                             >  1:840800/1‑150 (MQ=255)
                          aaGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTAtt                                                                                                                           <  2:136650/150‑1 (MQ=255)
                          aaGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTAtt                                                                                                                           >  1:1160274/1‑150 (MQ=255)
                           aGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtg                                                                                                                          <  2:1632633/150‑1 (MQ=255)
                            tgtagtagCCACAGCACGGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAATTTAAGATCCGCGTTCGTTTCGGCGCACCGAAACAGGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCAGTTTATtga                                                                                                                         <  2:623121/148‑1 (MQ=255)
                              tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtgatg                                                                                                                       >  1:1582793/1‑150 (MQ=255)
                              tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtgatg                                                                                                                       <  2:716646/150‑1 (MQ=255)
                              tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtgatg                                                                                                                       <  1:553947/150‑1 (MQ=255)
                              tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtgatg                                                                                                                       <  1:103891/150‑1 (MQ=255)
                              tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtgatg                                                                                                                       <  1:390853/150‑1 (MQ=255)
                              tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATACGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtgatg                                                                                                                       <  1:764198/150‑1 (MQ=255)
                                gtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGaa                                                                                                                     <  1:1315743/150‑1 (MQ=255)
                                accagcgcaGCAACACGCGGTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGaa                                                                                                                     >  2:478605/8‑150 (MQ=255)
                                   gCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGAc                                                                                                                  <  2:342507/150‑1 (MQ=255)
                                   gCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGAc                                                                                                                  >  1:803142/1‑150 (MQ=255)
                                    ccACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACt                                                                                                                 <  2:461410/150‑1 (MQ=255)
                                     cacaGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACtt                                                                                                                >  1:1180224/1‑150 (MQ=255)
                                      acaGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACttt                                                                                                               <  1:32473/150‑1 (MQ=255)
                                        aGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTaa                                                                                                             <  1:168736/150‑1 (MQ=255)
                                            ccGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCg                                                                                                        >  2:85014/1‑150 (MQ=255)
                                              gcgTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGa                                                                                                       >  1:839943/1‑150 (MQ=255)
                                                 tgtTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATg                                                                                                    >  2:168486/1‑150 (MQ=255)
                                                    tgctgcGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCa                                                                                                 <  1:1524844/150‑1 (MQ=255)
                                                     gctgcGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAg                                                                                                <  1:758569/150‑1 (MQ=255)
                                                       tgcGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGAt                                                                                              >  1:288821/1‑150 (MQ=255)
                                                        gcgcTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATg                                                                                             >  2:69911/1‑150 (MQ=255)
                                                         cgcTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGAGAATGGCAGATGt                                                                                            >  2:1024132/1‑150 (MQ=255)
                                                              gTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGc                                                                                                                                    <  2:937250/105‑1 (MQ=255)
                                                              attGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGc                                                                                                                                    >  1:937250/2‑105 (MQ=255)
                                                                   gTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCcgcg                                                                                  <  1:83836/150‑1 (MQ=255)
                                                                     aTGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTc                                                                                >  2:239371/1‑150 (MQ=255)
                                                                        aaGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGATCGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCgg                                                                             <  2:1278588/150‑1 (MQ=255)
                                                                          gATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCa                                                                           >  1:63776/1‑150 (MQ=255)
                                                                          gATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCa                                                                           <  1:1247506/150‑1 (MQ=255)
                                                                           aTCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAg                                                                          <  1:1146946/150‑1 (MQ=255)
                                                                                 cGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATgcgc                                                                    <  2:1624866/150‑1 (MQ=255)
                                                                                       tAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCcgcgcg                                                              <  1:419439/150‑1 (MQ=255)
                                                                                        aGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCgg                                                             <  1:1552051/150‑1 (MQ=255)
                                                                                        aGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCgg                                                             <  2:839943/150‑1 (MQ=255)
                                                                                           tAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCAt                                                          >  2:776042/1‑150 (MQ=255)
                                                                                           tAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCAt                                                          >  2:785434/1‑150 (MQ=255)
                                                                                             aTGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTg                                                        <  1:1024132/150‑1 (MQ=255)
                                                                                                aagtgaagATCCGCGTGCGTGTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCGCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCAATCAGGGTCCGGCAGATGCTCGGCGCGTCATTTAGt                                                     <  2:386811/150‑1 (MQ=255)
                                                                                                 agtgaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGaa                                                                                                                     >  2:396272/1‑85 (MQ=255)
                                                                                                 agtgaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGaa                                                                                                                     <  1:396272/85‑1 (MQ=255)
                                                                                                 agtgaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGc                                                                            >  2:1390271/1‑126 (MQ=255)
                                                                                                 agtgaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGc                                                                            <  1:1390271/126‑1 (MQ=255)
                                                                                                  gtgaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGttt                                                   >  1:1402890/1‑150 (MQ=255)
                                                                                                   tgaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTg                                                  >  1:748843/1‑150 (MQ=255)
                                                                                                    gaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTgg                                                 >  1:1357072/1‑150 (MQ=255)
                                                                                                    gaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTgg                                                 >  1:434102/1‑150 (MQ=255)
                                                                                                     aagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTAACTCCGCGTCCGGCAGATGTGCCGCGCGGCATTGAGTTTGGc                                                >  2:1101115/1‑150 (MQ=255)
                                                                                                      agATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCa                                               >  1:1173310/1‑150 (MQ=255)
                                                                                                           gttccgtgtTTCGGTGCGCCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTc                                          >  1:272791/8‑150 (MQ=255)
                                                                                                            gcgtgcgtTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTcc                                         <  2:803142/150‑1 (MQ=255)
                                                                                                             cgtgcgtTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCa                                        <  1:548689/150‑1 (MQ=255)
                                                                                                               tgcgtTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGc                                      >  1:354157/1‑150 (MQ=255)
                                                                                                                 cgtTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTcc                                                                               <  2:711284/107‑1 (MQ=255)
                                                                                                                 cgtTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTcc                                                                               >  1:711284/1‑107 (MQ=255)
                                                                                                                 cgtTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCtt                                    <  2:443723/150‑1 (MQ=255)
                                                                                                                 cgtTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCtt                                    >  1:516994/1‑150 (MQ=255)
                                                                                                                  gtTTCGGCGCACCGAAACACGGAACGGAAGTGGGCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTg                                   <  1:277049/150‑1 (MQ=255)
                                                                                                                  gtTTCGGCGCACCGAAACACGGAACGGAAGTGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTgg                                  <  1:85014/150‑1 (MQ=255)
                                                                                                                    ttCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGc                                 <  2:84442/150‑1 (MQ=255)
                                                                                                                          gcACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACAtg                           >  2:470134/1‑150 (MQ=255)
                                                                                                                           cACCGAAACACGGAACGGAAGTGGTCGACTACCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATACCAGCTTGGCAACAagg                          >  2:507878/1‑150 (MQ=255)
                                                                                                                            aCCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCTAATTCCAGCTTGGCAACAtggt                         <  1:1061732/150‑1 (MQ=255)
                                                                                                                            aCCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACAtggt                         <  2:233182/150‑1 (MQ=255)
                                                                                                                            aCCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACAtggt                         <  2:655931/150‑1 (MQ=255)
                                                                                                                             ccGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACAtggtg                        >  2:961628/1‑150 (MQ=255)
                                                                                                                              cGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACAtggtgg                       >  2:1052429/1‑150 (MQ=255)
                                                                                                                                 aaCACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTg                    <  1:239371/150‑1 (MQ=255)
                                                                                                                                 aaCACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTg                    >  2:1404429/1‑150 (MQ=255)
                                                                                                                                   cacGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGc                  <  1:72968/150‑1 (MQ=255)
                                                                                                                                    acggaacggaAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCg                 <  2:1368930/150‑1 (MQ=255)
                                                                                                                                     cggaacggaaGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCGa                <  2:288821/150‑1 (MQ=255)
                                                                                                                                     cggaacggaaGTGGTCGACTTCCGCGCCGCCAGCCCGCTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCGa                >  1:1205206/1‑150 (MQ=255)
                                                                                                                                       gaacggaaGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCGACa              >  1:822104/1‑150 (MQ=255)
                                                                                                                                        aacggaaGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGCCGACAt             >  1:944855/1‑150 (MQ=255)
                                                                                                                                           ggaaGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCGACATTAt          >  2:388847/1‑150 (MQ=255)
                                                                                                                                           ggaaGTCGACCACTTCCGTTCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCGACATTAt          >  1:148957/1‑150 (MQ=255)
                                                                                                                                                 ggTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCGACATTATCGAGTg    <  1:1082402/150‑1 (MQ=255)
                                                                                                                                                  gTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCGACATTATCGAGTGc   <  2:1147983/150‑1 (MQ=255)
                                                                                                                                                   tCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCGACATTATCGAGTGCg  >  2:613183/1‑150 (MQ=255)
                                                                                                                                                   tCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCGACATTATCGAGTGCg  <  2:272791/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTCTTTACAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATGGTGGTTGGCGACATTATCGAGTGCG  >  NC_000913/1374943‑1375239

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