Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,778,364 C→T 6.2% intergenic (+17/+26) ydiO → / ← ydiP putative acyl‑CoA dehydrogenase/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,778,3640CT 6.2% 254.1 / 4.1 81intergenic (+17/+26)ydiO/ydiPputative acyl‑CoA dehydrogenase/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (39/36);  new base T (2/3);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GGCTATACCGATGAGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTT  >  NC_000913/1778216‑1778509
                                                                                                                                                    |                                                                                                                                                 
ggCTAAACCGATGAGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTAAGTAGCAGCTCGGCAGATCCTGAAAGACTAACAGAACAAATAATCAGCACGCGGCGCAGct                                                                                                                                                  >  2:41/1‑150 (MQ=255)
    ataCCGATGAGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGcttctt                                                                                                                                              >  1:1000623/1‑150 (MQ=255)
      aCCGATGAGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCCGAACAAATAATCTGCAGGCGGCGCAGCTTCTTaa                                                                                                                                            >  1:1154570/1‑150 (MQ=255)
      aCCGATGAGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTaa                                                                                                                                            >  1:1005153/1‑150 (MQ=255)
      aCCGATGAGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTaa                                                                                                                                            >  2:133781/1‑150 (MQ=255)
      aCCGATGAGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTaa                                                                                                                                            <  2:1077989/150‑1 (MQ=255)
      aCCGATGAGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTAAAGACGAAATTATGATTTACGTCGCAGGTCGGCAGATCCTGAAAGACTATCAGAACATATTTTCTGACGGCGGCGCAGctgcgtcc                                                                                                                                            >  2:1083787/1‑144 (MQ=255)
        cGATGAGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACa                                                                                                                                          <  2:1296769/150‑1 (MQ=255)
           tGAGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAAc                                                                                                                                       <  1:688579/150‑1 (MQ=255)
             aGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTg                                                                                                                                     <  2:1269237/150‑1 (MQ=255)
             aGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTg                                                                                                                                     >  2:778760/1‑150 (MQ=255)
             aGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTg                                                                                                                                     <  2:857253/150‑1 (MQ=255)
              ggCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTgc                                                                                                                                    >  1:153635/1‑150 (MQ=255)
                 gcgcgTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGcgcc                                                                                                                                 <  2:1511177/150‑1 (MQ=255)
                 gcgcgTCTCCCGCTTCTGGCGGGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGcgcc                                                                                                                                 <  2:777454/150‑1 (MQ=255)
                  cgcgTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGcgccg                                                                                                                                >  1:1671436/1‑150 (MQ=255)
                   gcgTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGcgccgc                                                                                                                               <  1:1422854/150‑1 (MQ=255)
                    cgTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCgccgcc                                                                                                                              <  1:159765/150‑1 (MQ=255)
                      gctcCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcctg                                                                                                                            <  2:1716918/149‑3 (MQ=255)
                       ctcCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGa                                                                                                                           >  2:768219/1‑150 (MQ=255)
                         cccGCTTCTGGCGGGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGAtt                                                                                                                         <  2:828683/150‑1 (MQ=255)
                          ccGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGAttt                                                                                                                        >  2:1219582/1‑150 (MQ=255)
                           cGCTTCTGGAGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTa                                                                                                                       <  1:919752/150‑1 (MQ=255)
                                  ggCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaaca                                                                                                                >  1:599295/1‑150 (MQ=255)
                                    cGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGATTCTTAACAAACTGCGACACCAGATATATCCAACAAg                                                                                                              >  1:233351/1‑150 (MQ=255)
                                    cGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATGTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAAAAATCTGCAGGCGGCGCAGCGTCTTAACCAACTGCCCCGCCAGAATTATCcaacaac                                                                                                              >  1:1014727/1‑149 (MQ=255)
                                      tGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGAc                                                                                                            <  2:866292/150‑1 (MQ=255)
                                       gATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACt                                                                                                           >  1:197438/1‑150 (MQ=255)
                                           tCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTAcc                                                                                                       <  1:1316212/150‑1 (MQ=255)
                                           tCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTAcc                                                                                                       <  1:778760/150‑1 (MQ=255)
                                            ccGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCg                                                                                                      >  1:1413271/1‑150 (MQ=255)
                                              gTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGt                                                                                                    <  1:768219/150‑1 (MQ=255)
                                              gTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGt                                                                                                    <  2:1295904/150‑1 (MQ=255)
                                                   gAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAgg                                                                                               >  1:629198/1‑150 (MQ=255)
                                                    taCGTATCGGCGGCGGTACAGAGGAAATTATTATTTACGTAGCAGGTCGGCAGATCCAGAAAGACTATCAGAACAAATTATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCTGTTTAATCCAACAAGAGTTACCGGTTGAGGa                                                                                              <  2:302206/149‑1 (MQ=255)
                                                       gTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAAt                                                                                           <  1:1672400/150‑1 (MQ=255)
                                                         aTCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTc                                                                                         >  2:484482/1‑150 (MQ=255)
                                                          tCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCa                                                                                        <  1:133781/150‑1 (MQ=255)
                                                            ggcggcggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGa                                                                                      <  1:507503/150‑1 (MQ=255)
                                                              cggcggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATa                                                                                    >  1:922019/1‑150 (MQ=255)
                                                                   gTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGAtt                                                                               <  2:353821/150‑1 (MQ=255)
                                                                   gTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGAtt                                                                               >  1:1022713/1‑150 (MQ=255)
                                                                             aaTTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCg                                                                     >  1:712877/1‑150 (MQ=255)
                                                                               ttATGATTTACGTAGCAGGGCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCa                                                                   <  2:153635/150‑1 (MQ=255)
                                                                                 aTGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAATACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCaaa                                                                 >  2:1087420/1‑150 (MQ=255)
                                                                                       tACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtt                                                           >  2:755428/1‑150 (MQ=255)
                                                                                         cGTCGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCGAGATTTATCCAACATGACTTACCGGTTGTGGTAATTCAGATAGGATTTCTTGTTCCGCAAAACTGtttt                                                         <  1:41/150‑1 (MQ=255)
                                                                                         cGTAGCAGGTCGGCAGTTCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtttt                                                         >  2:1116398/1‑150 (MQ=255)
                                                                                         cGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAAATCAGATAGGATTTCTTGTTCCGCAAAACAGtttt                                                         >  2:1622852/1‑150 (MQ=255)
                                                                                          gTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttttt                                                        >  1:1547879/1‑150 (MQ=255)
                                                                                          gTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttttt                                                        <  1:917785/150‑1 (MQ=255)
                                                                                           tAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTa                                                       <  1:1219582/150‑1 (MQ=255)
                                                                                           tAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTa                                                       <  2:1671436/150‑1 (MQ=255)
                                                                                            aGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTaa                                                      >  2:15802/1‑150 (MQ=255)
                                                                                            aGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTAGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTATaa                                                      >  2:9826/1‑150 (MQ=255)
                                                                                             gCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTaaa                                                     <  2:1109571/150‑1 (MQ=255)
                                                                                             gCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTaaa                                                     <  2:1550324/150‑1 (MQ=255)
                                                                                              cAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAAt                                                    >  2:431359/1‑150 (MQ=255)
                                                                                              cAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAAt                                                    <  2:197438/150‑1 (MQ=255)
                                                                                                gtgcGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCTCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTg                                                  <  1:947900/147‑1 (MQ=255)
                                                                                                  tCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTc                                                >  1:936791/1‑150 (MQ=255)
                                                                                                           ccTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAAt                                       <  2:1084428/150‑1 (MQ=255)
                                                                                                            cTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATc                                      <  1:1239913/150‑1 (MQ=255)
                                                                                                             tGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCg                                     >  2:1187978/1‑150 (MQ=255)
                                                                                                              gAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGTTAATCGc                                    >  1:818380/1‑150 (MQ=255)
                                                                                                                    cTATCAGAACAAATACTCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCg                              <  1:1629744/150‑1 (MQ=255)
                                                                                                                       tCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTACCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTTGTCGaac                           >  2:262925/1‑150 (MQ=255)
                                                                                                                             cAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCATGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGa                     >  1:847143/1‑150 (MQ=255)
                                                                                                                               aataatCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACa                   >  1:279169/1‑150 (MQ=255)
                                                                                                                               aataatCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACa                   >  1:1503229/1‑150 (MQ=255)
                                                                                                                               aataatCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGGTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACa                   >  1:837916/1‑150 (MQ=255)
                                                                                                                                  aatCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATgg                >  2:1231621/1‑150 (MQ=255)
                                                                                                                                   atCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGgc               <  1:1187978/150‑1 (MQ=255)
                                                                                                                                    tCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACTGACATGgcg              <  1:1438502/150‑1 (MQ=255)
                                                                                                                                    tCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGgcg              <  1:1175683/150‑1 (MQ=255)
                                                                                                                                     cTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTggg                                >  1:1228544/1‑131 (MQ=255)
                                                                                                                                     cTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTggg                                <  2:1228544/131‑1 (MQ=255)
                                                                                                                                     cTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGgcgc             >  2:899492/1‑150 (MQ=255)
                                                                                                                                     cTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGgcgc             <  2:922019/150‑1 (MQ=255)
                                                                                                                                       gCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGGCGCaa           >  1:1271423/1‑150 (MQ=255)
                                                                                                                                         aGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGGCGCaaaa         >  2:585510/1‑150 (MQ=255)
                                                                                                                                            cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCTCTGGGCGCACAACCGACATGGCGCAAAAACa      >  2:111837/1‑150 (MQ=255)
                                                                                                                                              gcgcAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGt    >  1:15818/1‑150 (MQ=255)
                                                                                                                                                gcAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGttt  >  1:409390/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCTATACCGATGAGGCGCGCGTCTCCCGCTTCTGGCGTGATGTCCGTTGTGAACGTATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTT  >  NC_000913/1778216‑1778509

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: