Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,060,382 A→T 5.4% Q376L (CAG→CTG)  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,060,3820AT 5.4% 204.3 / 2.4 74Q376L (CAG→CTG) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base T (2/2);  total (38/36)
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.74e-01

TGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCT  >  NC_000913/1060235‑1060526
                                                                                                                                                   |                                                                                                                                                
tGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCac                                                                                                                                                 >  2:910724/1‑148 (MQ=255)
   cGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGTCTGAACAACTGACCGGCAGAGCTGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGCCCGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGt                                                                                                                                              >  1:972484/1‑149 (MQ=255)
    gaaTCTCCTTGGTGCGAAAGACGGGCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTa                                                                                                                                             <  2:671915/147‑1 (MQ=255)
    gTATCAGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCAGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTa                                                                                                                                             >  2:57498/1‑149 (MQ=255)
       tCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGCAAAACTGACCGGCATTGATGCCGAAACCATTCGTTGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATgc                                                                                                                                          <  2:165079/149‑1 (MQ=255)
       tCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATgc                                                                                                                                          <  1:477759/149‑1 (MQ=255)
         tgctgGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATgcag                                                                                                                                        <  1:706637/149‑1 (MQ=255)
           ctgGTTGAGAAAGCCGGTCAGCCTAAAGATGCCGCTTGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGGGGGGGGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCTGCATACgctgca                                                                                                                                      <  2:207285/149‑1 (MQ=255)
             ggttGAGAAAGACGGTCAGCCGAAAGATGCCTCTTTGGCTGAAAACCTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGATCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACg                                                                                                                                    <  2:195740/146‑1 (MQ=255)
             gggTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACg                                                                                                                                    >  2:891193/1‑149 (MQ=255)
              ggTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCTGCATACGCTGCACgc                                                                                                                                   <  2:696431/149‑2 (MQ=255)
               gTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGt                                                                                                                                  <  2:674542/149‑1 (MQ=255)
                tGAGAAAGACGGTCAGCCGAACGATGCCGCATGGGCAGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTg                                                                                                                                 <  2:456791/149‑1 (MQ=255)
                    aaaGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACa                                                                                                                             >  2:163204/1‑149 (MQ=255)
                    aaaGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACa                                                                                                                             >  1:685998/1‑149 (MQ=255)
                       gACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTg                                                                                                                          <  1:1113988/149‑1 (MQ=255)
                            tCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGt                                                                                                                     <  2:236115/149‑1 (MQ=255)
                              aGCCGAAAGATGCCGCATGGTCTGCAACACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGTCAGATGGCGGCGAACAGAACGAAAATTCTTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTgg                                                                                                                   <  1:818286/149‑1 (MQ=255)
                               gCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTggt                                                                                                                  >  1:833004/1‑148 (MQ=255)
                                ccGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGgat                                                                                                                 >  1:146441/1‑149 (MQ=255)
                                ccGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTCGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGgat                                                                                                                 >  1:530716/1‑149 (MQ=255)
                                  gAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGgatga                                                                                                               >  2:156562/1‑149 (MQ=255)
                                   aaaGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGgatgat                                                                                                              >  1:1083334/1‑149 (MQ=255)
                                          ccGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGtt                                                                                                       >  2:401468/1‑149 (MQ=255)
                                            gCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCt                                                                                                     <  1:94905/149‑1 (MQ=255)
                                                    tGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGa                                                                                             >  1:922329/1‑149 (MQ=255)
                                                      aaaaaCTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATg                                                                                           >  1:182263/1‑149 (MQ=255)
                                                        caaCTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGGGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGAGGATTGTGGTTCTGGCGGCGATGCt                                                                                         <  2:1117009/148‑1 (MQ=255)
                                                           cTGACCGGCATTGATGCCGAACCCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGCACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTggg                                                                                      >  1:1025057/1‑149 (MQ=255)
                                                            tGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGGGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTgggg                                                                                     <  1:534778/149‑1 (MQ=255)
                                                            tGACCGGCATTGATGCCGAAACCATTCGTCGGCGGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTgggg                                                                                     <  1:240117/149‑1 (MQ=255)
                                                                  gcatcaaTGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGATGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGCCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGCCAAttt                                                                               >  2:521369/7‑149 (MQ=255)
                                                                     ttGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGc                                                                            <  1:626818/149‑1 (MQ=255)
                                                                      tGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGcc                                                                           >  1:182265/1‑149 (MQ=255)
                                                                         tGCCGCAACCAGTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGGGCGTGCAGCGTATGCAGCACGGGGAACAGTGGGCGTGGATGATTGTGGTGCTGGCGGCGATGCTGGGGCAAATTGGCCTGc                                                                        <  1:1078386/149‑1 (MQ=255)
                                                                           ccGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGCACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGCACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGCTGCTGGGGCAAATTGGCCTGCCa                                                                      >  1:65113/1‑149 (MQ=255)
                                                                           ccGAAACCATACGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCa                                                                      >  2:1047346/1‑149 (MQ=255)
                                                                             gAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAgg                                                                    >  1:815000/1‑149 (MQ=255)
                                                                              aaaCCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAggt                                                                   <  2:558817/149‑1 (MQ=255)
                                                                              aaaCCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAggt                                                                   >  1:682902/1‑149 (MQ=255)
                                                                               aaCCATGCGGGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATGATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGACGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAggtg                                                                  <  1:92882/149‑1 (MQ=255)
                                                                                aCCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGAAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGAGTTTGGTTCTGGCGGCGATGCTGGGGGAAATTGTCCTGCCAggtgg                                                                 >  2:642226/1‑149 (MQ=255)
                                                                                 ccATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGGGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAggtggt                                                                <  1:864506/149‑1 (MQ=255)
                                                                                   aTTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAggtggtgg                                                              <  1:401468/149‑1 (MQ=255)
                                                                                        tGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTTGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGtggttttg                                                         >  2:1162769/1‑149 (MQ=255)
                                                                                        tGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCTGCGATGCTGGGGCAAATTTGACTGCAAGGTGGtggttttg                                                         >  2:1003916/1‑149 (MQ=255)
                                                                                        tGGGCTGCCGCGGCAGATGGCGGCGAACAGAACGCAAATAATTGCTGGCTGGAGCGTACAGCGTATGCAGCACGGTGAACAGAGGGCGTGGATGATTGTGGATCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGCTGGtggatttg                                                         >  2:342081/1‑149 (MQ=255)
                                                                                        ggggCTGGCGCGGCAGATTGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATGGTGGTTCTGGCGGCGATGCGGGGGCAAATTGGCCTGCCAGGTGGtggttttg                                                         <  1:384435/148‑1 (MQ=255)
                                                                                              ggCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTCCGGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCGGGGGCAAATTGGCCTGCCAGGTGGTGgttttggttttg                                                   <  2:923990/149‑1 (MQ=255)
                                                                                              ggCGCGGCAGATGGCGGAGAACAGAACGCAAATTATTGCTGGCTGGTGCCTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGAGGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGgttttggttttg                                                   <  2:655317/149‑1 (MQ=255)
                                                                                               gcgcGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGttttggttttgg                                                  >  1:900877/1‑149 (MQ=255)
                                                                                                    gCAGATGGCGGCGAACAGAACGCAAATTAATGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCGGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTtggctggc                                             <  2:99274/149‑1 (MQ=255)
                                                                                                      aGATGGCTGAGACCAGCACGCAAATTTCTGCTGGCTGGTGCGTGCAGCGTCTGCAGCAGGGTGAACAGTGGGCGTGGATGATTGTGTTTCTCGCGGCGATGCTGGGGCATATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCAc                                           <  2:805154/149‑1 (MQ=255)
                                                                                                        aTGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTa                                         >  2:478804/1‑149 (MQ=255)
                                                                                                              gcgAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACgg                                   >  1:968707/1‑149 (MQ=255)
                                                                                                              gcgAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACgg                                   >  1:791492/1‑149 (MQ=255)
                                                                                                                    aGACCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCGGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAgg                             <  1:1053595/149‑1 (MQ=255)
                                                                                                                     gAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGc                            <  2:682902/149‑1 (MQ=255)
                                                                                                                      aaCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCa                           >  2:460918/1‑149 (MQ=255)
                                                                                                                       aCGCAAATTATTGCTGGCTGGGGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCAc                          <  1:982710/149‑1 (MQ=255)
                                                                                                                             attattGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCggg                    >  2:609382/1‑149 (MQ=255)
                                                                                                                              ttattGCTGACTGGTGCGTGAAGCGTATGCAGCACGGTGAACAGTTGGCGTGGATGATTGTGGTTCTGGCGGCCATCCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTCGTTTTGGCTGGCACTACAACGGCGCAGGCACGACgggg                   >  1:71863/1‑149 (MQ=255)
                                                                                                                                 ttGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGt                >  1:208919/1‑149 (MQ=255)
                                                                                                                                    ctggctggTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGGGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTaaa             >  1:1099189/1‑149 (MQ=255)
                                                                                                                                     gggctggTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGGTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGGGGTGGTGTTGGTTTTGGCGGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAg            <  2:556916/148‑1 (MQ=255)
                                                                                                                                        gttcaccgtgcTGCATACGCTGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCg         >  1:1090081/7‑149 (MQ=255)
                                                                                                                                        ctggTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGACACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCg         <  1:163204/149‑1 (MQ=255)
                                                                                                                                          ggcgcgggcAGCGGATGCAGCACGGTGAACAGTGGGCGTGGATGATAGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCGGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGtt       <  1:568535/146‑1 (MQ=255)
                                                                                                                                           gtgcgtgcAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTa      >  2:632258/1‑149 (MQ=255)
                                                                                                                                            tgcgtgcAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACaa                                      <  1:437508/116‑1 (MQ=255)
                                                                                                                                            tgcgtgcAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACaa                                      >  2:437508/1‑116 (MQ=255)
                                                                                                                                             gcgcgcAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGCTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTAtt    <  1:123240/149‑1 (MQ=255)
                                                                                                                                              cgtgcTGCATACGCTGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTc   >  1:497078/1‑149 (MQ=255)
                                                                                                                                               gtgcAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCt  <  1:156562/149‑1 (MQ=255)
                                                                                                                                               atgcAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCt  >  1:878676/2‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCT  >  NC_000913/1060235‑1060526

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