Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,755,946 T→C 6.0% N189S (AAC→AGC)  aldB ← aldehyde dehydrogenase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,755,9460TC 6.0% 197.2 / 3.0 67N189S (AAC→AGC) aldBaldehyde dehydrogenase B
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base C (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

CGTTTCGAGGTCGCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAACTTCACTGATCCCACCTTCCTGC  >  NC_000913/3755799‑3756092
                                                                                                                                                   |                                                                                                                                                  
cGTTTCGAGGTCGCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGt                                                                                                                                                    <  2:244552/148‑1 (MQ=255)
 gTTTCGAGGTCGCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTg                                                                                                                                                  >  8:147681/1‑149 (MQ=255)
  tttCGAGGTCGCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGc                                                                                                                                                 >  1:9432/1‑149 (MQ=255)
      gAGGTCGCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGAACCACACAGTTGCccgc                                                                                                                                             >  7:139698/1‑149 (MQ=255)
         gTCGCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCcgccgc                                                                                                                                          >  1:199117/1‑149 (MQ=255)
         gTCGCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCcgccgc                                                                                                                                          >  8:62191/1‑149 (MQ=255)
          tCGCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAACTCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGGTTCAGCACCACACAGTTGCCCgccgcc                                                                                                                                         <  7:219559/149‑1 (MQ=255)
          tCGCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGt                                                                                                                                                    <  8:144888/138‑1 (MQ=255)
          tCGCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGt                                                                                                                                                    >  7:144888/1‑138 (MQ=255)
            gCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACAACGTTCACCACGCCCGGCGTCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTAGCCCGCCGCCAg                                                                                                                                       >  1:209869/1‑149 (MQ=255)
                 atatTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCCCCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCgg                                                                                                                                  <  1:103790/149‑1 (MQ=255)
                 atatTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTACATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCgg                                                                                                                                  <  2:106555/149‑1 (MQ=255)
                    ttCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAg                                                                                                                               >  3:178749/1‑149 (MQ=255)
                         cAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCCGCGCGGGAGCCAtt                                                                                                                          >  3:167278/1‑149 (MQ=255)
                          aaTTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCAttt                                                                                                                         <  4:178749/149‑1 (MQ=255)
                          aaTTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCAttt                                                                                                                         >  3:50983/1‑149 (MQ=255)
                               ccccACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAg                                                                                                                    >  2:33889/1‑149 (MQ=255)
                                cccACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCAGTTTCCAGc                                                                                                                   >  5:75482/1‑149 (MQ=255)
                                     cTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGcc                                                                                                              <  1:53331/149‑1 (MQ=255)
                                           cATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCACCAGTACAGAAAGCGCGGTAAGACGTGCGGGTTTCAGCACCACACAGTAGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATcag                                                                                                         >  5:63972/1‑148 (MQ=255)
                                            aTTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATcagca                                                                                                       <  7:73157/149‑1 (MQ=255)
                                                aCCCCGTTCACCACGCCCCGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTTAGACGTGCGGGGTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCCTCAGCAGCg                                                                                                    <  6:191775/148‑1 (MQ=255)
                                                aCCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCg                                                                                                    <  4:50983/148‑1 (MQ=255)
                                                  cACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCggg                                                                                                  >  5:154058/1‑148 (MQ=255)
                                                    cgtgctccacGCCCGGCGGCAGTATATCACCGACAATTTCGATTAGCAGCAGTTCAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCGATTTTCCAGCTCGCCATCAGCAGCGGGAAg                                                                                               <  3:146193/143‑1 (MQ=255)
                                                    cGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAg                                                                                               <  7:164671/149‑1 (MQ=255)
                                                           cacGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACg                                                                                        >  7:55112/1‑149 (MQ=255)
                                                             cGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACggg                                                                                      <  2:9432/149‑1 (MQ=255)
                                                                ccGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGAt                                                                                    >  2:178123/1‑148 (MQ=255)
                                                                   gcggcAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATc                                                                                <  3:33562/149‑1 (MQ=255)
                                                                   gcggcAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGGTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATc                                                                                <  3:30362/149‑1 (MQ=255)
                                                                       cAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGcc                                                                            <  2:199117/149‑1 (MQ=255)
                                                                         gTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCAtt                                                                                                                          >  1:59914/1‑101 (MQ=255)
                                                                         gTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCAtt                                                                                                                          <  2:59914/101‑1 (MQ=255)
                                                                            aaTCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCccac                                                                        >  2:60068/1‑148 (MQ=255)
                                                                                aCCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGc                                                                   <  5:205624/149‑1 (MQ=255)
                                                                                  cGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGcc                                                                  >  6:149421/1‑148 (MQ=255)
                                                                                   gACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCt                                                                 >  4:178066/1‑148 (MQ=255)
                                                                                    aCAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTaa                                                               <  6:154058/149‑1 (MQ=255)
                                                                                     cAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAAc                                                              <  8:55112/149‑1 (MQ=255)
                                                                                        tttCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGt                                                           <  8:17265/149‑1 (MQ=255)
                                                                                           ccATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACgg                                                            <  6:111527/145‑1 (MQ=255)
                                                                                           ccATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACgg                                                            >  5:111527/1‑145 (MQ=255)
                                                                                            cATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCa                                                        <  7:147681/148‑1 (MQ=255)
                                                                                              ttAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATg                                                      >  5:80700/1‑148 (MQ=255)
                                                                                                 gcagcaGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGaaa                                                  >  4:15108/1‑149 (MQ=255)
                                                                                                    gcagTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGa                                               >  7:13560/1‑149 (MQ=255)
                                                                                                        tACAGATAGCGGGGTAAGACTTGCGGGTTTCAGCACCAGACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGAGTGAAGTTCCACGGTATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAgg                                           <  2:209869/149‑1 (MQ=255)
                                                                                                        tACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAgg                                           <  5:149421/149‑1 (MQ=255)
                                                                                                        tACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGCTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCccacca                                                                      >  8:35352/1‑122 (MQ=255)
                                                                                                        tACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGCTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCccacca                                                                      >  8:31367/1‑122 (MQ=255)
                                                                                                        tACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGCTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCccacca                                                                      <  7:31367/122‑1 (MQ=255)
                                                                                                        tACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGCTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCccacca                                                                      <  7:35352/122‑1 (MQ=255)
                                                                                                                gCGGGCAAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCAAATTCAAGCTCGCCATCAGCAGCGGGAAGATCCACGGGATAATCTGCCCCACCACGCCTAACGGATCATGGAAATGATAGGCCCCggat                                   >  2:58057/1‑147 (MQ=255)
                                                                                                                  ggggTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTc                                 <  3:15108/149‑1 (MQ=255)
                                                                                                                    ggTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGc                               <  6:80700/149‑1 (MQ=255)
                                                                                                                            gTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAACt                       >  3:35429/1‑149 (MQ=255)
                                                                                                                                  gTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAACTTCACt                  >  4:208002/1‑148 (MQ=255)
                                                                                                                                   tttCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAACTTCACTGa                <  1:60068/149‑1 (MQ=255)
                                                                                                                                     tCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAACTTCACTGAt               <  1:33889/148‑1 (MQ=255)
                                                                                                                                        gCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAAAGGTTCATGGAAATGATAGGCCACGGTTTCGCTATAAACTACACGCATcccc           >  8:239208/1‑148 (MQ=255)
                                                                                                                                         caccacACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAACTTCACTGATCCCAc          <  8:13560/149‑1 (MQ=255)
                                                                                                                                         caccacACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAACTTCACTGATCCCAc          >  3:184566/1‑149 (MQ=255)
                                                                                                                                            cacacaGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAACTTCACTGATCCCACCtt       >  7:14948/1‑149 (MQ=255)
                                                                                                                                            cacacaGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAACTTCACTGATCCCACCtt       >  1:70966/1‑149 (MQ=255)
                                                                                                                                             acacaGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAAc                        <  1:69140/131‑1 (MQ=255)
                                                                                                                                             acacaGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAAc                        >  2:69140/1‑131 (MQ=255)
                                                                                                                                              cacaGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAACTTCACTGATCCCCCCTTcc     <  6:63972/149‑1 (MQ=255)
                                                                                                                                                caGTTGCCCGCCGCCAGCGCGGGCGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATCGGCCACGGTTTCGCTATCAACTTCACTGATCCCACCTTCCTg   <  7:239208/149‑1 (MQ=255)
                                                                                                                                                 aGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCCCGCCTAACGGTTCATGGAAATGATAGGCCCCGGTTTCGCTATCACCTTCAAGGATCCCACCTTCATgc  >  2:209112/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGTTTCGAGGTCGCCAGATATTCGCCAATTACCCCACCTGCGCCATTGACCACGTTCACCACGCCCGGCGGCAGTAAATCACCGACAATTTCCATTAGCAGCAGTACAGAAAGCGGGGTAAGACGTGCGGGTTTCAGCACCACACAGTTGCCCGCCGCCAGCGCGGGAGCCATTTTCCAGCTCGCCATCAGCAGCGGGAAGTTCCACGGGATAATCTGCCCCACCACGCCTAACGGTTCATGGAAATGATAGGCCACGGTTTCGCTATCAACTTCACTGATCCCACCTTCCTGC  >  NC_000913/3755799‑3756092

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: ‑