Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,257,223 T→A 7.6% F37I (TTC→ATC) ‡ lexA → transcriptional repressor of SOS regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,257,2230TA 7.6% 199.0 / 6.4 79F37I (TTC→ATC) ‡lexAtranscriptional repressor of SOS regulon
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base A (4/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

ATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCGTGTGGCTGCCGGTGA  >  NC_000913/4257076‑4257371
                                                                                                                                                   |                                                                                                                                                     
ataCTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGgttc                                                                                                                                                     <  1:553035/150‑1 (MQ=255)
ataCTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGgttc                                                                                                                                                     <  1:1280634/150‑1 (MQ=255)
 taCTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGgttcc                                                                                                                                                    <  1:606999/150‑1 (MQ=255)
 taCTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGAATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGgttcc                                                                                                                                                    >  1:111761/1‑150 (MQ=255)
     cacaGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGgttccgttc                                                                                                                                                >  2:1118883/1‑150 (MQ=255)
     cacaGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGgttccgttc                                                                                                                                                >  1:553050/1‑150 (MQ=255)
       caGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTccc                                                                                                                                              <  2:1167856/150‑1 (MQ=255)
         gCATAACTGTATATACACCCAGGGGGCTGAATGAAAGCGTTAACGGCCAGGCAACAAGATGTGTTTGATCTCATCCGTGATAACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCa                                                                                                                                            <  2:1305395/150‑1 (MQ=255)
         gCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCa                                                                                                                                            <  2:1310501/150‑1 (MQ=255)
            tAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAAc                                                                                                                                         >  1:841414/1‑150 (MQ=255)
            aaaCTGTATATACACCCAGGGGGCGTAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGAACGGAACCCCAAAc                                                                                                                                         <  1:387224/149‑1 (MQ=255)
              aCTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGATATCGCGCAGCGTTTGGGGTTCCGTTCCCCAACcgc                                                                                                                                       >  1:225466/1‑150 (MQ=255)
              aCTGTATAAACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAAcgc                                                                                                                                       >  2:493705/1‑150 (MQ=255)
                tGTATAAACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGAACGGAACCCCAAAcgctg                                                                                                                                     <  1:866449/150‑3 (MQ=255)
                 gTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGc                                                                                                                                    <  1:624053/150‑1 (MQ=255)
                    aataCACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTga                                                                                                                                 <  1:929090/149‑1 (MQ=255)
                       acacCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTgaaga                                                                                                                              >  1:708011/1‑150 (MQ=255)
                                   cGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCACACGCGGCTGAAGAACATCTGAAGgc                                                                                                                  >  2:895535/1‑150 (MQ=255)
                                      aaTGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCt                                                                                                               <  2:722763/150‑1 (MQ=255)
                                        tGAAAGCTTTAACGGCCAGGCTACAAGAGGGGGTTGATCTCATCCGTTATCACAGCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGGTCCGTTCCCCAAACGCGGCTGAAGATCATCTGAAGGCGCTgg                                                                                                             <  1:708265/150‑1 (MQ=255)
                                            aGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACg                                                                                                         <  2:391133/150‑1 (MQ=255)
                                            aGCGTCAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACg                                                                                                         >  1:692533/1‑150 (MQ=255)
                                             gCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGc                                                                                                        >  1:1095567/1‑150 (MQ=255)
                                               gTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCCAACGCGGCTGAAGACCATCTGAAGGCGCTGGCACGCaa                                                                                                      >  2:967819/1‑150 (MQ=255)
                                                    cGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCg                                                                                                 <  1:9818/150‑1 (MQ=255)
                                                      gCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCACACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGtt                                                                                               >  1:180817/1‑150 (MQ=255)
                                                      gCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACTCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGtt                                                                                               >  2:994152/1‑150 (MQ=255)
                                                        cAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTAt                                                                                             >  1:86226/1‑150 (MQ=255)
                                                          ggCAACAAGAGGTGTTTGATCTCATCCGTTATCACATCAGCCAGACAGGAATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCACAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTg                                                                                           <  1:831142/150‑1 (MQ=255)
                                                          ggCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCAGAAGGCGCTGGCCCGCAAAGGCGTTATTg                                                                                           >  1:1240248/1‑150 (MQ=255)
                                                           gCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGa                                                                                          >  1:67091/1‑150 (MQ=255)
                                                            caacaaGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGaa                                                                                         >  2:125631/1‑150 (MQ=255)
                                                            caacaaGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCACAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGaa                                                                                         >  2:1051332/1‑150 (MQ=255)
                                                                 agagGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑t                                                                                   >  1:417203/1‑150 (MQ=255)
                                                                     gtgtTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTcc                                                                               <  2:553050/150‑1 (MQ=255)
                                                                     gtgtTTGATCTCATCCGTGATCAAATCAGCCAGACATGTTAGCCGCCGACGCGTGCTGAAATCGCGCATCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTcc                                                                               <  2:368162/150‑1 (MQ=255)
                                                                       gtTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCgg                                                                             >  2:1324156/1‑150 (MQ=255)
                                                                       gtTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCgg                                                                             >  1:1103061/1‑150 (MQ=255)
                                                                       gtTTGATCTCATCCGTGATCACATCAGCCAAACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCgg                                                                             <  2:892835/150‑1 (MQ=255)
                                                                        tttGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGgc                                                                            <  2:1266329/150‑1 (MQ=255)
                                                                        tttGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGgc                                                                            <  2:121532/150‑1 (MQ=255)
                                                                          tGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTGTTTTCCGgcg                                                                           >  2:65004/1‑150 (MQ=255)
                                                                          tGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGgcgc                                                                          >  2:1123278/1‑150 (MQ=255)
                                                                             tctcATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATc                                                                       <  1:1118883/150‑1 (MQ=255)
                                                                                 aTCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCAcgc                                                                   >  2:487170/1‑150 (MQ=255)
                                                                                   ccGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCgg                                                                 <  1:497362/150‑1 (MQ=255)
                                                                                    cGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCggg                                                                <  2:52087/150‑1 (MQ=255)
                                                                                          cacaTCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGt                                                          >  1:752780/1‑150 (MQ=255)
                                                                                           acaTCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTc                                                         <  2:1103061/150‑1 (MQ=255)
                                                                                           acaTCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTc                                                         >  2:469802/1‑150 (MQ=255)
                                                                                            caTCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTTTTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCt                                                        >  1:928723/1‑150 (MQ=255)
                                                                                             aTCAGCGAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCAGGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTg                                                       <  1:255218/150‑1 (MQ=255)
                                                                                              tCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGt                                                      <  2:441706/150‑1 (MQ=255)
                                                                                              tCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGt                                                      <  2:718836/150‑1 (MQ=255)
                                                                                                aGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTg                                                    <  2:493862/150‑1 (MQ=255)
                                                                                                aGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTg                                                    >  2:885464/1‑150 (MQ=255)
                                                                                                    aGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCagg                                                <  2:559010/150‑1 (MQ=255)
                                                                                                     gACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCagga                                               >  2:153370/1‑150 (MQ=255)
                                                                                                     gACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCagga                                               >  2:966404/1‑150 (MQ=255)
                                                                                                       cAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCaggaag                                             >  1:554372/1‑150 (MQ=255)
                                                                                                           tATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCaggaagagga                                         <  1:1254871/150‑1 (MQ=255)
                                                                                                             tGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACTCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCaggaagaggaag                                       <  1:994152/150‑1 (MQ=255)
                                                                                                                cgccgACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACCTCTGAAGGCGCTGGCACGCAAAGGGGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGaagaag                                    >  2:1234257/1‑150 (MQ=255)
                                                                                                                      cgcgTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGc                              <  1:446120/150‑1 (MQ=255)
                                                                                                                       gcgTGCGGAAATCGCGCGGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGcc                             >  1:1253791/1‑150 (MQ=255)
                                                                                                                        cgTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCg                            >  1:564307/1‑150 (MQ=255)
                                                                                                                        cgTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCg                            >  2:61000/1‑150 (MQ=255)
                                                                                                                            cGGAAATCGCTCAGCGTTTGGGTTTCCGTTCACCAAACGCGGCTGAAGAACATCTGAAGGCGCTTGCACGCAAAGGCGTTATTGAAATTG‑TTTCCTGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGATGAAGAAGGGTTGCAGCTgg                        <  2:1010934/150‑1 (MQ=255)
                                                                                                                              gAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCGGAAGAACAGCTGAAGGCGCTGGCGCGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATGCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTa                      >  1:175608/1‑150 (MQ=255)
                                                                                                                                 aTCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGt                   <  1:125173/150‑1 (MQ=255)
                                                                                                                                 aTCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGt                   <  1:889409/150‑1 (MQ=255)
                                                                                                                                   cgcgcgGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCg                 <  2:1253791/150‑1 (MQ=255)
                                                                                                                                    gcgcAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCGGAAGAACAGCTGAAGGTGCTGGCACGCAAAGGCGGTATTGAAATTG‑TTTCCGGCGCATCACGCGGGAGTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGGAGGTCgt                >  1:121660/1‑150 (MQ=255)
                                                                                                                                     agccgCGTTTGGGGAACGGAACCCCAAACGCTGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCgtg               >  1:339680/5‑150 (MQ=255)
                                                                                                                                      gccgCGTTTGGGGAACGGAACCCCAAACGCTGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCgtgt              >  2:22645/4‑150 (MQ=255)
                                                                                                                                      gccgCGTTTGGGGAACGGAACCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCgtgt              >  2:404515/4‑150 (MQ=255)
                                                                                                                                        aGCGGTTGGGGGTCCGTTCCCCAATCGGGGCGGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTCTTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCGTGTgg            <  1:126021/150‑1 (MQ=255)
                                                                                                                                         gCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCGTGTGGc           <  2:374766/150‑1 (MQ=255)
                                                                                                                                          cGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCGTGTGGCt          <  2:708011/150‑1 (MQ=255)
                                                                                                                                          cGTTTGGGGAACGGAACCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCGTGTGGCt          >  1:1343367/1‑150 (MQ=255)
                                                                                                                                            tttGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCGTGTGGCTGc        <  2:180817/150‑1 (MQ=255)
                                                                                                                                             ttGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCGTGTGGCTGcc       <  2:1341728/150‑1 (MQ=255)
                                                                                                                                              tGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAATTTG‑TTTCCGGCGCATCACGCGGGCTTCGTCTGTTGCAGGAAGAGGAAGAAGTGTTGCCGCTGGTAGGTCGTGTGGCTGCCg      >  2:643386/1‑150 (MQ=255)
                                                                                                                                              tGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCGTGTGGCTGCCg      >  1:1282313/1‑150 (MQ=255)
                                                                                                                                                gggTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGCTTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCGTGTGGCTGCCGGt    >  1:375534/1‑150 (MQ=255)
                                                                                                                                                  gttccgttacCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGACGGGTTGCCGCTGGTAGGTCGTGTGGCTGCCGGTGa  >  1:190543/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGCCGCCGACGCGTGCGGAAATCGCGCAGCGTTTGGGGTTCCGTTCCCCAAACGCGGCTGAAGAACATCTGAAGGCGCTGGCACGCAAAGGCGTTATTGAAATTG‑TTTCCGGCGCATCACGCGGGATTCGTCTGTTGCAGGAAGAGGAAGAAGGGTTGCCGCTGGTAGGTCGTGTGGCTGCCGGTGA  >  NC_000913/4257076‑4257371

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