Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,016 T→C 6.9% intergenic (+32/‑28) secA → / → mutT preprotein translocase subunit, ATPase/dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,0160TC 6.9% 419.5 / 16.4 146intergenic (+32/‑28)secA/mutTpreprotein translocase subunit, ATPase/dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (61/74);  new base C (7/3);  total (68/78)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

Only 100 of 147 total aligned reads displayed.

ACTCTGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCG  >  NC_000913/110866‑111155
                                                                                                                                                      |                                                                                                                                           
acTCTGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAGGCGCAGGAtt                                                                                                                                             >  1:1060835/1‑150 (MQ=255)
      cAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgcg                                                                                                                                        <  1:164098/150‑1 (MQ=255)
        gCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTGTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCGAACTGTGGAAGTAAAAGGCGCAGGATTCTGCGcc                                                                                                                                      >  1:1073971/1‑150 (MQ=255)
          cGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCtt                                                                                                                                    >  1:356279/1‑150 (MQ=255)
          cGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCtt                                                                                                                                    >  2:726707/1‑150 (MQ=255)
           gCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTATGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttt                                                                                                                                   <  1:1510567/150‑1 (MQ=255)
           gCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCttt                                                                                                                                   <  2:1320218/150‑1 (MQ=255)
             aGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCttttt                                                                                                                                 <  2:74720/150‑1 (MQ=255)
              gCTGCACTGGCGGCGCAACCCGGAGAGCGCAAAGTAGGACGTAACGATCCGTGCCGGTGGGGTTCTGGTAAAAAATACAAGCTGTGCCATGGCCGCGTGCAATAAAAGCTTGGTGTTGAAGTAAAAGGCGCAGGATTCTGCGCGtgtttt                                                                                                                                >  2:670094/1‑150 (MQ=255)
               cTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATGCAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTta                                                                                                                               >  2:302356/1‑150 (MQ=255)
               cTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttttact                                                                                                                               <  1:459193/150‑4 (MQ=255)
                tGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtat                                                                                                                              <  2:420825/150‑1 (MQ=255)
                tGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCATTGCCCTTCCGGATCTGGTAAAAAATACAAGCAGTGCCATGCGCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtat                                                                                                                              >  2:326845/1‑150 (MQ=255)
                 gCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttttacttc                                                                                                                             <  2:1063896/150‑6 (MQ=255)
                  cACTGGCGGCGGAAACCGGACAGCGCAAAGTAGGACGCTACGATCCTTGCACGTGAGGTTCTGCTAAAAACTACCAGCAGGGCCATGGCCGCCTGCAATAAAAGCTAACTGCTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAg                                                                                                                            <  1:272876/150‑1 (MQ=255)
                   aCTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGAGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAgg                                                                                                                           >  1:1203365/1‑150 (MQ=255)
                      ggcggcGCAAACCGGATAGCGCAAAGTATGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCAAGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGttt                                                                                                                        <  2:862118/150‑1 (MQ=255)
                      ggcggcGCAAACCGGAGATCGCAAAGTAGGACGTACAGATCCTTGCCCGTGAGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGttt                                                                                                                        <  2:820646/150‑1 (MQ=255)
                      ggcggcGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGttt                                                                                                                        <  2:909313/150‑1 (MQ=255)
                      ggcggcGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTa                                                                                                                       >  1:386046/1‑150 (MQ=255)
                       gcggcgCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTa                                                                                                                       <  1:807191/150‑1 (MQ=255)
                         ggcgCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAg                                                                                                                     >  1:1173015/1‑150 (MQ=255)
                          gcgcAAACCGGAGAGCGCAAAGTAGGACGTAACTATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAACGCTAAGTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGa                                                                                                                    >  1:552419/1‑150 (MQ=255)
                           cgcAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAc                                                                                                                   >  1:1342257/1‑150 (MQ=255)
                               aaCCGGAGAGCGCAAAGTAGGACGTAACGCTCCTTGCCCGTGCGGTGCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATg                                                                                                               <  2:1091281/150‑1 (MQ=255)
                               aaCCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATg                                                                                                               <  2:792275/150‑1 (MQ=255)
                               aaCCGGAGAGCGCAAAGTAGGACGAAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATg                                                                                                               <  2:836107/150‑1 (MQ=255)
                                    gagaGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaaaaa                                                                                                          <  1:587442/150‑1 (MQ=255)
                                     agagCGCAAAGTAGGCCGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAACGCTAACGGTTGAAGTAAAAGGCGCACGATTCTGCGCGTTTTTTATAGGTTTAAGACAATGAAAAAg                                                                                                         >  2:89277/1‑150 (MQ=255)
                                         cgcAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGc                                                                                                     <  2:356279/150‑1 (MQ=255)
                                          gcAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCa                                                                                                    <  2:797369/150‑1 (MQ=255)
                                               gtagtaCGTAACGATCCTTGCCAGTGCGGATCTTGTAAAAAATACAAGCAGTTCCATGGCATACTGAAATAAAATCTAACTGTTGAAGTAAAAGTCGCAGTATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTg                                                                                               <  1:1239770/150‑1 (MQ=255)
                                               gTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTg                                                                                               >  2:1356218/1‑150 (MQ=255)
                                                 aGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                                                             >  2:489781/1‑150 (MQ=255)
                                                  ggACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTTCAATAAATGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCgg                                                                                            <  1:426273/150‑1 (MQ=255)
                                                  ggACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCgg                                                                                            <  1:1491709/150‑1 (MQ=255)
                                                  ggACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCgg                                                                                            <  2:202725/150‑1 (MQ=255)
                                                  ggACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGATGCAAATTGCgg                                                                                            >  1:169671/1‑150 (MQ=255)
                                                     cGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGTAGCAGGATTCTGCGCCTTTTTTATAGGGTTAAGAAAATGAAATAGCTTCAAATTGCggtag                                                                                         <  1:1542259/150‑1 (MQ=255)
                                                     cGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtag                                                                                         >  1:812301/1‑150 (MQ=255)
                                                     cGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGACTTTTTTAAAGGTTTAAGACAATGAAAAAGATGCAATTAGCggaat                                                                                         >  2:1344022/1‑147 (MQ=255)
                                                      gTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtagg                                                                                        <  1:190925/150‑1 (MQ=255)
                                                       tAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtaggt                                                                                       <  1:1091611/150‑1 (MQ=255)
                                                       tAACGATCCTTGCCCGAGCTGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtaggt                                                                                       <  2:770808/150‑1 (MQ=255)
                                                       tAACGAGCCTTGCCCGTTCGTGTCTGGTAAAAAATACAATCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGGTGACGTAAAAGTCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtaggt                                                                                       <  2:1216270/150‑1 (MQ=255)
                                                         aCGAGCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGtat                                                                                     <  1:402824/150‑1 (MQ=255)
                                                           gATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGtatta                                                                                   >  2:1255876/1‑150 (MQ=255)
                                                             gccatGCCCGTGCGGTTCTGGTAAAAAAAACAAGCAGTGCCATGGCCGCCAGCAATAATAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTAATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTattatt                                                                                 <  2:282127/146‑1 (MQ=255)
                                                              ccTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTc                                                                                >  2:1610413/1‑150 (MQ=255)
                                                              ccTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGATAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTc                                                                                <  2:1479463/150‑1 (MQ=255)
                                                                         cGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAAc                                                                     <  1:158178/150‑1 (MQ=255)
                                                                         cGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAAc                                                                     <  2:1527463/150‑1 (MQ=255)
                                                                            ttCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAAt                                                                  >  2:1036700/1‑150 (MQ=255)
                                                                            ttCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAAt                                                                  >  2:1017163/1‑150 (MQ=255)
                                                                              cTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGa                                                                <  1:1544132/150‑1 (MQ=255)
                                                                              cTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGCCGCAGGATTCTGCGCCGTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGa                                                                >  1:1089697/1‑150 (MQ=255)
                                                                                ggTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGaaa                                                              >  1:1574076/1‑150 (MQ=255)
                                                                                ggTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGaaa                                                              >  1:693848/1‑150 (MQ=255)
                                                                                   aaaaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCt                                                           <  2:161563/150‑1 (MQ=255)
                                                                                   aaaaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGATATTATTCGCAACGTAACAATGAAATCtt                                                          >  2:1518664/1‑150 (MQ=255)
                                                                                    aaaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCtt                                                          <  2:976113/150‑1 (MQ=255)
                                                                                     aaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                         <  2:975099/150‑1 (MQ=255)
                                                                                      aaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTta                                                        <  1:1043547/150‑1 (MQ=255)
                                                                                      aaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTta                                                        <  1:866839/150‑1 (MQ=255)
                                                                                       attACAAGCAGTGCCATGGCCGACTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtat                                                       <  1:352431/148‑1 (MQ=255)
                                                                                        atgcAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtata                                                      <  1:302356/147‑1 (MQ=255)
                                                                                        aTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtata                                                      <  1:699278/150‑1 (MQ=255)
                                                                                               cAGTGCAATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTATTATAGGTTTAAGACTATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATATTTATAACGCGTc                                               >  2:72499/1‑150 (MQ=255)
                                                                                                aGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTcg                                              <  1:1668393/150‑1 (MQ=255)
                                                                                                aGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTcg                                              <  1:1674863/150‑1 (MQ=255)
                                                                                                   gCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCgcgc                                           >  2:212700/1‑150 (MQ=255)
                                                                                                      aTGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTATTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACTCGTCGCgcagc                                        >  2:1407392/1‑150 (MQ=255)
                                                                                                       tGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCgcagca                                       <  2:1485554/150‑1 (MQ=255)
                                                                                                         gccgccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGa                                     >  1:1212559/1‑150 (MQ=255)
                                                                                                              cTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCATTATAACGCGTCGCGCAGCAGATgcgc                                >  2:1388306/1‑150 (MQ=255)
                                                                                                               tGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGAAGCGca                               >  1:669462/1‑150 (MQ=255)
                                                                                                               tGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTATTTTATAGGGTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAAGGACAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGca                               >  2:660952/1‑150 (MQ=255)
                                                                                                                gCAATAAAAGCTAACACTTGAAGTAAAAGGCGCAGGAATCTGCGCCTTTTTAATAGGTTTAAGACAATGAAAAAGCTGCAAATTCCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACTCGACGAGCAGCAGAACCGcac                              >  2:1734140/1‑150 (MQ=255)
                                                                                                                 cAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGcaca                             <  1:673612/150‑1 (MQ=255)
                                                                                                                     aaaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGc                         <  1:373893/150‑1 (MQ=255)
                                                                                                                           tAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACACGTCGCGCAGCAGATGCGCACATGGCtaataa                   >  2:848134/1‑150 (MQ=255)
                                                                                                                              cTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACt                >  2:1678578/1‑150 (MQ=255)
                                                                                                                              cTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTAGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACt                >  2:1672059/1‑150 (MQ=255)
                                                                                                                               tGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTg               <  1:489781/150‑1 (MQ=255)
                                                                                                                                gTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTgg              <  2:466203/150‑1 (MQ=255)
                                                                                                                                 cctataaaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa             >  1:23212/8‑150 (MQ=255)
                                                                                                                                  tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTATGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAg            <  2:713461/150‑1 (MQ=255)
                                                                                                                                  tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAg            <  2:707351/150‑1 (MQ=255)
                                                                                                                                  tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAg            <  1:1167413/150‑1 (MQ=255)
                                                                                                                                  tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAg            <  2:1075872/150‑1 (MQ=255)
                                                                                                                                  tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAg            <  1:477467/150‑1 (MQ=255)
                                                                                                                                   gAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGt           >  1:374610/1‑150 (MQ=255)
                                                                                                                                   gAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGt           >  2:710303/1‑150 (MQ=255)
                                                                                                                                     aGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGttt         <  2:1404957/150‑1 (MQ=255)
                                                                                                                                        aaaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTccc      >  2:605614/1‑150 (MQ=255)
                                                                                                                                         aaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcg     >  1:453839/1‑150 (MQ=255)
                                                                                                                                         aaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcg     >  1:865171/1‑150 (MQ=255)
                                                                                                                                         aaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcg     >  2:626290/1‑150 (MQ=255)
                                                                                                                                           aGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCTTCTCGCAGCAGATGCGCACATGTCGAATAAACTGGAGTTTCAcggc   >  1:224674/1‑150 (MQ=255)
                                                                                                                                            ggCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcggcg  <  2:416963/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                           
ACTCTGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCG  >  NC_000913/110866‑111155

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