Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,697,242 Δ1 bp 6.1% intergenic (‑190/+31) uidR ← / ← hdhA transcriptional repressor/7‑alpha‑hydroxysteroid dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,697,2420A. 6.1% 187.9 / 6.5 66intergenic (‑190/+31)uidR/hdhAtranscriptional repressor/7‑alpha‑hydroxysteroid dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base . (2/2);  total (35/31)
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.97e-01

GTCAGGTGAGCGCCTATCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCGGTTGGCCCA  >  NC_000913/1697094‑1697388
                                                                                                                                                    |                                                                                                                                                   
gTCAGGTGAGCGCCTATCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAacac                                                                                                                                                    >  1:204017/1‑150 (MQ=255)
  cAGGTGAGCGCCTATCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAAcacgc                                                                                                                                                  >  1:514302/1‑150 (MQ=255)
           gcCTATCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑gttgtt                                                                                                                                        <  2:262532/150‑1 (MQ=255)
            ccTATCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGttt                                                                                                                                       >  2:384409/1‑150 (MQ=255)
            ccTATCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGttt                                                                                                                                       >  2:419886/1‑150 (MQ=255)
              tATCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTAc                                                                                                                                     >  2:1084332/1‑150 (MQ=255)
              tATCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGACAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCAAG‑GGTGTTTAc                                                                                                                                     >  1:985274/1‑150 (MQ=255)
               agcATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTAcc                                                                                                                                    <  1:734945/148‑1 (MQ=255)
               aTCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTAcc                                                                                                                                    <  1:559205/150‑1 (MQ=255)
                tCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCg                                                                                                                                   <  1:1244841/150‑1 (MQ=255)
                tCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAAC‑CGTGCGTGTTGTTTACCg                                                                                                                                   <  1:900229/150‑1 (MQ=255)
                 cATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAAC‑CGTGCGTGTTGTTTACCgg                                                                                                                                  <  2:847777/150‑1 (MQ=255)
                 cATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCCTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCgg                                                                                                                                  <  1:887346/150‑1 (MQ=255)
                   ttGCTGTTGAATAGTGCGAATGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTc                                                                                                                                <  1:250375/150‑1 (MQ=255)
                       tGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGtt                                                                                                                            >  2:19315/1‑150 (MQ=255)
                       tGTAGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGATATTGGTTAACCTATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGGTTACCGGTCCGtt                                                                                                                            >  1:80084/1‑150 (MQ=255)
                          tGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACTACACGCACG‑GTTGTTTACCGGTCCGTTAgt                                                                                                                         >  1:312182/1‑150 (MQ=255)
                           gAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTGATTTAGCAGAGACGTTGATGAGTTTATCCCGGTCTTCTACACGCACG‑GTTGTTTACCGGTCCGTTtgtg                                                                                                                        >  2:824409/1‑150 (MQ=255)
                             aTAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTa                                                                                                                      >  2:471394/1‑150 (MQ=255)
                              tAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTAt                                                                                                                     >  2:31502/1‑150 (MQ=255)
                                  gCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTTATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAAt                                                                                                                 <  1:1310247/150‑1 (MQ=255)
                                    gAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTg                                                                                                               >  1:1054591/1‑150 (MQ=255)
                                    gAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTg                                                                                                               <  1:597855/150‑1 (MQ=255)
                                       ggCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGc                                                                                                            <  2:1247379/150‑1 (MQ=255)
                                            aatctATTCGTATACTAGTCACTCTGAAATTGGTGAACAAATAGCCTTGCTTTTGCTGCCAGTTATTTAGCAGAGACGATGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGATTACCGGTCCGTTAGTGTATTAATTGAGCTCCTg                                                                                                       <  2:61748/148‑1 (MQ=255)
                                                    cGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCccac                                                                                               >  1:1266261/1‑150 (MQ=255)
                                                     gTATACTTTTCACTCTGAAATTGGTTAACCAAAAGCCGTGATTTGGCTGCAAGTTATTTAGCAGAGACGTTGATGAGTTAATCACGGAAAACAACACGAACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCcaacc                                                                                              >  1:945567/1‑150 (MQ=255)
                                                        taCTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTGGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCaccacc                                                                                           >  2:93585/1‑150 (MQ=255)
                                                        taCTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGGTATTTAGCAGAGACGTTGATGAGTTTATCGCGTTACACAACCCGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGATCCTGTACCCCacctcc                                                                                           >  2:87096/1‑150 (MQ=255)
                                                        taCTATTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCaccacc                                                                                           <  1:361577/150‑1 (MQ=255)
                                                           ttttCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGga                                                                                        <  2:739856/150‑1 (MQ=255)
                                                             ttCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGgaga                                                                                      >  2:1157993/1‑150 (MQ=255)
                                                                    tGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTgag                                                                               >  2:44301/1‑150 (MQ=255)
                                                                         ttGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTTATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAAttt                                                                          <  1:660507/150‑1 (MQ=255)
                                                                          tGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTg                                                                         <  1:1247976/150‑1 (MQ=255)
                                                                           ggTTAACCAATAGCCGGGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGATCTCCTGTAACCCACCACCGGAGACGGTGAGAATTTGt                                                                        >  1:1274910/1‑150 (MQ=255)
                                                                             ttAAGCTATAGCCGTGTTTTGGCTGCGAGTTATTTAGCAGTGACGTTGATGAGTTTATCCCGGTAAACAACAGGCATG‑GTTGTTTTCCGGTCCGTTAGTGTATTAATTTAGCTCCTGTACCCCTCCACCAGTGACGGTGAGATTTTGtac                                                                      <  2:1090828/150‑3 (MQ=255)
                                                                               aaCCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGc                                                                    >  1:1023319/1‑150 (MQ=255)
                                                                                  cAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTa                                                                 <  2:49835/150‑1 (MQ=255)
                                                                                   aaTAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTAc                                                                <  1:949961/150‑1 (MQ=255)
                                                                                    aTAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTAcc                                                               <  1:63137/150‑1 (MQ=255)
                                                                                       gCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCATCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAg                                                            >  2:1299583/1‑150 (MQ=255)
                                                                                       gCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAg                                                            >  2:515498/1‑150 (MQ=255)
                                                                                        ccGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTAATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCAGTACCCCACCACCGGAGACGGTGAGAATTTTTCCGCTTACCCAGc                                                           >  1:124419/1‑150 (MQ=255)
                                                                                           tGATTTGCCTGCCAGTTATTTAGCTTAGACGTTGATGAGTTTATCCCGTTATAAAACACGCACG‑GTTGTTTACCGGTCCGTTTGTCTATTAATTCAGCTCCTGTACCACACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCg                                                        <  1:963778/150‑1 (MQ=255)
                                                                                            gATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCCCCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCgc                                                       >  1:297195/1‑150 (MQ=255)
                                                                                                  gCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCTGTAATCAACACGCACG‑GGTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAgg                                                 <  1:940211/150‑1 (MQ=255)
                                                                                                    tGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCg                                               >  1:910897/1‑150 (MQ=255)
                                                                                                     gCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCgc                                              >  2:232980/1‑149 (MQ=255)
                                                                                                      ccAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCCCG‑GTTGTTTACCGGTCCGTTAGAGTATTAATTGAGCTCCTGTACCCCACCAACGGAGACGGTGAGAATTTGTCCGCTAACCCAGCTCGCAGCAGGCGAg                                             >  2:823801/1‑150 (MQ=255)
                                                                                                      ccAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAg                                             <  1:603777/150‑1 (MQ=255)
                                                                                                       cAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGc                                            >  1:370548/1‑150 (MQ=255)
                                                                                                       cAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGc                                            >  1:303378/1‑150 (MQ=255)
                                                                                                        aGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAATCAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCCCCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCa                                           >  2:1252790/1‑150 (MQ=255)
                                                                                                             tttAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGa                                      <  2:514302/150‑1 (MQ=255)
                                                                                                               tAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAAc                                    <  2:265048/150‑1 (MQ=255)
                                                                                                                aGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACa                                   >  1:1077129/1‑150 (MQ=255)
                                                                                                                   agagACGTTGACGAGTTAATCCCGGTAACCAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTCTTAATTGAGCTCCTGTACCCCACCACCGTAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAgcg                                <  1:191929/150‑1 (MQ=255)
                                                                                                                         gTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGt                          >  2:1020625/1‑150 (MQ=255)
                                                                                                                             atgaGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGc                      >  1:398401/1‑150 (MQ=255)
                                                                                                                                aGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAat                   <  2:1049611/150‑1 (MQ=255)
                                                                                                                                 gTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAata                  <  1:102828/150‑1 (MQ=255)
                                                                                                                                        ccGGTAAACAAC‑CGTGCGTGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTA‑CCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCgg          >  1:839488/1‑150 (MQ=255)
                                                                                                                                         cGGTAAACAACACGCACG‑GTTTTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGTTGAGAATTTGTCCGCTTACCCATCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCgg          <  1:1100976/150‑1 (MQ=255)
                                                                                                                                         cGGTAAACAAC‑CGTGCGTGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCgg          >  2:822304/1‑150 (MQ=255)
                                                                                                                                         cGGAAAACAACACGCACG‑GTTGTTTACCGGTCCGTTCGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGACCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCgg          <  2:1223830/150‑1 (MQ=255)
                                                                                                                                           gTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCGGtt        >  1:1237706/1‑150 (MQ=255)
                                                                                                                                             aaaCAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTCAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCGGTTgg      <  2:204017/150‑1 (MQ=255)
                                                                                                                                               acaacaCGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCGGTTGGcc    <  2:811250/150‑1 (MQ=255)
                                                                                                                                               acaacaCGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCGGTTGGcc    >  1:1137607/1‑150 (MQ=255)
                                                                                                                                               acaacaCGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCGGTTGGcc    <  2:973972/150‑1 (MQ=255)
                                                                                                                                                caacaCGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCGGTTGGccc   <  2:910897/150‑1 (MQ=255)
                                                                                                                                                 aacaCGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCGGTTGGCCCa  <  2:465209/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCAGGTGAGCGCCTATCATTGCTGTTGAATAGTGCGAAGGCACACTCTATTCGTATACTTTTCACTCTGAAATTGGTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACG‑GTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCTTGCGGTTGGCCCA  >  NC_000913/1697094‑1697388

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