Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,579,152 A→C 5.2% G62G (GGT→GGG)  mcrB ← 5‑methylcytosine‑specific restriction enzyme subunit McrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,579,1520AC 5.2% 414.0 / 11.9 136G62G (GGT→GGG) mcrB5‑methylcytosine‑specific restriction enzyme subunit McrB
Reads supporting (aligned to +/- strand):  ref base A (68/61);  new base C (3/4);  total (71/65)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

Only 100 of 142 total aligned reads displayed.

ACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAACGTTGTTGCTGTGCTTGCT  >  NC_000913/4579003‑4579300
                                                                                                                                                     |                                                                                                                                                    
acacCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGAATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGAtatata                                                                                                                                                      <  1:649121/150‑1 (MQ=255)
 cacCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCTTGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATAc                                                                                                                                                     <  2:19373/150‑1 (MQ=255)
       aaGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGTATATATCCCGTTat                                                                                                                                               <  1:803427/150‑2 (MQ=255)
       aaGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGTATATATCCCGTTat                                                                                                                                               <  1:704150/150‑2 (MQ=255)
        aGTTGCCTGAAAATACTCTGCGATTTTTTTATGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTAAACCATAATGCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGa                                                                                                                                              <  2:1626678/150‑1 (MQ=255)
        aGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGa                                                                                                                                              <  1:1534318/150‑1 (MQ=255)
           tGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGTATATATCCCGTTattctc                                                                                                                                           <  2:143958/150‑6 (MQ=255)
           tGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGTATATATCCCGTTattctc                                                                                                                                           <  1:1632105/150‑6 (MQ=255)
           tGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGc                                                                                                                                           <  2:1005078/150‑1 (MQ=255)
           tGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCAATGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGc                                                                                                                                           >  2:1611510/1‑150 (MQ=255)
            gCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCt                                                                                                                                          >  1:1104896/1‑150 (MQ=255)
               tGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTcc                                                                                                                                       >  1:911298/1‑150 (MQ=255)
                gAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCt                                                                                                                                      >  1:1069918/1‑150 (MQ=255)
                gAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCt                                                                                                                                      <  1:950863/150‑1 (MQ=255)
                   aaTACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGAc                                                                                                                                   >  2:237847/1‑150 (MQ=255)
                          tGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTcc                                                                                                                            >  1:1451249/1‑150 (MQ=255)
                          tGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTcc                                                                                                                            >  2:752972/1‑150 (MQ=255)
                           gCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCa                                                                                                                           >  2:1074626/1‑150 (MQ=255)
                            cGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCaa                                                                                                                          >  1:466540/1‑150 (MQ=255)
                             gATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAg                                                                                                                         >  1:1524923/1‑150 (MQ=255)
                                 gTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAAt                                                                                                                     <  1:275154/150‑1 (MQ=255)
                                   tttAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGc                                                                                                                   >  1:1504530/1‑150 (MQ=255)
                                   tttAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGc                                                                                                                   <  1:414495/150‑1 (MQ=255)
                                       ggTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAAc                                                                                                               >  1:855187/1‑150 (MQ=255)
                                        gTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAAcc                                                                                                              <  1:124462/150‑1 (MQ=255)
                                            gTCTGAAGAGAACTGGCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAggg                                                                                                          <  2:553842/150‑1 (MQ=255)
                                              cTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAACCCAGGGaa                                                                                                        >  1:171481/1‑150 (MQ=255)
                                                    agaACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACg                                                                                                  <  2:900749/150‑1 (MQ=255)
                                                    agaACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACg                                                                                                  >  1:710413/1‑150 (MQ=255)
                                                    agaACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACg                                                                                                  <  1:1304889/150‑1 (MQ=255)
                                                        cTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTaaa                                                                                              >  2:110481/1‑150 (MQ=255)
                                                        cTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTaaa                                                                                              <  2:1469997/150‑1 (MQ=255)
                                                        cTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTaaa                                                                                              >  1:1436036/1‑150 (MQ=255)
                                                         tGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTaaaa                                                                                             <  2:1451249/150‑1 (MQ=255)
                                                            cATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTa                                                                                          <  2:185992/150‑1 (MQ=255)
                                                            cATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTa                                                                                          >  1:1038316/1‑150 (MQ=255)
                                                             aTTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTAc                                                                                         <  2:1379008/150‑1 (MQ=255)
                                                               tGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCa                                                                                       <  1:1304897/150‑1 (MQ=255)
                                                                     tgtgGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCg                                                                                 <  2:1385922/150‑1 (MQ=255)
                                                                     tgtgGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCg                                                                                 <  2:1455952/150‑1 (MQ=255)
                                                                      gtgGTTGATTCGTGTCGCTAATACCATAAGCCAAAACCAACTCCTCAAAATCTATATAATAGAGAATAACGGGATATATACCGTTAGAAGATACCTGACGTTCTCCAAGAAATGCAAAACAGGGAAAAGACGTAAAATTACCATAAACGa                                                                                >  2:1289898/1‑150 (MQ=255)
                                                                       tgGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGaa                                                                               <  1:1603111/150‑1 (MQ=255)
                                                                         gTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGATGCTTCCTGACCTTCTCCAATAAATTCAAACCAGTGAATAGACTTAAAATTAACATAACCGaaag                                                                             <  2:1227217/150‑2 (MQ=255)
                                                                              ttCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAAt                                                                        <  1:1158440/150‑1 (MQ=255)
                                                                              ttCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAAt                                                                        <  2:1197453/150‑1 (MQ=255)
                                                                                cGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAAttt                                                                      <  2:466540/150‑1 (MQ=255)
                                                                                 gtgtCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAAATCAAtttt                                                                     >  1:703057/1‑150 (MQ=255)
                                                                                  tgtCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTa                                                                    <  1:1213321/150‑1 (MQ=255)
                                                                                   gtCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTAc                                                                   <  1:367950/150‑1 (MQ=255)
                                                                                     cGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTc                                                                 <  2:344453/150‑1 (MQ=255)
                                                                                      gCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCg                                                                >  2:512267/1‑150 (MQ=255)
                                                                                        ttATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCa                                                              <  2:1512203/150‑1 (MQ=255)
                                                                                             ccATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCCAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGttt                                                         >  2:589613/1‑149 (MQ=255)
                                                                                             ccATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTa                                                         >  1:665237/1‑150 (MQ=255)
                                                                                             ccATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACAGTTAGAAGCTTCCTGACCTTCTCCAACAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAAATTTACTCGCAGGTTa                                                         >  2:230683/1‑150 (MQ=255)
                                                                                              cATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACGATAACCGAAACTCAATTTTACTCGCAGGTTAc                                                        >  2:973864/1‑150 (MQ=255)
                                                                                               aTAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACg                                                       >  2:1629895/1‑150 (MQ=255)
                                                                                                tAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACgg                                                      <  2:742677/150‑1 (MQ=255)
                                                                                                    ccAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAg                                                  >  1:505367/1‑150 (MQ=255)
                                                                                                     cAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGa                                                 <  1:752972/150‑1 (MQ=255)
                                                                                                      aaaaCCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGAc                                                >  2:524573/1‑150 (MQ=255)
                                                                                                       aaaCCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGa                                                 <  2:1419337/148‑1 (MQ=255)
                                                                                                       aaaCCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGa                                                 >  1:1419337/1‑148 (MQ=255)
                                                                                                        aaCCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGa                                                                                                                        <  1:475345/76‑1 (MQ=255)
                                                                                                        aaCCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGa                                                                                                                        >  2:475345/1‑76 (MQ=255)
                                                                                                          ccAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTg                                            >  2:571803/1‑150 (MQ=255)
                                                                                                          ccAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTg                                            <  1:1611510/150‑1 (MQ=255)
                                                                                                            aaCTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATTGACGTAATATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGa                                          <  2:1069918/150‑1 (MQ=255)
                                                                                                             aCTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGAt                                         <  2:1038316/150‑1 (MQ=255)
                                                                                                             aCTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGAt                                         <  2:1375150/150‑1 (MQ=255)
                                                                                                             aCTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGAt                                         <  1:571803/150‑1 (MQ=255)
                                                                                                              cTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATa                                        >  1:59812/1‑150 (MQ=255)
                                                                                                                catcaAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAAt                                      >  2:278152/1‑150 (MQ=255)
                                                                                                                catcaAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAAt                                      >  2:1436417/1‑150 (MQ=255)
                                                                                                                 atcaAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTgg                                           >  1:542794/1‑144 (MQ=255)
                                                                                                                 atcaAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATc                                     >  2:1590863/1‑150 (MQ=255)
                                                                                                                 atcaAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATc                                     <  2:1529939/150‑1 (MQ=255)
                                                                                                                 atcaAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTAAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATc                                     >  2:830487/1‑150 (MQ=255)
                                                                                                                   caAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCtt                                   <  2:406374/150‑1 (MQ=255)
                                                                                                                    aaaaTCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCttt                                  >  2:1021151/1‑150 (MQ=255)
                                                                                                                    aaaaTCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATAttt                                  >  2:1180744/1‑150 (MQ=255)
                                                                                                                       aTCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGt                               <  2:855187/150‑1 (MQ=255)
                                                                                                                       aTCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGt                               >  2:69719/1‑150 (MQ=255)
                                                                                                                         cTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTgg                             >  2:962980/1‑150 (MQ=255)
                                                                                                                         cTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTgg                             >  2:1240041/1‑150 (MQ=255)
                                                                                                                           ttATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGAt                           >  2:1063593/1‑150 (MQ=255)
                                                                                                                              taataGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGc                        >  2:808494/1‑150 (MQ=255)
                                                                                                                                ataGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGa                      >  1:326057/1‑150 (MQ=255)
                                                                                                                                ataGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGa                      >  2:691051/1‑150 (MQ=255)
                                                                                                                                  agagAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAAc                    >  2:832600/1‑150 (MQ=255)
                                                                                                                                   gagaATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAACg                   >  1:1215824/1‑150 (MQ=255)
                                                                                                                                     ttctAACGGTATATATCCCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAACgtt                 >  1:1344766/4‑150 (MQ=255)
                                                                                                                                     gaATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAACgtt                 <  2:1072096/150‑1 (MQ=255)
                                                                                                                                       aTAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAAGgttgt               >  1:562329/1‑150 (MQ=255)
                                                                                                                                         aaCGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAACGttgttg             <  1:1073127/150‑1 (MQ=255)
                                                                                                                                          aCGGGATATATACCGTTAGAAGCTTCCTGACCTACTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAACGTTGTTGc            >  1:415972/1‑150 (MQ=255)
                                                                                                                                             ggATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAACGTTGTTGCtgt         >  2:204633/1‑150 (MQ=255)
                                                                                                                                                  tataCCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAACGTTGTTGCTGtgcttg    >  2:105037/1‑150 (MQ=255)
                                                                                                                                                   ataCCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAACGTTGTTGCTGtgcttgc   >  1:1312959/1‑150 (MQ=255)
                                                                                                                                                    taCCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAACGTTGTTGCTGtgcttgct  >  2:1512394/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAATCTTTATAATAGAGAATAACGGGATATATACCGTTAGAAGCTTCCTGACCTTCTCCAAGAAATGCAAACCAGGGAATAGACGTAAAATTACCATAACCGAAACTCAATTTTACTCGCAGGTTACGGTAAGACGTTGGATAATCTTTAGTGGATTGCGAACGTTGTTGCTGTGCTTGCT  >  NC_000913/4579003‑4579300

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