Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,264 A→C 5.2% intergenic (+34/‑183) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2640AC 5.2% 1594.5 / 59.8 480intergenic (+34/‑183)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (221/233);  new base C (13/12);  total (234/246)
Fisher's exact test for biased strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 487 total aligned reads displayed.

GTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGA  >  NC_000913/4060116‑4060412
                                                                                                                                                    |                                                                                                                                                     
     tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                                                                 >  1:5290790/1‑150 (MQ=255)
     tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                                                                 >  1:1883559/1‑150 (MQ=255)
       gatgatGACGAACTGGTAGAAGTGACTCCGCCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑t                                                                                                                                              <  1:3576843/150‑1 (MQ=255)
        atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tt                                                                                                                                             >  2:4695262/1‑150 (MQ=255)
         tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttt                                                                                                                                            <  2:2694708/150‑1 (MQ=255)
         tgaggaCGAACTGGGAGAAGTGACTCCGACCTCTCTCCGTATTCGTAAACGTCACCGGCCGGAACACGATCGTCGCCGCGCCAACCGCGCACCGAACGACGATTACTTTCGTTTTTTAGTCATAAAAAACCTGCCCGCGATGGCAGG‑ttt                                                                                                                                            <  2:3219190/150‑1 (MQ=255)
            tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttt                                                                                                                                         <  1:3224146/150‑1 (MQ=255)
            tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttt                                                                                                                                         <  1:2621128/150‑1 (MQ=255)
             gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑tttttta                                                                                                                                        <  2:749847/150‑2 (MQ=255)
             gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttt                                                                                                                                        <  1:6691516/150‑1 (MQ=255)
             gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttt                                                                                                                                        >  2:5948805/1‑150 (MQ=255)
               cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCTAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                                                                      <  1:2106462/150‑1 (MQ=255)
               cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                                                                      <  2:832136/150‑1 (MQ=255)
               cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                                                                      <  2:2444151/150‑1 (MQ=255)
                gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTg                                                                                                                                     <  1:1119171/150‑1 (MQ=255)
                 aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGc                                                                                                                                    <  2:4130036/150‑1 (MQ=255)
                  aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGcc                                                                                                                                   >  2:3596725/1‑150 (MQ=255)
                  aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTGTTTGcc                                                                                                                                   >  2:4794402/1‑150 (MQ=255)
                      gTAGAAGTGACTCCGACCTCTATCCGTATTGGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTc                                                                                                                               >  1:3706125/1‑150 (MQ=255)
                      gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTc                                                                                                                               >  1:5829657/1‑150 (MQ=255)
                       tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTct                                                                                                                              >  1:2935216/1‑149 (MQ=255)
                       tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCg                                                                                                                              >  1:842696/1‑150 (MQ=255)
                       tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCg                                                                                                                              >  2:628087/1‑150 (MQ=255)
                          aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGa                                                                                                                           >  2:544058/1‑150 (MQ=255)
                            gTGACTCCGACCTCTATCCGGATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAgc                                                                                                                         <  2:3798269/150‑1 (MQ=255)
                               aCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCTca                                                                                                                      <  1:2825661/150‑1 (MQ=255)
                                 tCCGACCTCTCTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGcaca                                                                                                                    <  2:6383936/150‑1 (MQ=255)
                                  ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCacac                                                                                                                   <  2:2009158/150‑1 (MQ=255)
                                      ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCCACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAAt                                                                                                               <  2:371768/150‑1 (MQ=255)
                                          tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                           <  2:5028486/150‑1 (MQ=255)
                                          gaTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                           <  2:3360390/149‑1 (MQ=255)
                                           aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                                                                          >  2:5616967/1‑150 (MQ=255)
                                               gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                                                                      >  2:4920321/1‑150 (MQ=255)
                                                 aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                                                                    >  2:6086334/1‑150 (MQ=255)
                                                  ttCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                                                                   <  1:893475/150‑1 (MQ=255)
                                                   tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                                                                  <  1:6323721/150‑1 (MQ=255)
                                                       aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                                                                              <  2:3949202/150‑1 (MQ=255)
                                                       aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                                                                              <  1:6561708/150‑1 (MQ=255)
                                                        aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                                                                             <  2:295926/150‑1 (MQ=255)
                                                        aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                                                                             >  2:4480305/1‑150 (MQ=255)
                                                         aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                                                                            <  1:288969/150‑1 (MQ=255)
                                                             cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                                                                        <  1:1046523/150‑1 (MQ=255)
                                                             cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                                                                        <  1:603622/150‑1 (MQ=255)
                                                               ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                                                                      >  1:4398402/1‑150 (MQ=255)
                                                                cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                     <  2:4657278/150‑1 (MQ=255)
                                                                cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                     <  1:3401858/150‑1 (MQ=255)
                                                                 tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttttt                                                                                                                                       >  2:1328816/1‑99 (MQ=255)
                                                                 tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                    >  2:5291362/1‑150 (MQ=255)
                                                                   aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGtt                                                                                  <  2:3883173/150‑1 (MQ=255)
                                                                    cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                                 <  2:6472939/150‑1 (MQ=255)
                                                                        aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                                                                             <  1:5200554/150‑1 (MQ=255)
                                                                        aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                                                                             <  1:4513785/150‑1 (MQ=255)
                                                                         aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                                                                            <  2:5149755/150‑1 (MQ=255)
                                                                          aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCaa                                                                           <  2:3093177/150‑1 (MQ=255)
                                                                           cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                                                                          <  1:4756762/150‑1 (MQ=255)
                                                                             aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                                                                        >  1:1360321/1‑150 (MQ=255)
                                                                             aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                                                                        >  1:5408046/1‑150 (MQ=255)
                                                                              tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                                                                       <  1:5716416/150‑1 (MQ=255)
                                                                              tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                                                                       >  2:4539727/1‑150 (MQ=255)
                                                                              tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTACGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                                                                       <  2:2242902/150‑1 (MQ=255)
                                                                                gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                                                                     >  1:4128187/1‑150 (MQ=255)
                                                                                 tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                                                                    <  1:1001104/150‑1 (MQ=255)
                                                                                    ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                                                                 >  2:1821334/1‑150 (MQ=255)
                                                                                         ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca                                                            >  2:5504368/1‑150 (MQ=255)
                                                                                          cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac                                                           >  1:3966181/1‑150 (MQ=255)
                                                                                           aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt                                                                                                                                                      >  1:973199/1‑59 (MQ=255)
                                                                                           aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac                                                           >  2:673016/1‑150 (MQ=255)
                                                                                               gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa                                                      >  2:1705973/1‑150 (MQ=255)
                                                                                               gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa                                                      >  2:586543/1‑150 (MQ=255)
                                                                                                 gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc                                                    <  2:6636902/150‑1 (MQ=255)
                                                                                                  cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt                                                   <  1:2316302/150‑1 (MQ=255)
                                                                                                   aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTc                                                                                                                               <  1:2756388/73‑1 (MQ=255)
                                                                                                   aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg                                                  >  2:1244266/1‑150 (MQ=255)
                                                                                                    ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                                                                             <  2:1163387/90‑1 (MQ=255)
                                                                                                             gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATa                                        <  1:4935898/150‑1 (MQ=255)
                                                                                                             gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATa                                        >  1:1280051/1‑150 (MQ=255)
                                                                                                                tAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGc                                     <  1:1521818/150‑1 (MQ=255)
                                                                                                                 aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGcc                                    <  2:283562/150‑1 (MQ=255)
                                                                                                                   tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGt                                  <  1:2502695/150‑1 (MQ=255)
                                                                                                                    ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGtt                                 >  1:6184006/1‑150 (MQ=255)
                                                                                                                       gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAAt                              >  2:2585706/1‑150 (MQ=255)
                                                                                                                           tttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGct                          <  1:4846198/150‑1 (MQ=255)
                                                                                                                           tttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctt                         <  1:3862882/150‑1 (MQ=255)
                                                                                                                            ttAGTCATACAAAACCTGCCCGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctt                         <  1:6765792/150‑1 (MQ=255)
                                                                                                                               gTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttct                      <  1:6086334/150‑1 (MQ=255)
                                                                                                                               gTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttct                      >  2:1419822/1‑150 (MQ=255)
                                                                                                                                  aTAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttt                   >  2:4316077/1‑150 (MQ=255)
                                                                                                                                   tAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCttttt                  >  2:573146/1‑150 (MQ=255)
                                                                                                                                   tAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCttttt                  >  1:464926/1‑150 (MQ=255)
                                                                                                                                   tAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTGTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCttttt                  >  1:1733341/1‑150 (MQ=255)
                                                                                                                                    aaaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt                 <  2:348236/150‑1 (MQ=255)
                                                                                                                                    aaaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt                 >  1:5492120/1‑150 (MQ=255)
                                                                                                                                     aaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg                >  1:1899515/1‑150 (MQ=255)
                                                                                                                                     aaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGGTAATTGCTTTCTTTTTTg                >  1:203418/1‑150 (MQ=255)
                                                                                                                                       aaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc              >  2:1499326/1‑150 (MQ=255)
                                                                                                                                         aCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgt            >  1:5938884/1‑150 (MQ=255)
                                                                                                                                          ccTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgta           >  1:6485537/1‑150 (MQ=255)
                                                                                                                                            tGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaag         <  2:5805247/150‑1 (MQ=255)
                                                                                                                                                 gCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgtaa    <  1:4272817/150‑1 (MQ=255)
                                                                                                                                                  cGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGcgtaagcgtaag   <  1:2264921/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGA  >  NC_000913/4060116‑4060412

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑