Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,020,777 T→G 5.4% G184G (GGT→GGG ubiB → ubiquinone biosynthesis protein UbiB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,020,7770TG 5.4% 252.9 / 4.8 88G184G (GGT→GGGubiBubiquinone biosynthesis protein UbiB
Reads supporting (aligned to +/- strand):  ref base T (43/40);  new base G (3/2);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CGCCCAGGTTCATACCGCGCGATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAGTTTACCCTG  >  NC_000913/4020630‑4020925
                                                                                                                                                   |                                                                                                                                                    
cGCCCAGGTTCATACCGCGCGATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCg                                                                                                                                                    >  1:129953/1‑150 (MQ=255)
 gCCCAGGTTCATACCGCGCGATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGc                                                                                                                                                   >  1:518722/1‑150 (MQ=255)
 gCCCAGGTTCATACCGCGCGATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCAACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGc                                                                                                                                                   >  2:1463489/1‑150 (MQ=255)
 gCCCAGGTTCATACCGCGCGATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATATCCCGTATGGAACCATGGGTGCCGCGTTTGATGCCGGATGTTCGc                                                                                                                                                   >  1:385451/1‑150 (MQ=255)
   ccAGGTTCATACCGCGCGATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCg                                                                                                                                                 <  2:993226/150‑1 (MQ=255)
     aGGTTCATACCGCGCGATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTc                                                                                                                                               <  2:617167/150‑1 (MQ=255)
      ggTTCATACCGCGCGATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCt                                                                                                                                              <  2:471826/150‑1 (MQ=255)
            tACCGCGCGATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGccc                                                                                                                                        >  2:1325500/1‑150 (MQ=255)
             aCCGCGCGATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCa                                                                                                                                       <  2:1468041/150‑1 (MQ=255)
               cgcgcgATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAAc                                                                                                                                     >  1:74437/1‑150 (MQ=255)
               cccgcgATTGAAATCGTACGCTTAAGAGGTGGTGATTAATGTCATCCGCGCGGTTATTTTGCCGGCTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAAc                                                                                                                                     <  2:1452360/148‑1 (MQ=255)
                gcgcgATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAAcc                                                                                                                                    >  1:956769/1‑150 (MQ=255)
                gcccgATTCAAATCGAATGGTATAGAGTTCGTCATTTAACTCATCCGCCCCGATTTTTTTCCGGTATTTAACGCGCATCTGAAACTTATCTACCGACTGGCTCGCACGGTGCCGCGTATGCTGCCGGATGGTCGCCGACTGCGCCCAAcc                                                                                                                                    >  2:2906/4‑150 (MQ=255)
                     aTTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAgt                                                                                                                               <  1:214608/150‑1 (MQ=255)
                       tGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAgtgg                                                                                                                             >  1:836926/1‑150 (MQ=255)
                         aaaTCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAgtggtg                                                                                                                           <  1:1345286/150‑1 (MQ=255)
                             cGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGcgcg                                                                                                                       <  2:1456957/150‑1 (MQ=255)
                              gAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGGGGTGCGCGa                                                                                                                      >  2:240086/1‑150 (MQ=255)
                               aaTGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAg                                                                                                                     <  1:1311077/150‑1 (MQ=255)
                               aaTGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAg                                                                                                                     <  1:514138/150‑1 (MQ=255)
                                 tGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTa                                                                                                                   <  1:1117706/150‑1 (MQ=255)
                                      aaGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGaaa                                                                                                              <  2:956769/150‑1 (MQ=255)
                                          ggtggtGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGAc                                                                                                          >  2:1278975/1‑150 (MQ=255)
                                          ggtggtGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGAc                                                                                                          >  2:378340/1‑150 (MQ=255)
                                          ggtggtGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGAc                                                                                                          <  2:518722/150‑1 (MQ=255)
                                          ggtggtGATTAAAGTCATCCGCACGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGAc                                                                                                          <  2:673548/150‑1 (MQ=255)
                                            tggtgATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACat                                                                                                        >  1:642575/1‑150 (MQ=255)
                                              gtgATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAAGTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACattg                                                                                                      >  2:372257/1‑150 (MQ=255)
                                               tgATTAAAGTCATCCGCCCGGATATTTTGCAGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATCGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACattga                                                                                                     >  2:331779/1‑150 (MQ=255)
                                                gATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACattgat                                                                                                    >  1:1396049/1‑150 (MQ=255)
                                                 aTTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACattgatt                                                                                                   <  2:920369/150‑1 (MQ=255)
                                                  ttAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACattgattg                                                                                                  >  1:1183100/1‑150 (MQ=255)
                                                   tAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACAttgattga                                                                                                 <  1:79215/150‑1 (MQ=255)
                                                    aaaGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGTATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATtgattgat                                                                                                <  2:1021448/150‑1 (MQ=255)
                                                     aaGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATtgatg                                                                                               >  1:114388/1‑150 (MQ=255)
                                                         cATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACt                                                                                           <  1:1555847/150‑1 (MQ=255)
                                                         cATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACt                                                                                           >  1:385384/1‑150 (MQ=255)
                                                         cATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACt                                                                                           >  2:1055663/1‑150 (MQ=255)
                                                          aTCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGGCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTg                                                                                          <  2:841897/150‑1 (MQ=255)
                                                             cGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCTACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAAt                                                                                       <  2:1592831/150‑1 (MQ=255)
                                                             cGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAAt                                                                                       >  2:593185/1‑150 (MQ=255)
                                                             cGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAAt                                                                                       >  1:344955/1‑150 (MQ=255)
                                                               cccGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAAttt                                                                                     >  2:1043130/1‑150 (MQ=255)
                                                               cccGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAAttt                                                                                     <  2:788761/150‑1 (MQ=255)
                                                               cccGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAAttt                                                                                     <  2:794620/150‑1 (MQ=255)
                                                                 cGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTtgc                                                                                   <  1:1325500/150‑1 (MQ=255)
                                                                    atagTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGGGCCGCGTTTGCTGCCGGATGGTCGGCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTtgctgc                                                                                <  2:1608829/150‑1 (MQ=255)
                                                                       ttttGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCggg                                                                             <  2:1411847/150‑1 (MQ=255)
                                                                         ttGCCGGTTATTAAAGCGGAGCTGAAACTTATCTACCGTCTGGCTCGCTGGGGGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGaa                                                                           <  2:1598641/150‑1 (MQ=255)
                                                                          tGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGTGGGAAt                                                                          >  2:537851/1‑150 (MQ=255)
                                                                             cGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTg                                                                       <  2:476194/150‑1 (MQ=255)
                                                                                   ttaAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACg                                                                 >  1:20900/1‑150 (MQ=255)
                                                                                   ttaAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACg                                                                 <  1:340938/150‑1 (MQ=255)
                                                                                     aaaGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGcc                                                               >  2:1353121/1‑150 (MQ=255)
                                                                                     aaaGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGcc                                                               >  1:1322533/1‑150 (MQ=255)
                                                                                     aaaGCGGATCTGAAACTTATCTACCGTATGGCTCGCTGGGTGCCGAGTTTGCTGCCGGATGGTCGCCGTCTGCGACCAACAGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGcc                                                               >  2:64394/1‑150 (MQ=255)
                                                                                       aGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCAt                                                             >  2:1507915/1‑150 (MQ=255)
                                                                                       aGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCAt                                                             <  1:965965/150‑1 (MQ=255)
                                                                                        gCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCAtt                                                            <  1:331779/150‑1 (MQ=255)
                                                                                          ggATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCa                                                          >  2:1626591/1‑150 (MQ=255)
                                                                                            aTCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGc                                                        <  1:378340/150‑1 (MQ=255)
                                                                                               tGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTc                                                           >  1:987083/1‑144 (MQ=255)
                                                                                               tGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTc                                                           <  2:987083/144‑1 (MQ=255)
                                                                                                 aaaCTTATCTACCGTCTGGCGCGCTGGGCGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAACTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCgg                                                   <  2:265642/150‑1 (MQ=255)
                                                                                                   aCTTATCAACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGgcg                                                 <  1:1463489/150‑1 (MQ=255)
                                                                                                              gTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAg                                      <  2:1396049/150‑1 (MQ=255)
                                                                                                                  ggCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAg                                  >  1:400748/1‑150 (MQ=255)
                                                                                                                    cTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGcc                                >  1:429224/1‑150 (MQ=255)
                                                                                                                    cTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCAGCGGCGCAATGTTGAAGACAGcc                                >  1:1427604/1‑150 (MQ=255)
                                                                                                                      cGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCg                              <  2:344955/150‑1 (MQ=255)
                                                                                                                        cTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGAt                            >  1:1352529/1‑150 (MQ=255)
                                                                                                                           ggTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTtgctgc                                                                                <  2:1412156/95‑1 (MQ=255)
                                                                                                                           ggTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTtgctgc                                                                                >  1:1412156/1‑95 (MQ=255)
                                                                                                                            gTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGctc                        >  2:1164070/1‑150 (MQ=255)
                                                                                                                            gTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGctc                        >  1:775521/1‑150 (MQ=255)
                                                                                                                                cgcgTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACa                    >  1:266645/1‑150 (MQ=255)
                                                                                                                                       gctgcCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTAc                     <  2:1132213/142‑1 (MQ=255)
                                                                                                                                       gctgcCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTAc                     >  1:1132213/1‑142 (MQ=255)
                                                                                                                                         tgcCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAg           <  1:1507915/150‑1 (MQ=255)
                                                                                                                                         tgcCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAg           <  2:642575/150‑1 (MQ=255)
                                                                                                                                         tgcCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAg           <  2:129953/150‑1 (MQ=255)
                                                                                                                                         cttcGGTTGGGCGCAGACGGCGACCATCCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAg           <  2:265639/147‑1 (MQ=255)
                                                                                                                                         cttcGGTTGGGCGCAGACGGCGACCATCCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAg           >  2:1420534/4‑150 (MQ=255)
                                                                                                                                         cttcGGTTGGGCGCAGACGGCGACCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAg           >  1:780548/4‑150 (MQ=255)
                                                                                                                                         cttcGGTTGGGCGCAGACGGCGACCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAg           >  1:774969/4‑150 (MQ=255)
                                                                                                                                                  gTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAGTTTACCCTg  >  2:1444477/1‑150 (MQ=255)
                                                                                                                                                  gTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAGTTTACCCTg  >  2:303458/1‑150 (MQ=255)
                                                                                                                                                  gTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAGTTTACCCTg  >  2:784261/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCCCAGGTTCATACCGCGCGATTGAAATCGAATGGTAAAGAGGTGGTGATTAAAGTCATCCGCCCGGATATTTTGCCGGTTATTAAAGCGGATCTGAAACTTATCTACCGTCTGGCTCGCTGGGTGCCGCGTTTGCTGCCGGATGGTCGCCGTCTGCGCCCAACCGAAGTGGTGCGCGAGTACGAAAAGACATTGATTGATGAACTGAATTTGCTGCGGGAATCTGCCAACGCCATTCAGCTTCGGCGCAATTTTGAAGACAGCCCGATGCTCTACATCCCGGAAGTTTACCCTG  >  NC_000913/4020630‑4020925

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: