Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,060,388 T→C 5.5% M378T (ATG→ACG)  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,3880TC 5.5% 208.9 / 2.9 73M378T (ATG→ACG) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (35/34);  new base C (2/2);  total (37/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 = 7.96e-01

TCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTT  >  NC_000913/1060242‑1060534
                                                                                                                                                  |                                                                                                                                                  
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                                                                                                                                     >  1:685998/1‑149 (MQ=255)
             aaaGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACa                                                                                                                                     >  2:163204/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                                                                           >  1:682902/1‑149 (MQ=255)
                                                                       aaaCCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAggt                                                                           <  2:558817/149‑1 (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                                              >  2:437508/1‑116 (MQ=255)
                                                                                                                                     tgcgtgcAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACaa                                              <  1:437508/116‑1 (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)
                                                                                                                                               gTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTGCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGt   <  2:1025057/149‑1 (MQ=255)
                                                                                                                                                tATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGtt  >  1:628024/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGCTGGCACTACAACGGCGCAGGCACGCCGGGGCGTAAAGGCGTTATTCTGAGTGGTT  >  NC_000913/1060242‑1060534

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