Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,061,193 G→A 5.5% T646T (ACG→ACA)  torA → trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,061,1930GA 5.5% 231.2 / 2.4 72T646T (ACG→ACA) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base A (2/2);  total (34/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 = 4.86e-01

TGTACAGCAAGGCAAAGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGT  >  NC_000913/1061046‑1061339
                                                                                                                                                   |                                                                                                                                                  
tgtACAGCAAGGCAAAGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGc                                                                                                                                                   <  5:142221/149‑1 (MQ=255)
    cAGCAAGGCAAAGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGa                                                                                                                                                <  8:383104/148‑1 (MQ=255)
     aGCAAGGCAAAGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACACCGAGt                                                                                                                                              >  4:168851/1‑149 (MQ=255)
     aGCAAGGCAAAGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACACCGAGt                                                                                                                                              >  4:163697/1‑149 (MQ=255)
       cacggCAAAGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTgg                                                                                                                                            <  6:122328/146‑1 (MQ=255)
       cAAGGCAAAGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTgg                                                                                                                                            >  1:179085/1‑149 (MQ=255)
       cAAGGCAAAGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTgg                                                                                                                                            >  5:309598/1‑149 (MQ=255)
          ggCAAAGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCt                                                                                                                                         >  2:231484/1‑149 (MQ=255)
           gCAAAGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCt                                                                                                                                         >  7:190079/1‑148 (MQ=255)
             aaaGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTg                                                                                                                                             <  1:94867/142‑1 (MQ=255)
             aaaGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTg                                                                                                                                             >  2:94867/1‑142 (MQ=255)
                  aCGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGAtt                                                                                                                                     >  6:292245/1‑145 (MQ=255)
                  aCGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGAtt                                                                                                                                     <  5:292245/145‑1 (MQ=255)
                  aCGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTgaga                                                                                                                                 <  7:382126/149‑1 (MQ=255)
                    gcggcgTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTga                                                                                                                                   >  4:321280/1‑145 (MQ=255)
                    gcggcgTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTga                                                                                                                                   <  3:321280/145‑1 (MQ=255)
                         gTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCCTCCGCACAGGTTTGGTCGCCACCAGGCATTCCCCGAAGATCCGGATCTCGAACCGCAGGGCACGCCGAGTGGCCTGATTTAGATCTACTc                                                                                                                          >  4:32564/1‑149 (MQ=255)
                         gTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTc                                                                                                                          >  7:169608/1‑149 (MQ=255)
                             aTCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGaa                                                                                                                       <  5:206648/148‑1 (MQ=255)
                                tGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAAc                                                                                                                    >  4:131697/1‑148 (MQ=255)
                                  ccAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACt                                                                                                                           >  4:178492/1‑139 (MQ=255)
                                  ccAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACt                                                                                                                           <  3:178492/139‑1 (MQ=255)
                                  ccAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTAt                                                                                                                 >  4:134270/1‑149 (MQ=255)
                                   cAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTAc                                                                                                                            >  2:184008/1‑137 (MQ=255)
                                   cAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTAc                                                                                                                            <  1:184008/137‑1 (MQ=255)
                                   cAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATc                                                                                                                >  2:229361/1‑149 (MQ=255)
                                   cAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATc                                                                                                                <  6:129033/149‑1 (MQ=255)
                                    aGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTg                                                                                                                                        >  8:289251/1‑124 (MQ=255)
                                    aGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTg                                                                                                                                        <  7:289251/124‑1 (MQ=255)
                                           gatgaCTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTgaga                                                                                                                                 <  2:109496/124‑1 (MQ=255)
                                           gatgaCTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTgaga                                                                                                                                 >  1:109496/1‑124 (MQ=255)
                                           gatgaCTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGata                                                                                                         <  3:134270/148‑1 (MQ=255)
                                                 ttCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTa                                                                                                  >  1:5221/1‑149 (MQ=255)
                                                 ttCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTa                                                                                                  >  2:382610/1‑149 (MQ=255)
                                                 ttCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTa                                                                                                  >  6:87337/1‑149 (MQ=255)
                                                   cTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTac                                                                                                 <  3:131697/148‑1 (MQ=255)
                                                         tAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGAtt                                                                                          <  5:87337/149‑1 (MQ=255)
                                                         tAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGAtt                                                                                          <  3:317460/149‑1 (MQ=255)
                                                           aCAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTg                                                                                         >  2:186160/1‑148 (MQ=255)
                                                            cAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTc                                                                                       <  7:230947/149‑1 (MQ=255)
                                                             aaaGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCTCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCa                                                                                      <  1:186160/149‑1 (MQ=255)
                                                             aaaGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGt                                                                                        >  8:334680/1‑147 (MQ=255)
                                                               agagTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAgg                                                                                    <  3:36399/149‑1 (MQ=255)
                                                                         gAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGa                                                                           <  7:334680/148‑1 (MQ=255)
                                                                           gTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGAtgt                                                                        <  2:179085/149‑1 (MQ=255)
                                                                              tGACCATCCGCAGATGTTTTTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGaaaa                                                                                                                     <  6:345273/101‑1 (MQ=255)
                                                                              tGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGaaaa                                                                                                                     >  5:345273/1‑101 (MQ=255)
                                                                               gACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACACCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGt                                                                     <  3:163697/148‑1 (MQ=255)
                                                                               gACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACACCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGt                                                                     <  3:168851/148‑1 (MQ=255)
                                                                                aCCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGata                                                                                                         <  4:324840/111‑1 (MQ=255)
                                                                                aCCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGata                                                                                                         >  3:324840/1‑111 (MQ=255)
                                                                                aCCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGttt                                                                   <  6:309598/149‑1 (MQ=255)
                                                                                  cATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTga                                                                 <  1:229361/149‑1 (MQ=255)
                                                                                            tgtttgttCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGaa                                                       <  8:169608/149‑1 (MQ=255)
                                                                                            tgtttgttCGCCACCAGGCATTCCGCGAAGATCCGGATCTAGAACCGCTGGGCACGCCGAGTCGCCTGATTTAGATCTACTCTAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGaa                                                       <  5:121046/149‑1 (MQ=255)
                                                                                                 gttCGCCACAAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCTAACGCTc                                                  >  6:108200/1‑149 (MQ=255)
                                                                                                     gCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCAc                                              >  8:231212/1‑149 (MQ=255)
                                                                                                      ccaccaGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACg                                             >  7:53678/1‑149 (MQ=255)
                                                                                                       caccaGGCATTCCGCGAAGATCCGTATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGttt                                                                   <  4:244225/126‑1 (MQ=255)
                                                                                                       caccaGGCATTCCGCGAAGATCCGTATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGttt                                                                   >  3:244225/1‑126 (MQ=255)
                                                                                                         ccaGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTg                                          <  1:382610/149‑1 (MQ=255)
                                                                                                          caGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCTGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTgg                                         <  1:231484/149‑1 (MQ=255)
                                                                                                                   cgcgAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTc                                <  8:53678/149‑1 (MQ=255)
                                                                                                                   cgcgAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTGTGAGAAAATCGCACGCTCCCACGGTGGGCCTGGc                                  >  1:381093/1‑147 (MQ=255)
                                                                                                                         gATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCaaaa                          <  8:190079/149‑1 (MQ=255)
                                                                                                                         gATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCaaaa                          <  5:108200/149‑1 (MQ=255)
                                                                                                                            ccGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTa                       >  5:167857/1‑149 (MQ=255)
                                                                                                                                    cGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTg                >  4:26481/1‑148 (MQ=255)
                                                                                                                                       aCCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATc            >  3:348691/1‑149 (MQ=255)
                                                                                                                                            tGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCaa       >  8:263327/1‑149 (MQ=255)
                                                                                                                                                 aCGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCtgt  <  6:167857/149‑1 (MQ=255)
                                                                                                                                                 aCGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCtgt  <  2:5221/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGTACAGCAAGGCAAAGGACGCGGCGTTCATCTGCCAGCGTTTGATGACTTCTGGAATAACAAAGAGTACGTCGAGTTTGACCATCCGCAGATGTTTGTTCGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGT  >  NC_000913/1061046‑1061339

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