Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,367,774 (C)6→7 intergenic (+132/‑53) yhcF → / → yhcG DUF1120 domain‑containing protein YhcF/DUF1016 domain‑containing protein YhcG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,367,7681.C94.6% 187.5 / 2.7 56intergenic (+126/‑59)yhcF/yhcGDUF1120 domain‑containing protein YhcF/DUF1016 domain‑containing protein YhcG
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (26/27);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCG  >  NC_000913/3367629‑3367896
                                                                                                                                            |                                                                                                                                   
tAGTCTACTTGTGATAACAATGTCGGACAGGAGAACTATCAGGCATTCTTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCa                                                                                                                                      <  2:400487/140‑1 (MQ=255)
 aGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAc                                                                                                                                    <  1:178587/140‑1 (MQ=255)
   tCTACTTGTGATAACAAGGCCGAACAGGAGTACTACCAGGCATTCGTTCTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAccc                                                                                                                                <  2:262486/140‑1 (MQ=255)
      aCTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAcccccc                                                                                                                            <  1:684419/140‑1 (MQ=255)
      aCTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAcccccc                                                                                                                            <  1:785573/140‑1 (MQ=255)
       cTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAccccccc                                                                                                                           >  1:547284/1‑139 (MQ=255)
       cTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAccccccc                                                                                                                           <  1:424914/140‑2 (MQ=255)
        ttGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAccccccca                                                                                                                          >  2:9343/1‑138 (MQ=255)
        ttGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAccccccca                                                                                                                          <  2:543490/140‑3 (MQ=255)
        ttGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAccccccca                                                                                                                          <  1:616018/140‑3 (MQ=255)
        ttGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAccccccca                                                                                                                          <  1:624135/140‑3 (MQ=255)
            gATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGt                                                                                                                       >  1:106392/1‑140 (MQ=255)
             aTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGta                                                                                                                      >  1:792839/1‑140 (MQ=255)
             aTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGta                                                                                                                      >  2:119238/1‑140 (MQ=255)
             aTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAccc                                                                                                                                >  1:201318/1‑130 (MQ=255)
             aTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCAccc                                                                                                                                <  2:201318/130‑1 (MQ=255)
                  aaTGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCt                                                                                                                 <  2:502658/140‑1 (MQ=255)
                         gACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCa                                                                                                          <  1:9343/140‑1 (MQ=255)
                         gACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCa                                                                                                          <  1:641543/140‑1 (MQ=255)
                         gACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCa                                                                                                          <  2:618958/140‑1 (MQ=255)
                              gAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCattatt                                                                                                     <  1:119238/140‑1 (MQ=255)
                              gAGTACTATCAGGCATTCGTTTCATACGATCTGTTTATGCTTCAGGAAGGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCattatt                                                                                                     <  2:232837/140‑1 (MQ=255)
                                   cTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACt                                                                                                <  1:708759/140‑1 (MQ=255)
                                        aGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAg                                                                                             >  1:502798/1‑138 (MQ=255)
                                        aGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAg                                                                                             <  2:502798/138‑1 (MQ=255)
                                        aGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCt                                                                                           >  1:40770/1‑140 (MQ=255)
                                        aGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCt                                                                                           >  2:801019/1‑140 (MQ=255)
                                        aGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCACCCCAGAGTATGGCGCATTCCGTTATTCAACTTAGCt                                                                                           >  1:232837/1‑140 (MQ=255)
                                           cATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGc                                                                                            >  1:652530/1‑136 (MQ=255)
                                           cATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGc                                                                                            <  2:652530/136‑1 (MQ=255)
                                           cATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCa                                                                                        <  2:40770/140‑1 (MQ=255)
                                           cATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGCGCTCACCCCCCCCGAGCATAGCTCATTCCATTATTCAACTTAGCTCCa                                                                                        >  2:296172/1‑140 (MQ=255)
                                             ttCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAAAACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCata                                                                                      >  2:519907/1‑140 (MQ=255)
                                                    tataCGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAggg                                                                               <  2:792839/140‑1 (MQ=255)
                                                           tCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATa                                                                        >  2:522502/1‑140 (MQ=255)
                                                               ttAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATgg                                                                    >  2:183106/1‑140 (MQ=255)
                                                                 aaTGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAg                                                                                 <  2:365828/125‑1 (MQ=255)
                                                                 aaTGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAg                                                                                 >  1:365828/1‑125 (MQ=255)
                                                                  aTGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAAt                                                                 <  1:801019/140‑1 (MQ=255)
                                                                     tttCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATaa                                                                       >  2:594887/1‑131 (MQ=255)
                                                                     tttCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATaa                                                                       <  1:594887/131‑1 (MQ=255)
                                                                       tCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGctc                                                            >  2:802794/1‑140 (MQ=255)
                                                                       tCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGctc                                                            >  1:170549/1‑140 (MQ=255)
                                                                       tCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGctc                                                            >  1:117749/1‑140 (MQ=255)
                                                                         aGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACGGATAATGGAATCGCtctc                                                          >  2:118168/1‑140 (MQ=255)
                                                                         aGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAAATTAGCTCGATATAAAGTGACACATAATGGAATCGCtctc                                                          >  2:318330/1‑140 (MQ=255)
                                                                            aaaGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGa                                                       <  1:522502/140‑1 (MQ=255)
                                                                             aaGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGaa                                                      <  2:106392/140‑1 (MQ=255)
                                                                             aaGCTCTTCGCATATCTTGGGGAAAGTAAATGTAAAACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGaa                                                      <  1:519907/140‑1 (MQ=255)
                                                                                             ttGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCt                                      >  1:785503/1‑140 (MQ=255)
                                                                                                 ggAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCa                                  <  1:802794/140‑1 (MQ=255)
                                                                                                     aGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCaa                              >  1:112814/1‑140 (MQ=255)
                                                                                                        aaaTGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATc                           <  2:112814/140‑1 (MQ=255)
                                                                                                           ggTAATACGCTGTTCTTAATTGCTCCGTGTTCACCCCGCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCac                        <  1:259556/139‑1 (MQ=255)
                                                                                                             tAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCacga                      <  2:547284/140‑1 (MQ=255)
                                                                                                             tAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCacga                      <  2:768434/140‑1 (MQ=255)
                                                                                                             tAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCacga                      <  2:785503/140‑1 (MQ=255)
                                                                                                             tAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCacga                      <  1:190220/140‑1 (MQ=255)
                                                                                                             tAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCacga                      <  1:183106/140‑1 (MQ=255)
                                                                                                                tACGCTGTTCTTAATTGTTCAGTGCTCACACCCCCAGCGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACgg                   <  1:296172/140‑1 (MQ=255)
                                                                                                                tACGCTGTTCTTAATTGTTCAGCGCTCACCCCCCCAGTGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACgg                   <  1:318330/140‑1 (MQ=255)
                                                                                                                 aCGCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGt                  >  1:571531/1‑140 (MQ=255)
                                                                                                                   gCTGTTCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGtat                <  1:825939/140‑1 (MQ=255)
                                                                                                                       ttCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTGTCAGAAGGGACCGCAGTAGGCTACCTGCAATTCCACGACGGtattatc            >  2:259556/1‑139 (MQ=255)
                                                                                                                       ttCTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTattatt            >  1:564170/1‑140 (MQ=255)
                                                                                                                         cTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCacga                      >  1:686769/1‑128 (MQ=255)
                                                                                                                         cTTAATTGTTCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCacga                      <  2:686769/128‑1 (MQ=255)
                                                                                                                                 ttCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCg  >  2:665828/1‑140 (MQ=255)
                                                                                                                                 ttCAGTGCTCACCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCg  >  2:806324/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                   
TAGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGTTAATGTTTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCG  >  NC_000913/3367629‑3367896

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: