Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,488,050 T→C 7.7% intergenic (+15/‑182) 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,0500TC 7.7% 366.3 / 17.0 130intergenic (+15/‑182)ydcF/aldADUF218 superfamily protein, SAM‑binding/aldehyde dehydrogenase A, NAD‑linked
Reads supporting (aligned to +/- strand):  ref base T (55/63);  new base C (5/5);  total (60/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

Only 100 of 136 total aligned reads displayed.

GCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCAC  >  NC_000913/1487902‑1488198
                                                                                                                                                    |                                                                                                                                                    
gCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGAGGTGCTCATCGAGGCGATGCAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTtgtg                                                                                                                                                     >  1:1776552/1‑150 (MQ=255)
  gATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGaa                                                                                                                                                   >  2:75175/1‑150 (MQ=255)
   aTGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAAc                                                                                                                                                  <  2:1719914/150‑1 (MQ=255)
    tGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAAcc                                                                                                                                                 <  1:1025561/150‑1 (MQ=255)
        tACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACtt                                                                                                                                             >  1:1925512/1‑150 (MQ=255)
         aCGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTg                                                                                                                                            <  1:720088/150‑1 (MQ=255)
           ggTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGGCATCCATGCATGGCAAACGCGGAACCATGATGCGGAGCTCATCGAGGCGATGGAAAGTCGCGCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGtt                                                                                                                                          <  2:122932/150‑1 (MQ=255)
             tCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTg                                                                                                                                        <  2:1564172/150‑1 (MQ=255)
             gccccGCGGGCGAGATTTTATCGTTCACCTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCGGAAACCAGATCCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTAACGTTAAAAATTGCCCGTATGTGAACCACTTGtttt                                                                                                                                        <  2:103077/149‑2 (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                                                                                                                              >  2:1154737/1‑150 (MQ=255)
                       cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                              <  1:419336/150‑1 (MQ=255)
                       cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                              >  1:770088/1‑150 (MQ=255)
                       cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCa                                                                                                                              <  1:978506/150‑1 (MQ=255)
                       cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCACGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                              >  1:532335/1‑150 (MQ=255)
                        gagaTTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCAt                                                                                                                             <  2:1691225/150‑1 (MQ=255)
                        gagaTTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCaa                                                                                                                             <  1:803388/150‑2 (MQ=39)
                         agaTTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGGTACGGTAAAAATTGCCCGTGTGTGAACCACTTGTTTGCAAACGGGCATg                                                                                                                            >  2:834662/1‑150 (MQ=255)
                           aTTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGAc                                                                                                                          <  1:1660298/150‑1 (MQ=255)
                            ttttATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACt                                                                                                                         >  2:671663/1‑150 (MQ=255)
                               tATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGTGCATGACTCCt                                                                                                                      >  1:1019152/1‑150 (MQ=255)
                               tATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCGCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCt                                                                                                                      <  2:1062310/150‑1 (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)
                                      cACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACtttt                                                                                                               >  2:864890/1‑150 (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)
                                              ttttCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGc                                                                                                       >  1:734620/1‑150 (MQ=255)
                                              ttttCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGc                                                                                                       <  2:1821590/150‑1 (MQ=255)
                                                   cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                  >  1:1927644/1‑150 (MQ=255)
                                                   cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                  <  2:2004374/150‑1 (MQ=255)
                                                   cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                  >  1:1033438/1‑150 (MQ=255)
                                                      aataaGTAATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCGCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAACCGGGCATGACTCCTGACTTTTATTTCTGCCTTTTAtt                                                                                               <  1:1233766/147‑1 (MQ=255)
                                                         aaGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGTTGGTAAGTCGCTCGTTACGTTAAAAATTGCCGGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCt                                                                                            <  1:1609657/150‑1 (MQ=255)
                                                            tcctCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGTTGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtttt                                                                                         <  1:1465679/147‑1 (MQ=255)
                                                            tCATCCATGCATGGCAAACGCTTAAACATTATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtttt                                                                                         <  1:1420647/150‑1 (MQ=255)
                                                            tCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtttt                                                                                         <  1:247987/150‑1 (MQ=255)
                                                            tCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtttt                                                                                         <  2:1546078/150‑1 (MQ=255)
                                                             cATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTa                                                                                        >  1:1707950/1‑150 (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)
                                                                        ggCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGttttt                                                                             >  2:1170409/1‑150 (MQ=255)
                                                                           aaaCGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATg                                                                          <  2:503030/150‑1 (MQ=255)
                                                                              cGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAg                                                                       <  2:424502/150‑1 (MQ=255)
                                                                                cTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGcc                                                                     >  2:366023/1‑150 (MQ=255)
                                                                                        tgatgCGGTGCTCTTCGTGGCGATGGAATGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAg                                                             <  1:1768472/150‑1 (MQ=255)
                                                                                        tgatgCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAg                                                             <  2:179582/150‑1 (MQ=255)
                                                                                            gCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTTCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCACTTCACAGAAtt                                                         >  2:2017116/1‑150 (MQ=255)
                                                                                                 gCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCt                                                    >  2:1381225/1‑150 (MQ=255)
                                                                                                   tCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTGCTTTTACACTTGTTTTTATGATGCCCTTCACAGAATTGTCCttt                                                  >  1:1792873/1‑150 (MQ=255)
                                                                                                   tCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCttt                                                  <  1:1908485/150‑1 (MQ=255)
                                                                                                    cATCGAGGCGATGGAACGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCCCTTGTTTGCAAACGGGCATGACTCCTGACTTTTAGTTCTGCGGTTTATTACTTTTACACTTGTTTTTATGAAGCCCTTCACATCATTGTCCTTTc                                                 >  1:1204837/1‑150 (MQ=255)
                                                                                                      tCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCAc                                               <  1:1547464/150‑1 (MQ=255)
                                                                                                        gAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGa                                             >  1:1670957/1‑150 (MQ=255)
                                                                                                           gCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTc                                          <  2:1655354/150‑1 (MQ=255)
                                                                                                             gATGGAAAGTGGCTAGTTACGTTAAAAATTGCCAGTTTGTGTACCACTTGTTTTCAATCGGTGATGACTCCTTACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCg                                        <  1:1208901/150‑1 (MQ=255)
                                                                                                               tGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                      >  2:1827172/1‑150 (MQ=255)
                                                                                                               tGGAAAGTCGCGCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTTCCTTTTATTCCTTTTACACTTTTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                      <  2:1391465/150‑1 (MQ=255)
                                                                                                                 gAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACGTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtctc                                    >  2:186477/1‑150 (MQ=255)
                                                                                                                  aaaGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtctct                                   >  2:421170/1‑150 (MQ=255)
                                                                                                                    aGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTctctct                                 <  2:1707293/150‑1 (MQ=255)
                                                                                                                    aGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTctctct                                 >  1:695421/1‑150 (MQ=255)
                                                                                                                     gTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGGTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtg                                <  1:846704/150‑1 (MQ=255)
                                                                                                                      tCGCTCGTTACGTTAAAAATTTCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtga                               <  1:851113/150‑1 (MQ=255)
                                                                                                                        gCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGATTTGTCCTTTCACGATTCCGTCTCTCtgatg                             <  1:1247664/150‑1 (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:264344/150‑1 (MQ=255)
                                                                                                                                    aaaaaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaat                 <  2:1911411/150‑1 (MQ=255)
                                                                                                                                     aaaaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaatt                >  1:1029872/1‑150 (MQ=255)
                                                                                                                                      aaaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaatta               <  1:990607/150‑1 (MQ=255)
                                                                                                                                      aaaTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACACCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaatta               <  2:1961093/150‑1 (MQ=255)
                                                                                                                                        catGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAAc             >  2:1510078/3‑150 (MQ=38)
                                                                                                                                        catGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAAc             >  2:1497431/3‑150 (MQ=38)
                                                                                                                                        catGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTAATACTTTTACACATGTTTTTTTGAAGCCCTTCACAGAATAGTCATTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAAc             >  2:1497433/3‑150 (MQ=37)
                                                                                                                                            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)
                                                                                                                                                  tgtgAACCAATTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCa   <  1:864890/150‑1 (MQ=255)
                                                                                                                                                   gtgAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTACTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTAAACGATTCCGTCTCTCTGATGATTGATGTTAATTAACATTGTATTCAc  >  2:338427/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCAC  >  NC_000913/1487902‑1488198

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: