Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,222,778 A→G 5.3% N49S (AAT→AGT)  ebgA → evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,222,7780AG 5.3% 239.2 / 3.2 76N49S (AAT→AGT) ebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (36/36);  new base G (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

CCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTC  >  NC_000913/3222630‑3222923
                                                                                                                                                    |                                                                                                                                                 
ccTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGaa                                                                                                                                                   >  5:186917/1‑149 (MQ=255)
ccTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGaa                                                                                                                                                   >  2:377/1‑149 (MQ=255)
 cTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAAt                                                                                                                                                  >  6:52989/1‑149 (MQ=255)
      aaTCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCa                                                                                                                                             <  5:88484/149‑1 (MQ=255)
          gCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACtt                                                                                                                                          <  5:52989/148‑1 (MQ=255)
             gggAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACtttttt                                                                                                                                      >  8:279523/1‑149 (MQ=255)
              ggAAAACATTCATCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCTCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGTAATTTCCACtttttt                                                                                                                                      <  1:348929/148‑1 (MQ=255)
              ggAAAACATTCAGCGCACCCCCGTAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGCTGCGCAAGCGCGTTCCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTACGCGGGCAGTGGAATTTCCACttttt                                                                                                                                       <  1:373506/147‑1 (MQ=255)
                aaaaCATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGAc                                                                                                                                   >  8:36228/1‑149 (MQ=255)
                aaaaCATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGAc                                                                                                                                   <  1:377/149‑1 (MQ=255)
                    cATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGATGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATc                                                                                                                               <  1:118195/149‑1 (MQ=255)
                      ttCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAGTTTCCACTTTTTTGACCATCCg                                                                                                                             >  2:152996/1‑149 (MQ=255)
                         aGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTg                                                                                                                          >  1:172826/1‑149 (MQ=255)
                          gCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAGTTTCCACTTTTTTGACCATCCGCTGc                                                                                                                         <  1:152996/149‑1 (MQ=255)
                              aCCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAg                                                                                                                      >  5:158461/1‑148 (MQ=255)
                                ccACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTAc                                                                                                                   >  6:277020/1‑149 (MQ=255)
                                    gAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGa                                                                                                                                    >  6:101999/1‑128 (MQ=255)
                                    gAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGa                                                                                                                                    <  5:101999/128‑1 (MQ=255)
                                    gAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAg                                                                                                                >  7:23659/1‑148 (MQ=255)
                                            aCTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCa                                                                                                       <  2:172826/149‑1 (MQ=255)
                                                 cgccgcGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACctc                                                                                                   <  7:279523/148‑1 (MQ=255)
                                                  gccgcGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACctct                                                                                                  >  1:68902/1‑148 (MQ=255)
                                                     gcgtgcgtACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGt                                                                                              <  7:337125/149‑1 (MQ=255)
                                                         gcgtACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACctct                                                                                                  >  2:263382/1‑141 (MQ=255)
                                                         gcgtACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACctct                                                                                                  <  1:263382/141‑1 (MQ=255)
                                                          cgtACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATg                                                                                         <  2:11333/149‑1 (MQ=255)
                                                             aCTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGCCCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCt                                                                                      >  8:349507/1‑149 (MQ=255)
                                                             aCTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGa                                                                                                >  5:30736/1‑139 (MQ=255)
                                                             aCTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGa                                                                                                <  6:30736/139‑1 (MQ=255)
                                                               tttttttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGa                                                                                    >  8:165328/1‑149 (MQ=255)
                                                               tttttttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGa                                                                                    >  8:323184/1‑149 (MQ=255)
                                                                 tttttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGa                                                                                    <  7:36228/147‑1 (MQ=255)
                                                                   tttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTg                                                                                 <  2:68902/148‑1 (MQ=255)
                                                                    ttCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTggg                                                                               >  3:266588/1‑149 (MQ=255)
                                                                     tCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTggg                                                                               <  7:165328/148‑1 (MQ=255)
                                                                         atGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCaa                                                                                                                       <  2:282608/104‑1 (MQ=255)
                                                                         atGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCaa                                                                                                                       >  1:282608/1‑104 (MQ=255)
                                                                           gATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCat                                                                          <  1:308426/147‑1 (MQ=255)
                                                                           gATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCat                                                                          >  2:308426/1‑147 (MQ=255)
                                                                                tGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCg                                                                   >  2:241764/1‑149 (MQ=255)
                                                                                  ttGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGt                                                                  <  6:186917/148‑1 (MQ=255)
                                                                                  ttGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTc                                                                 >  1:80085/1‑149 (MQ=255)
                                                                                            gcgTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAGTTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGa                                                                                                >  4:228636/1‑108 (MQ=255)
                                                                                            gcgTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAGTTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGa                                                                                                <  3:228636/108‑1 (MQ=255)
                                                                                              gTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATg                                                                                         >  3:288536/1‑113 (MQ=255)
                                                                                              gTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATg                                                                                         <  4:288536/113‑1 (MQ=255)
                                                                                                  cTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACTATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCaaa                                                 >  2:105819/1‑149 (MQ=255)
                                                                                                  cTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTggg                                                                               <  7:53170/119‑1 (MQ=255)
                                                                                                  cTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTggg                                                                               >  8:53170/1‑119 (MQ=255)
                                                                                                     tGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATg                                               <  2:80085/148‑1 (MQ=255)
                                                                                                       cccGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTTACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGaa                                            <  6:158461/149‑1 (MQ=255)
                                                                                                       cccGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATgg                                              >  1:299908/1‑147 (MQ=255)
                                                                                                       cccGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGaa                                            <  8:23659/149‑1 (MQ=255)
                                                                                                       cccGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGaa                                            <  8:16058/149‑1 (MQ=255)
                                                                                                        ccGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAg                                           >  4:149579/1‑149 (MQ=255)
                                                                                                          gcgAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGt                                         >  4:13410/1‑149 (MQ=255)
                                                                                                             aaaCCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACTATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCAc                                      <  1:105819/149‑1 (MQ=255)
                                                                                                               aCCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACgg                                    <  7:323184/149‑1 (MQ=255)
                                                                                                                       ccTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATgg                                                                                        >  4:39274/1‑89 (MQ=255)
                                                                                                                       ccTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATgg                                                                                        <  3:39274/89‑1 (MQ=255)
                                                                                                                        cTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCAGCCATGTGGCAAATGGAAGGTCACGGCAAACTGc                            >  3:25555/1‑148 (MQ=255)
                                                                                                                           tttCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTggggg                                                                             <  8:273016/96‑1 (MQ=255)
                                                                                                                           tttCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTggggg                                                                             >  7:273016/1‑96 (MQ=255)
                                                                                                                            ttCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAata                        >  2:89707/1‑148 (MQ=255)
                                                                                                                              cTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAAtata                      <  1:241764/148‑1 (MQ=255)
                                                                                                                                 cccTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACcga                  >  8:74804/1‑149 (MQ=255)
                                                                                                                                  ccTTAAGCTGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACcgac                 >  2:194687/1‑149 (MQ=255)
                                                                                                                                   cTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGaa                                            >  7:349107/1‑121 (MQ=255)
                                                                                                                                   cTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGaa                                            <  8:349107/121‑1 (MQ=255)
                                                                                                                                        gCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCTTCCCCGCCATGTGGAAAATGGAAGGACACGGCAAACTGCGAAAAACCGACGATGGt           >  3:53485/1‑149 (MQ=255)
                                                                                                                                            tCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTc       <  3:13410/149‑1 (MQ=255)
                                                                                                                                               gTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCaaa                                <  2:39548/121‑1 (MQ=255)
                                                                                                                                               gTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCaaa                                >  1:39548/1‑121 (MQ=255)
                                                                                                                                                tGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAat                         <  2:231432/127‑1 (MQ=255)
                                                                                                                                                tGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAat                         >  1:231432/1‑127 (MQ=255)
                                                                                                                                                tGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACcg                   <  8:160196/133‑1 (MQ=255)
                                                                                                                                                tGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACcg                   >  7:160196/1‑133 (MQ=255)
                                                                                                                                                  gAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAATTGCAATATACCGACGAAGGTTTTCCGTTc  <  3:149579/148‑1 (MQ=255)
                                                                                                                                                    aTTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCaa                          <  6:74162/122‑1 (MQ=255)
                                                                                                                                                    aTTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCaa                          >  5:74162/1‑122 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCTATGAATCGCTGGGAAAACATTCAGCTCACCCACGAAAACCGACTTGCGCCGCGTGCGTACTTTTTTTCATATGATTCTGTTGCGCAAGCGCGTACCTTTGCCCGCGAAACCAGCAGCCTGTTTCTGCCCTTAAGCGGTCAGTGGAATTTCCACTTTTTTGACCATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTC  >  NC_000913/3222630‑3222923

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