Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,608,704 A→G 10.5% S107P (TCG→CCG)  yneE ← bestrophin family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,608,7040AG10.5% 244.2 / 11.4 83S107P (TCG→CCG) yneEbestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base G (5/4);  total (42/41)
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

CTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGC  >  NC_000913/1608557‑1608850
                                                                                                                                                   |                                                                                                                                                  
cTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGaa                                                                                                                                                   <  3:511266/149‑1 (MQ=255)
  ggAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAAt                                                                                                                                                  >  1:438086/1‑148 (MQ=255)
        ttCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaa                                                                                                                                           >  3:251370/1‑149 (MQ=255)
        ttCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAACGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaa                                                                                                                                           >  3:615157/1‑149 (MQ=255)
         tCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                                                                          <  4:251370/149‑2 (MQ=255)
         tCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCaa                                                                                                                                           >  2:184187/1‑148 (MQ=255)
         tCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAAt                                                                                                                                          <  4:641577/149‑1 (MQ=255)
          cGAAGCCAGTACACGCTGATGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaagt                                                                                                                                         <  4:615157/149‑3 (MQ=255)
            aaGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTc                                                                                                                                       >  4:645987/1‑149 (MQ=255)
              gCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCg                                                                                                                                      <  1:662922/148‑1 (MQ=255)
                 aGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtttt                                                                                                                                  <  3:627844/149‑1 (MQ=255)
                   tACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGttttt                                                                                                                                 <  2:445821/148‑1 (MQ=255)
                        gCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcg                                                                                                                           >  4:660904/1‑149 (MQ=255)
                        gCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTcg                                                                                                                           >  4:73517/1‑149 (MQ=255)
                             aGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                                                                          <  2:76342/129‑2 (MQ=255)
                             aGATCTTCACTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGGATTCGGCaag                                                                                                                                          >  1:76342/1‑128 (MQ=255)
                                  ttCAGTCTTGAGATAATGAGCCAGCACTTCCGCCAGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCCAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACg                                                                                                                 >  3:601698/1‑149 (MQ=255)
                                    cAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGa                                                                                                                >  1:691585/1‑148 (MQ=255)
                                     aGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGAc                                                                                                               >  1:361471/1‑148 (MQ=255)
                                     aGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGAc                                                                                                               <  1:184187/148‑1 (MQ=255)
                                        cTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGt                                                                                                            >  2:563781/1‑148 (MQ=255)
                                            agaTAATGAGCCAGCACTTCCGCCTGTGGCTGTTCGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGCCGTTTTTACCTCGCGCAGTAACGACCGTGAGGc                                                                                                       >  4:794998/1‑149 (MQ=255)
                                             gaTAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCTGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                                                                        <  4:523941/147‑1 (MQ=255)
                                             gaTAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAgg                                                                                                        >  3:523941/1‑147 (MQ=255)
                                                aaTGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAat                                                                                                    <  2:361471/148‑1 (MQ=255)
                                                 aTGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatc                                                                                                   <  2:438086/148‑1 (MQ=255)
                                                   gAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatcat                                                                                                 >  1:667295/1‑148 (MQ=255)
                                                    aGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAatcatc                                                                                                >  2:413394/1‑148 (MQ=255)
                                                     gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAAtcatca                                                                                               <  1:147924/148‑1 (MQ=255)
                                                     gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAAtcatca                                                                                               <  1:738137/148‑1 (MQ=255)
                                                     gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCaa                                                                                              <  3:7661/149‑1 (MQ=255)
                                                     gCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCaa                                                                                              >  4:112638/1‑149 (MQ=255)
                                                       cAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACt                                                                                            >  3:477016/1‑149 (MQ=255)
                                                        aGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGt                                                                                                                                     <  4:376780/107‑1 (MQ=255)
                                                        aGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGt                                                                                                                                     >  3:376780/1‑107 (MQ=255)
                                                        aGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACt                                                                                            >  1:40967/1‑148 (MQ=255)
                                                             ttCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGccccc                                                                                      <  4:638817/149‑1 (MQ=255)
                                                               ccGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtt                                                                                                                                    >  4:645718/1‑101 (MQ=255)
                                                               ccGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGtt                                                                                                                                    <  3:645718/101‑1 (MQ=255)
                                                                 gCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAg                                                                                  <  4:613789/149‑1 (MQ=255)
                                                                     gtgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                                                                              <  4:397532/149‑1 (MQ=255)
                                                                      tgGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGtttt                                                                              <  2:244137/148‑1 (MQ=255)
                                                                                cGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCaa                                                                    <  2:40967/148‑1 (MQ=255)
                                                                                  tAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACg                                                                  >  2:749313/1‑148 (MQ=255)
                                                                                    aTGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTa                                                                >  2:14623/1‑148 (MQ=255)
                                                                                       tCATGCGTAAACAGTGGTCGAAAGCGATGTGCAGCCGGGCAAACGCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGGTTTCGCGCTTCAACGTAACgg                                                            >  4:853115/1‑149 (MQ=255)
                                                                                       tCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACg                                                             <  2:125143/148‑1 (MQ=255)
                                                                                        cATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACggg                                                           <  3:73517/149‑1 (MQ=255)
                                                                                           gCGTAAACAGCGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCg                                                         >  2:287406/1‑148 (MQ=255)
                                                                                            cGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGt                                                        >  2:539439/1‑148 (MQ=255)
                                                                                               aaaCAGCGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTAcc                                                     <  1:287406/148‑1 (MQ=255)
                                                                                                 aCAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCg                                                   <  1:539439/148‑1 (MQ=255)
                                                                                                          gAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAATGACCGTCAGGCAATCATCAACAGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCCGCATTATTAc                                         <  4:586491/149‑1 (MQ=255)
                                                                                                               cGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCattatt                                           >  1:438592/1‑142 (MQ=255)
                                                                                                               cGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCattatt                                           <  2:438592/142‑1 (MQ=255)
                                                                                                               agATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGaaa                                     >  2:725268/2‑148 (MQ=255)
                                                                                                                gATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTCACCTCGCGCAGTAACGACCGTGAGGCAATCATCACCTGCCCCCAAAGTTTTCTCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAAc                                   <  4:601698/149‑1 (MQ=255)
                                                                                                                   ttGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATCACGAAAACCaa                                <  3:794998/149‑1 (MQ=255)
                                                                                                                         ccGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGaaaaa                          >  4:678812/1‑149 (MQ=255)
                                                                                                                          cGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGaaaaa                          <  2:687469/148‑1 (MQ=255)
                                                                                                                          cGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAat                         >  3:68412/1‑149 (MQ=255)
                                                                                                                           gggCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatg                        <  4:477016/149‑1 (MQ=255)
                                                                                                                           gggCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAat                         <  2:68956/148‑1 (MQ=255)
                                                                                                                            ggCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAat                         <  3:620029/147‑1 (MQ=255)
                                                                                                                            ggCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAat                         >  4:620029/1‑147 (MQ=255)
                                                                                                                               aaaCTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcg                     <  2:691585/148‑1 (MQ=255)
                                                                                                                                   tCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatggc                >  3:211821/1‑149 (MQ=255)
                                                                                                                                    cccTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCattat                                            <  2:446131/120‑1 (MQ=255)
                                                                                                                                    cccTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCattat                                            >  1:446131/1‑120 (MQ=255)
                                                                                                                                    cccTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatggc                >  2:44182/1‑148 (MQ=255)
                                                                                                                                      cTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGac             >  4:886204/1‑149 (MQ=255)
                                                                                                                                       ttACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAAc                                   >  4:451315/1‑126 (MQ=255)
                                                                                                                                       ttACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAAc                                   <  3:451315/126‑1 (MQ=255)
                                                                                                                                        tACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGacac           <  4:211559/149‑1 (MQ=255)
                                                                                                                                        tACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGaca            >  1:820300/1‑148 (MQ=255)
                                                                                                                                         gacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatggc                <  4:444794/139‑1 (MQ=255)
                                                                                                                                         gacattGCCGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatggc                >  3:444794/5‑143 (MQ=255)
                                                                                                                                            cTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCga        >  1:110353/1‑148 (MQ=255)
                                                                                                                                             ttGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCgaga      >  4:613838/1‑149 (MQ=255)
                                                                                                                                                ccGGATTCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAAt    >  2:736320/1‑148 (MQ=255)
                                                                                                                                                ccGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatg                  >  2:155607/1‑134 (MQ=255)
                                                                                                                                                ccGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAatggcgatg                  <  1:155607/134‑1 (MQ=255)
                                                                                                                                                 cGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGc  <  3:660904/149‑1 (MQ=255)
                                                                                                                                                  gAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGc  <  1:413394/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAACTCCCTTACACTTGCCGAATCCGGCAATGTCGTTTTTACCTCGCGCAGTAACGACCGTGAGGCAATCATCAACTGCCCCCAAAGTTTTCGCGCTTCAACGTAACGGGCGTACCCGGCATTATTACGAAAACCAAGAAAAATGGCGATGGCGACACCGAGAATGC  >  NC_000913/1608557‑1608850

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