Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,224,242 A→T 7.8% intergenic (+335/+37) ycgI → / ← minE pseudogene/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,224,2420AT 7.8% 356.9 / 16.4 103intergenic (+335/+37)ycgI/minEpseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base A (46/49);  new base T (6/2);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 106 total aligned reads displayed.

AACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACG  >  NC_000913/1224094‑1224387
                                                                                                                                                    |                                                                                                                                                 
aaCCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATaa                                                                                                                                                   >  3:395477/1‑149 (MQ=255)
aaCCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATaa                                                                                                                                                   >  3:844118/1‑149 (MQ=255)
   cAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTAtctc                                                                                                                                                <  4:448766/149‑1 (MQ=255)
   cAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAActc                                                                                                                                                >  4:313681/1‑149 (MQ=255)
      gCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGc                                                                                                                                             >  3:238728/1‑149 (MQ=255)
      gCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGc                                                                                                                                             <  3:434153/149‑1 (MQ=255)
      gCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGc                                                                                                                                             <  3:865828/149‑1 (MQ=255)
       cATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGc                                                                                                                                             >  1:451930/1‑148 (MQ=255)
        aTAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCt                                                                                                                                           >  4:420374/1‑149 (MQ=255)
         tAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCt                                                                                                                                           <  1:349487/148‑1 (MQ=255)
          aaaCATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                                                                                                                                         >  4:521185/1‑149 (MQ=255)
              aTTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGa                                                                                                                                     <  3:313681/149‑1 (MQ=255)
               ttttgttttgTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                                                                                                    >  3:75168/1‑149 (MQ=255)
                tttgttttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGAt                                                                                                                                    >  2:19954/1‑148 (MQ=255)
                tttgttttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATa                                                                                                                                   <  4:825446/149‑1 (MQ=255)
                   gttttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                                                                                <  4:844118/149‑1 (MQ=255)
                   gttttgtGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAAt                                                                                                                                <  3:94405/149‑1 (MQ=255)
                       tgtgCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcgc                                                                                                                            <  4:633079/149‑1 (MQ=255)
                       tgtgCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgcg                                                                                                                             <  2:451930/148‑1 (MQ=255)
                          gCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCtt                                                                                                                          <  1:59570/148‑1 (MQ=255)
                            gTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAGACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTa                                                                                                                       <  4:395477/149‑1 (MQ=255)
                                  aCAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCg                                                                                                                  <  2:344766/148‑1 (MQ=255)
                                      ccTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAATAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAACCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCt                                                                                                              <  1:468280/148‑1 (MQ=255)
                                       cTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCt                                                                                                                                           <  2:319376/118‑1 (MQ=255)
                                       cTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGTTATCTCTGCCt                                                                                                                                           >  1:319376/1‑118 (MQ=255)
                                       cTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCtt                                                                                                             >  1:435784/1‑148 (MQ=255)
                                                 aaaGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGcc                                                                                                                                            >  2:4087/1‑107 (MQ=255)
                                                 aaaGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGcc                                                                                                                                            <  1:4087/107‑1 (MQ=255)
                                                      gggAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTg                                                                                              >  1:846548/1‑148 (MQ=255)
                                                          ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTg                                                                                              >  2:832214/1‑144 (MQ=255)
                                                          ggTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTg                                                                                              <  1:832214/144‑1 (MQ=255)
                                                            tGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCg                                                                                        >  2:842253/1‑148 (MQ=255)
                                                              gAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCAGGCTTATTTCAGCTCTTCTGCTTCCGGTa                                                                                     <  3:538950/149‑1 (MQ=255)
                                                               aaaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                                                                                >  1:14046/1‑137 (MQ=255)
                                                               aaaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCtcttc                                                                                                <  2:14046/137‑1 (MQ=255)
                                                               aaaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTaa                                                                                    <  3:806906/149‑1 (MQ=255)
                                                                aaTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTaa                                                                                    <  1:376851/148‑1 (MQ=255)
                                                                 aTAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAg                                                                                   <  1:784129/148‑1 (MQ=255)
                                                                             agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTcc                                                                                         <  1:326658/130‑1 (MQ=255)
                                                                             agTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTcc                                                                                         >  2:326658/1‑130 (MQ=255)
                                                                              gTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTc                                                                     <  4:75168/149‑1 (MQ=255)
                                                                              gTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTc                                                                     >  3:782380/1‑149 (MQ=255)
                                                                               tATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCg                                                                                        <  3:271872/129‑1 (MQ=255)
                                                                               tATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCg                                                                                        >  4:271872/1‑129 (MQ=255)
                                                                               tATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCa                                                                    >  3:442912/1‑149 (MQ=255)
                                                                                    gCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGaa                                                                >  1:38978/1‑148 (MQ=255)
                                                                                           tACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa                                                                 >  1:610008/1‑140 (MQ=255)
                                                                                           tACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGa                                                                 <  2:610008/140‑1 (MQ=255)
                                                                                             cTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAatat                                                       <  2:846548/148‑1 (MQ=255)
                                                                                              tCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATc                                                      >  1:257486/1‑148 (MQ=255)
                                                                                               cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGt                                                                            <  1:231286/125‑1 (MQ=255)
                                                                                               cgccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGt                                                                            >  2:231286/1‑125 (MQ=255)
                                                                                                gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCa                                                                               >  1:376936/1‑121 (MQ=255)
                                                                                                gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCa                                                                               <  2:376936/121‑1 (MQ=255)
                                                                                                gccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGcc                                                   >  3:573088/1‑149 (MQ=255)
                                                                                                 ccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGcc                                                   >  1:743666/1‑148 (MQ=255)
                                                                                                 ccgcTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGcc                                                   <  1:648416/148‑1 (MQ=255)
                                                                                                       tactataaCGTGGATGAGAGGGATGAATAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCtttt                                            <  4:782380/143‑1 (MQ=255)
                                                                                                       tACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCttt                                             <  1:815818/148‑1 (MQ=255)
                                                                                                        aCCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCttttt                                           <  4:238728/149‑1 (MQ=255)
                                                                                                          cacaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTg                                          >  1:332153/1‑148 (MQ=255)
                                                                                                          cacaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGc                                         >  4:816133/1‑149 (MQ=255)
                                                                                                           acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCtttt                                            <  3:269509/145‑1 (MQ=255)
                                                                                                           acaACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCtttt                                            >  4:269509/1‑145 (MQ=255)
                                                                                                            caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGt                                                                                 >  2:671373/1‑107 (MQ=255)
                                                                                                            caACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGt                                                                                 <  1:671373/107‑1 (MQ=255)
                                                                                                              aCGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCaa                                     <  3:400008/149‑1 (MQ=255)
                                                                                                                     tGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTa                               <  1:266859/148‑1 (MQ=255)
                                                                                                                       agagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACg                             <  2:743666/148‑1 (MQ=255)
                                                                                                                        gagGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATGGCCATCTTTTTGCTCAAGCTGTACGGt                           <  3:420374/149‑1 (MQ=255)
                                                                                                                         agGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAAt                                                               >  1:425822/1‑112 (MQ=255)
                                                                                                                         agGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAAt                                                               <  2:425822/112‑1 (MQ=255)
                                                                                                                           ggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                                                                      <  3:69846/71‑1 (MQ=255)
                                                                                                                           ggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAg                                                                                                      >  4:69846/1‑71 (MQ=255)
                                                                                                                           ggATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTAc                        <  3:275146/149‑1 (MQ=255)
                                                                                                                            gATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATc                                                <  1:720065/124‑1 (MQ=255)
                                                                                                                            gATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATc                                                >  2:720065/1‑124 (MQ=255)
                                                                                                                            gATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTAcc                       <  4:573088/149‑1 (MQ=255)
                                                                                                                             aTGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTAcc                       >  2:911604/1‑148 (MQ=255)
                                                                                                                              tGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCAt                     >  3:673389/1‑149 (MQ=255)
                                                                                                                                aaaaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCAtct                   >  4:337056/1‑149 (MQ=255)
                                                                                                                                  aaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCAtctc                  <  1:331366/148‑1 (MQ=255)
                                                                                                                                  aaaCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCAtctc                  >  1:208508/1‑148 (MQ=255)
                                                                                                                                    aCTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAgg               >  4:256111/1‑149 (MQ=255)
                                                                                                                                      tCAAGGCAGAGTTATCTCTGCCTTGAGGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGa              >  2:406200/1‑148 (MQ=255)
                                                                                                                                      tCAAGGCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGa              >  2:48581/1‑148 (MQ=255)
                                                                                                                                      tCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGa              <  2:306679/148‑1 (MQ=255)
                                                                                                                                       cAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCa                                                                    >  2:681423/1‑93 (MQ=255)
                                                                                                                                       cAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCa                                                                    <  1:681423/93‑1 (MQ=255)
                                                                                                                                        aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGtt                          >  4:423418/1‑134 (MQ=255)
                                                                                                                                        aaGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCa           >  4:742177/1‑149 (MQ=255)
                                                                                                                                         aGGCACAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCa           <  1:572851/148‑1 (MQ=255)
                                                                                                                                          ggCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATaaa          >  1:333837/1‑148 (MQ=255)
                                                                                                                                          ggCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCaa          <  1:694309/148‑1 (MQ=255)
                                                                                                                                          ggCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAAt         <  3:816133/149‑1 (MQ=255)
                                                                                                                                           gCAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAAt         >  1:878336/1‑148 (MQ=255)
                                                                                                                                           gCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAAt         >  2:40150/1‑148 (MQ=255)
                                                                                                                                            cAGAGTTATCTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAAtt        >  2:222487/1‑148 (MQ=255)
                                                                                                                                                 aTAACTCTTCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCGTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTAc   >  1:109460/1‑148 (MQ=255)
                                                                                                                                                  tAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACg  >  1:753143/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AACCAGGCATAAACATTTTGTTTTGTGCGTGGGAACAGCCTTAAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTACAGCGGGCTTATTTCAGCTCTTCTGCTTCCGGTAAGGTCACGTTCAGCTCAAGAATAGAAATATCGCCATCTTTTTGCTCAAGCTGTACGGTTACCATCTCAGGATCAATTTGTACG  >  NC_000913/1224094‑1224387

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: