Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,967,705 T→A 10.0% intergenic (+29/‑211) rho → / → rfe transcription termination factor Rho/UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,7050TA10.0% 316.1 / 20.3 110intergenic (+29/‑211)rho/rfetranscription termination factor Rho/UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase
Reads supporting (aligned to +/- strand):  ref base T (46/53);  new base A (3/8);  total (49/61)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 116 total aligned reads displayed.

GTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCC  >  NC_000913/3967556‑3967853
                                                                                                                                                     |                                                                                                                                                    
gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgtgt                                                                                                                                                      >  1:530505/1‑150 (MQ=255)
 tgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGtt                                                                                                                                                     >  2:1150478/1‑150 (MQ=255)
  ggATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGttt                                                                                                                                                    >  1:318699/1‑150 (MQ=255)
   gATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTa                                                                                                                                                   >  1:743279/1‑150 (MQ=255)
   gATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTa                                                                                                                                                   <  1:533989/150‑1 (MQ=255)
     tCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACg                                                                                                                                                 <  1:636815/150‑1 (MQ=255)
      ccTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGt                                                                                                                                                <  2:29203/150‑1 (MQ=255)
      ccTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGt                                                                                                                                                >  1:1324210/1‑150 (MQ=255)
        tGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTgg                                                                                                                                              <  2:695668/150‑1 (MQ=255)
         gcgcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGc                                                                                                                                             >  1:1413689/1‑150 (MQ=255)
          cgcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTAAGTGGCg                                                                                                                                            >  2:1333944/1‑150 (MQ=255)
           gcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGt                                                                                                                                           <  1:1030178/150‑1 (MQ=255)
            cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCGtt                                                                                                                                          <  2:1051146/150‑1 (MQ=255)
             aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCGttt                                                                                                                                         <  1:1308189/150‑1 (MQ=255)
             aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCGttt                                                                                                                                         <  1:633995/150‑1 (MQ=255)
               aaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTg                                                                                                                                       <  2:155379/150‑1 (MQ=255)
               aaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCGttttt                                                                                                                                       <  1:879052/150‑2 (MQ=255)
               aaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGGGGCGttttt                                                                                                                                       <  2:604693/150‑2 (MQ=255)
                aTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGc                                                                                                                                      <  1:112219/150‑1 (MQ=255)
                aTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCGtttttg                                                                                                                                      <  1:1225961/150‑3 (MQ=255)
                    ttCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCtttt                                                                                                                                  >  1:1032953/1‑150 (MQ=255)
                    ttCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCtttt                                                                                                                                  >  2:107121/1‑150 (MQ=255)
                    ttCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCGtttttggcat                                                                                                                                  <  2:1079087/150‑7 (MQ=255)
                     tCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCttttt                                                                                                                                 >  2:833260/1‑149 (MQ=255)
                       aCCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTtata                                                                                                                               <  1:1040484/150‑1 (MQ=255)
                          cGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCt                                                                                                                            >  1:1246476/1‑150 (MQ=255)
                           gATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTg                                                                                                                           >  1:1155641/1‑150 (MQ=255)
                            aTGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGt                                                                                                                          >  1:207483/1‑150 (MQ=255)
                               ggCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAAt                                                                                                                       <  2:1103929/150‑1 (MQ=255)
                                 cGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCGATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCt                                                                                                                     <  2:1206053/150‑1 (MQ=255)
                                 cGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCt                                                                                                                     <  1:1328201/150‑1 (MQ=255)
                                      tCGTTGCAATTGAATTCCTCCTTAATAAACTTGCAATGAACAATACCAATGACGATTTCTTCGAAATGTTTTAACGCTCATAAATTTGGCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATg                                                                                                                <  2:870647/150‑1 (MQ=255)
                                       cGATGCAATGGAATTCCTCATTAATAACCTGGCAATGACCAAGACGAAGGACGATTTCTTCGAAATGATGAAACGCTGATAAATTTGTCTTATGCCAAAAACGCCACGAGTTTACGTGGGGTTTTGCTTTTATATCTGTAATCTTAATGc                                                                                                               >  2:865324/1‑150 (MQ=255)
                                        gATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGcc                                                                                                              >  1:324873/1‑150 (MQ=255)
                                        gATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGcc                                                                                                              >  2:340827/1‑150 (MQ=255)
                                         aTGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCcg                                                                                                             >  2:440486/1‑150 (MQ=255)
                                           gCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCcgcg                                                                                                           <  1:800532/150‑1 (MQ=255)
                                             aaTGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATTAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATTTCTGTAATCTTAATGCCgcgcg                                                                                                         >  2:1383814/1‑149 (MQ=255)
                                             aaTGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCt                                                                                                         >  1:355208/1‑150 (MQ=255)
                                                ggAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTTCGTTGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGc                                                                                                      <  2:928684/150‑1 (MQ=255)
                                                ggAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGc                                                                                                      <  1:1045886/150‑1 (MQ=255)
                                                    ttCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCGCGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATg                                                                                                  <  2:1147587/150‑1 (MQ=255)
                                                      ccTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGtt                                                                                                <  2:1483712/150‑1 (MQ=255)
                                                       cTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTa                                                                                               >  1:1436051/1‑150 (MQ=255)
                                                         cATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAgg                                                                                             <  1:675703/150‑1 (MQ=255)
                                                         cATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAgg                                                                                             >  2:518405/1‑150 (MQ=255)
                                                         cATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTAAACACGTGGCGtttttgg                                                                                                                                     <  2:894495/110‑4 (MQ=38)
                                                               taaACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAAtt                                                                                       >  1:291297/1‑150 (MQ=255)
                                                                  aCTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCt                                                                                    >  1:146715/1‑150 (MQ=255)
                                                                   cTGGCAATGACCAAGTCCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTg                                                                                   <  2:73510/150‑1 (MQ=255)
                                                                   cTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTg                                                                                   <  1:1241092/150‑1 (MQ=255)
                                                                    tGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTgg                                                                                  <  2:955803/150‑1 (MQ=255)
                                                                      gCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGaa                                                                                >  1:164099/1‑150 (MQ=255)
                                                                         aTGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAAttt                                                                             >  1:19668/1‑150 (MQ=255)
                                                                         aTGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAAttt                                                                             >  2:1159854/1‑150 (MQ=255)
                                                                           gaccaagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGATTTTGc                                                                           <  1:440486/150‑1 (MQ=255)
                                                                            accaagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCt                                                                          <  1:205535/150‑1 (MQ=255)
                                                                             ccaagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTg                                                                         <  1:1389169/150‑1 (MQ=255)
                                                                             ccaagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTg                                                                         >  2:616570/1‑150 (MQ=255)
                                                                               aagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGc                                                                       <  2:445741/150‑1 (MQ=255)
                                                                                   ccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGt                                                                   >  2:681036/1‑150 (MQ=255)
                                                                                      aTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTAt                                                                <  1:926000/150‑1 (MQ=255)
                                                                                        gaagATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCa                                                             >  1:633913/4‑150 (MQ=255)
                                                                                         aCGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCa                                                             >  1:370657/1‑150 (MQ=255)
                                                                                         aCGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCa                                                             <  2:69431/150‑1 (MQ=255)
                                                                                            aTTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAAtt                                                          <  2:530505/150‑1 (MQ=255)
                                                                                            aTTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAAtt                                                          <  2:1032953/150‑1 (MQ=255)
                                                                                              ttcttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTg                                                        <  1:930911/150‑1 (MQ=255)
                                                                                              ttcttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGc                                                       <  2:633913/150‑1 (MQ=255)
                                                                                               tcttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGc                                                       <  1:1159854/150‑1 (MQ=255)
                                                                                               tcttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGc                                                       <  1:616570/150‑1 (MQ=255)
                                                                                                 ttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCat                                                     <  2:633053/150‑1 (MQ=255)
                                                                                                 ttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCat                                                     >  2:1031019/1‑150 (MQ=255)
                                                                                                  tcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCata                                                    <  1:168909/150‑1 (MQ=255)
                                                                                                   cGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCatat                                                   >  1:824011/1‑150 (MQ=255)
                                                                                                    gAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTAAGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATc                                                  <  1:1333944/150‑1 (MQ=255)
                                                                                                      aaTGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCaa                                                >  1:753204/1‑150 (MQ=255)
                                                                                                       atgatgAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCAAGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCaaa                                               >  2:1002708/1‑150 (MQ=255)
                                                                                                         gatgaAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATg                                             >  2:1140719/1‑150 (MQ=255)
                                                                                                                 gCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAAttt                                     >  2:1206780/1‑150 (MQ=255)
                                                                                                                             gTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACt                         <  1:107121/150‑1 (MQ=255)
                                                                                                                             gTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACt                         >  1:1490586/1‑150 (MQ=255)
                                                                                                                             gTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACt                         >  2:928435/1‑150 (MQ=255)
                                                                                                                                 tATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACtggtg                     <  1:518405/150‑1 (MQ=255)
                                                                                                                                   tGCCAAAACCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGAAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTggg                   <  2:1158552/150‑1 (MQ=255)
                                                                                                                                    gCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGt                  >  2:1492998/1‑150 (MQ=255)
                                                                                                                                    gCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGt                  <  2:1277199/150‑1 (MQ=255)
                                                                                                                                     ccAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGtt                 <  2:19668/150‑1 (MQ=255)
                                                                                                                                      cAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGttt                <  2:479226/150‑1 (MQ=255)
                                                                                                                                       caaaaCGCCACGTAAACACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTg               >  2:593949/2‑150 (MQ=255)
                                                                                                                                        aaaaCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTgg              <  2:1111654/150‑1 (MQ=255)
                                                                                                                                         aaaCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGa             <  1:749054/150‑1 (MQ=255)
                                                                                                                                         aaaCGCCACGTAAACACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGa             >  2:61793/1‑150 (MQ=255)
                                                                                                                                            cGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACg          <  1:464386/150‑1 (MQ=255)
                                                                                                                                            cGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACg          <  1:634372/150‑1 (MQ=255)
                                                                                                                                             gCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACgg         >  1:464952/1‑150 (MQ=255)
                                                                                                                                              ccACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGa        >  1:1109376/1‑150 (MQ=255)
                                                                                                                                                   tgtTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTc   <  2:531325/150‑1 (MQ=255)
                                                                                                                                                   tgtTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTc   >  1:104537/1‑150 (MQ=255)
                                                                                                                                                    gtTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTcc  >  2:815555/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCC  >  NC_000913/3967556‑3967853

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: