Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,488,057 C→G 7.5% intergenic (+22/‑175) ydcF → / → aldA DUF218 superfamily protein, SAM‑binding/aldehyde dehydrogenase A, NAD‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,488,0570CG 7.5% 388.9 / 21.8 134intergenic (+22/‑175)ydcF/aldADUF218 superfamily protein, SAM‑binding/aldehyde dehydrogenase A, NAD‑linked
Reads supporting (aligned to +/- strand):  ref base C (56/67);  new base G (5/5);  total (61/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-02

Only 100 of 135 total aligned reads displayed.

TACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGAAAACA  >  NC_000913/1487910‑1488206
                                                                                                                                                   |                                                                                                                                                     
tACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACtt                                                                                                                                                     >  1:1925512/1‑150 (MQ=255)
   gggccccGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGtt                                                                                                                                                  <  2:376720/147‑1 (MQ=255)
   ggTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGGCATCCATGCATGGCAAACGCGGAACCATGATGCGGAGCTCATCGAGGCGATGGAAAGTCGCGCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGtt                                                                                                                                                  <  2:122932/150‑1 (MQ=255)
     gccccGCGGGCGAGATTTTATCGTTCACCTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCGGAAACCAGATCCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTAACGTTAAAAATTGCCCGTATGTGAACCACTTGtttt                                                                                                                                                <  2:103077/149‑2 (MQ=255)
     gacccGAGGGCGAGATTTTATCGTTCACGTTGATTATCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTg                                                                                                                                                <  1:175786/148‑1 (MQ=255)
         cgcgGGCGAGATGTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGATCATCGAGTCGATGGCAAGTCGCTCATTACGTTAAAAATTGCCCGTTTGTGAACAACTTGTTTTCaaa                                                                                                                                            >  2:712833/1‑150 (MQ=255)
             ggCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACggg                                                                                                                                        <  1:951014/150‑1 (MQ=255)
              gCGAGCTTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACGGGc                                                                                                                                       <  2:857597/150‑1 (MQ=255)
              gCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGc                                                                                                                                       >  2:247987/1‑150 (MQ=255)
              gCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGc                                                                                                                                       >  1:1756924/1‑150 (MQ=255)
              gCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGGCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACGGGc                                                                                                                                       <  2:1302025/150‑1 (MQ=255)
               cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                                      <  1:419336/150‑1 (MQ=255)
               cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                                      >  1:770088/1‑150 (MQ=255)
               cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                                      >  2:1525494/1‑150 (MQ=255)
               cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                                      >  2:1154737/1‑150 (MQ=255)
               cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCa                                                                                                                                      <  1:978506/150‑1 (MQ=255)
               cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCACGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                                      >  1:532335/1‑150 (MQ=255)
               cGAGATTTTATCGTTCACGTTGATTTTCAGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                                      >  1:1391465/1‑150 (MQ=255)
                gagaTTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCAt                                                                                                                                     <  2:1691225/150‑1 (MQ=255)
                gagaTTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCaa                                                                                                                                     <  1:803388/150‑2 (MQ=39)
                 agaTTTTAGCGTTCACGTTGATTTTCCGGCAGAATTCATCCATTAATGGCAAACGCTGAAACATGTTGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATg                                                                                                                                    <  2:167390/150‑1 (MQ=255)
                    ttttATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACt                                                                                                                                 >  2:671663/1‑150 (MQ=255)
                     tttATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCGCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTc                                                                                                                                <  2:1567750/150‑1 (MQ=255)
                       tATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGTGCATGACTCCt                                                                                                                              >  1:1019152/1‑150 (MQ=255)
                        aTCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTg                                                                                                                             >  2:1565603/1‑150 (MQ=255)
                          cGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGAc                                                                                                                           >  2:23828/1‑150 (MQ=255)
                          cGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGAc                                                                                                                           >  2:1624303/1‑150 (MQ=255)
                          cGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGAc                                                                                                                           >  2:1700435/1‑150 (MQ=255)
                           gTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCCTCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACt                                                                                                                          <  2:358873/150‑1 (MQ=255)
                           gTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACt                                                                                                                          <  2:844711/150‑1 (MQ=255)
                            ttCCCGTGGATTTTCCGGCTGACGGCCTCCCTGCCTGGCACACGCTGAAACATGCGGCGTGGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGGTTGGGAACCCCTTGTTTGCAAACGGGCATGACTCCTGACtt                                                                                                                         <  1:110090/150‑1 (MQ=255)
                               aCGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTCCAAACGGGCATGACTCCTGACTTTTa                                                                                                                      <  1:1229034/150‑1 (MQ=255)
                                 gTTGGTTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTAtt                                                                                                                    >  2:816176/1‑150 (MQ=255)
                                 gTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTAtt                                                                                                                    <  1:1894444/150‑1 (MQ=255)
                                  ttGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTAttt                                                                                                                   >  2:1716094/1‑150 (MQ=255)
                                   tGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTc                                                                                                                  <  2:381930/150‑1 (MQ=255)
                                     tttttCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTg                                                                                                                <  2:458278/149‑1 (MQ=255)
                                      ttttCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGc                                                                                                               >  1:734620/1‑150 (MQ=255)
                                      ttttCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGc                                                                                                               <  2:1821590/150‑1 (MQ=255)
                                           cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                          <  2:2004374/150‑1 (MQ=255)
                                           cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                          >  1:1927644/1‑150 (MQ=255)
                                           cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                          >  1:1033438/1‑150 (MQ=255)
                                                 aaGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGTTGGTAAGTCGCTCGTTACGTTAAAAATTGCCGGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCt                                                                                                    <  1:1609657/150‑1 (MQ=255)
                                                    tcctCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGTTGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtttt                                                                                                 <  1:1465679/147‑1 (MQ=255)
                                                    tCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtttt                                                                                                 <  2:1546078/150‑1 (MQ=255)
                                                     cATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTa                                                                                                >  2:1247664/1‑150 (MQ=255)
                                                      aTCCATGCATGGCAAACGCTGAAACATTATGCGGTGCTCATCGAGGCGATGTAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTac                                                                                               <  1:1021185/150‑1 (MQ=255)
                                                      aTCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTac                                                                                               >  1:1655354/1‑150 (MQ=255)
                                                             catgGCAAACGCTGAAACATGATGCGGTTCTCTTCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGtt                                                                                        <  2:1039777/150‑1 (MQ=255)
                                                             catgGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGtt                                                                                        >  1:1911411/1‑150 (MQ=255)
                                                                ggCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGttttt                                                                                     >  2:1170409/1‑150 (MQ=255)
                                                                   aaaCGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATg                                                                                  <  2:503030/150‑1 (MQ=255)
                                                                   aaaCGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATg                                                                                  >  1:679079/1‑150 (MQ=255)
                                                                                tgatgCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAg                                                                     <  2:179582/150‑1 (MQ=255)
                                                                                    gCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTTCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCACTTCACAGAAtt                                                                 >  2:2017116/1‑150 (MQ=255)
                                                                                     cGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTg                                                                >  2:862226/1‑150 (MQ=255)
                                                                                         gCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCt                                                            >  2:1381225/1‑150 (MQ=255)
                                                                                           tCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTGCTTTTACACTTGTTTTTATGATGCCCTTCACAGAATTGTCCttt                                                          >  1:1792873/1‑150 (MQ=255)
                                                                                           tCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCttt                                                          >  1:1784786/1‑150 (MQ=255)
                                                                                                gAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGa                                                     >  1:1670957/1‑150 (MQ=255)
                                                                                                   gCGAGGGAAAGTCGCTCGTTACTTTAAAAATTGACAGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTATTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTc                                                  <  2:327258/150‑1 (MQ=255)
                                                                                                     gATGGAAAGTGGCTAGTTACGTTAAAAATTGCCAGTTTGTGTACCACTTGTTTTCAATCGGTGATGACTCCTTACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCg                                                <  1:1208901/150‑1 (MQ=255)
                                                                                                      aTGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGt                                               >  1:1616995/1‑150 (MQ=255)
                                                                                                       tGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                              >  2:1827172/1‑150 (MQ=255)
                                                                                                       tGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                              <  1:1025909/150‑1 (MQ=255)
                                                                                                       tGGAAAGTCGCGCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTTCCTTTTATTCCTTTTACACTTTTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                              <  2:1391465/150‑1 (MQ=255)
                                                                                                         gAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACGTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtctc                                            >  2:186477/1‑150 (MQ=255)
                                                                                                          aaaGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtctct                                           >  2:421170/1‑150 (MQ=255)
                                                                                                            aGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTctctct                                         >  1:695421/1‑150 (MQ=255)
                                                                                                             gTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtg                                        <  2:1019152/150‑1 (MQ=255)
                                                                                                             gTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGGTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtg                                        <  1:846704/150‑1 (MQ=255)
                                                                                                              tCGCTCGTTACGTTAAAAATTTCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtga                                       <  1:851113/150‑1 (MQ=255)
                                                                                                              tCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtga                                       >  2:1708835/1‑150 (MQ=255)
                                                                                                                gCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGATTTGTCCTTTCACGATTCCGTCTCTCtgatg                                     <  1:1247664/150‑1 (MQ=255)
                                                                                                                gCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgatg                                     <  1:1446968/150‑1 (MQ=255)
                                                                                                                gCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgatg                                     >  1:167730/1‑150 (MQ=255)
                                                                                                                 gtCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgatga                                    <  2:558865/149‑1 (MQ=255)
                                                                                                                      tacgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGAtgattgat                               <  1:1227181/150‑1 (MQ=255)
                                                                                                                       acgttaAAAATTGCCCTTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGAATTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATg                              <  2:48851/150‑1 (MQ=255)
                                                                                                                       acgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATg                              <  2:770088/150‑1 (MQ=255)
                                                                                                                       acgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTATTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATg                              >  1:1445908/1‑150 (MQ=255)
                                                                                                                        cgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTGGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGt                             >  1:1757431/1‑150 (MQ=255)
                                                                                                                         gttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGtt                            >  1:1096591/1‑150 (MQ=255)
                                                                                                                          ttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGtta                           >  2:824107/1‑150 (MQ=255)
                                                                                                                            aaaaaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaat                         <  2:1911411/150‑1 (MQ=255)
                                                                                                                              aaaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACACCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaatta                       <  2:1961093/150‑1 (MQ=255)
                                                                                                                                catGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAAc                     >  2:1497431/3‑150 (MQ=38)
                                                                                                                                catGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTAATACTTTTACACATGTTTTTTTGAAGCCCTTCACAGAATAGTCATTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAAc                     >  2:1497433/3‑150 (MQ=37)
                                                                                                                                 ttGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACa                    <  2:1082044/150‑1 (MQ=255)
                                                                                                                                  tGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACTATTCCGTCTCTCTGATGATTGATGTTAATTAACaa                   >  1:948089/1‑150 (MQ=255)
                                                                                                                                  tGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAActa                   >  1:362986/1‑148 (MQ=255)
                                                                                                                                    cccGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATg                 >  1:1384632/1‑150 (MQ=255)
                                                                                                                                     ccGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGt                <  2:1756924/150‑1 (MQ=255)
                                                                                                                                      cGTTTGGGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTa               <  2:887121/150‑1 (MQ=255)
                                                                                                                                       gTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCATTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTAt              >  1:616630/1‑150 (MQ=255)
                                                                                                                                          tgtgAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCa           >  2:1385897/1‑150 (MQ=255)
                                                                                                                                           gtgAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTACTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTAAACGATTCCGTCTCTCTGATGATTGATGTTAATTAACATTGTATTCAc          >  2:338427/1‑150 (MQ=255)
                                                                                                                                             gAACCACTTGTTTGCAAACGGGAATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCg        <  2:1843684/150‑1 (MQ=255)
                                                                                                                                                  aCTTTTTTGCTAACGGGCATGACTCCTGACTTTTATTTCTTCCTTTTATTCCTTTTAAACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGAaaac   <  2:1497831/150‑1 (MQ=255)
                                                                                                                                                   cTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGAaaaca  >  2:908629/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCGAAAACA  >  NC_000913/1487910‑1488206

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