Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,057,226 G→T 10.0% E839* (GAA→TAA)  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
*CP0092731,057,2260GT10.0% 488.3 / 22.8 140E839* (GAA→TAA) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (52/74);  new base T (3/11);  total (55/85)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 140 total aligned reads displayed.

CAACGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAAT  >  CP009273/1057081‑1057372
                                                                                                                                                 |                                                                                                                                                  
cAACGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGaa                                                                                                                                                  >  1:952986/1‑148 (MQ=255)
cAACGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGaa                                                                                                                                                  >  1:1489494/1‑148 (MQ=255)
 accGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAat                                                                                                                                                 <  1:1380957/146‑1 (MQ=255)
 aaCGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAat                                                                                                                                                 <  1:1380916/148‑1 (MQ=255)
 aaCGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAat                                                                                                                                                 >  2:573840/1‑148 (MQ=255)
 aaCGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAACTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAat                                                                                                                                                 <  1:244376/148‑1 (MQ=255)
   cGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAatat                                                                                                                                               >  1:450754/1‑148 (MQ=255)
   cGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAatat                                                                                                                                               >  1:3140641/1‑148 (MQ=255)
   cGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAatat                                                                                                                                               >  1:960733/1‑148 (MQ=255)
     tgtTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGt                                                                                                                                             <  1:1093878/148‑1 (MQ=255)
       ttGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTc                                                                                                                                           <  2:3318302/148‑1 (MQ=255)
        tGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTcc                                                                                                                                          <  1:2548485/148‑1 (MQ=255)
        tGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTcc                                                                                                                                          >  1:1709991/1‑148 (MQ=255)
        tGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTcc                                                                                                                                          >  1:1709987/1‑148 (MQ=255)
            cATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCcgcg                                                                                                                                      <  2:3569381/148‑1 (MQ=255)
             aTCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGt                                                                                                                                     <  2:3488112/148‑1 (MQ=255)
                   aTCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAgg                                                                                                                               <  1:371461/148‑1 (MQ=255)
                   aTCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAgg                                                                                                                               <  1:3416630/148‑1 (MQ=255)
                   aTCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAgg                                                                                                                               >  1:571152/1‑148 (MQ=255)
                   aTCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAgg                                                                                                                               <  1:371924/148‑1 (MQ=255)
                        tACATCGCAGCTGGCGCAGGCGACCAGTGCGCACAATCCGCTGGTGTAAATTGTGAAGTACAACGGAACAGTGTAGCAGGTGACGGCGGTTAAAGGACCAGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGaaa                                                                                                                          >  1:334759/1‑148 (MQ=255)
                         aCATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAAt                                                                                                                         >  1:2210241/1‑148 (MQ=255)
                          cATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATc                                                                                                                        <  2:407101/148‑1 (MQ=255)
                           aTCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCa                                                                                                                       <  1:1367436/148‑1 (MQ=255)
                           aTCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCa                                                                                                                       <  2:2982326/148‑1 (MQ=255)
                           aTCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCa                                                                                                                       <  2:2982331/148‑1 (MQ=255)
                              gCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGa                                                                                                                    >  2:806148/1‑148 (MQ=255)
                                 gCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCa                                                                                                                 <  1:2023369/148‑1 (MQ=255)
                                 gCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCa                                                                                                                 <  2:3005560/148‑1 (MQ=255)
                                    ggCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGc                                                                                                              >  1:1504863/1‑148 (MQ=255)
                                     gcgcAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCt                                                                                                             >  1:3530181/1‑148 (MQ=255)
                                     gcgcAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCt                                                                                                             >  1:3406852/1‑148 (MQ=255)
                                      cgcAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTg                                                                                                            <  2:1010766/148‑1 (MQ=255)
                                      cgcAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTg                                                                                                            >  2:2273637/1‑148 (MQ=255)
                                        cAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGAc                                                                                                          <  1:2628476/148‑1 (MQ=255)
                                              aCCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGcaca                                                                                                    >  1:2829536/1‑148 (MQ=255)
                                              aCCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGcaca                                                                                                    >  1:544643/1‑148 (MQ=255)
                                               ccAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGTTGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCacaa                                                                                                   <  2:1117066/148‑1 (MQ=255)
                                               ccAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCacaa                                                                                                   <  2:963165/148‑1 (MQ=255)
                                                 aGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCacaaca                                                                                                 <  1:2320728/148‑1 (MQ=255)
                                                 aGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCacaaca                                                                                                 <  1:950021/148‑1 (MQ=255)
                                                  gTGCGCACACTACGCTGGTGGAAATTGAGAAGTACACCGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAg                                                                                                <  1:1025390/148‑1 (MQ=255)
                                                         cacTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCt                                                                                         >  2:1675511/1‑148 (MQ=255)
                                                          acTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCtg                                                                                        <  1:3306235/148‑1 (MQ=255)
                                                           cTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCtgt                                                                                       >  1:1060519/1‑148 (MQ=255)
                                                           cTACGCTGGTGGAAAATGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGATTATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCtgt                                                                                       <  2:2202976/148‑1 (MQ=255)
                                                            tACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCtgtg                                                                                      <  1:3153747/148‑1 (MQ=255)
                                                            tACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTAGCAGGTGAAATCATGAACACGCTGACAGCACAACAGATTGCCtgtg                                                                                      <  1:3154668/148‑1 (MQ=255)
                                                             aCGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTgtgt                                                                                     <  2:3141943/148‑1 (MQ=255)
                                                              cGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGtt                                                                                    <  2:1718638/148‑1 (MQ=255)
                                                                   gtggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGc                                                                               <  1:3212280/148‑1 (MQ=255)
                                                                   gtggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGc                                                                               <  2:3683862/148‑1 (MQ=255)
                                                                   gtggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGc                                                                               <  1:1656362/148‑1 (MQ=255)
                                                                    tggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGcc                                                                              <  1:529974/148‑1 (MQ=255)
                                                                    tggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGcc                                                                              <  2:3018083/148‑1 (MQ=255)
                                                                     ggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCt                                                                             >  1:2283075/1‑148 (MQ=255)
                                                                      gAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTg                                                                            >  1:1627963/1‑148 (MQ=255)
                                                                      gAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTg                                                                            >  1:2746787/1‑148 (MQ=255)
                                                                           tGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAg                                                                       >  2:2043422/1‑148 (MQ=255)
                                                                             agaAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTAAGCCTGGCTAgcg                                                                     >  2:2652553/1‑148 (MQ=255)
                                                                              gaAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAgcgc                                                                    <  2:453180/148‑1 (MQ=255)
                                                                                   aCAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAgttg                                                               >  1:66607/1‑148 (MQ=255)
                                                                                    cAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAgttgt                                                              <  2:317449/148‑1 (MQ=255)
                                                                                    cAACGGAACAGTGGAGCAGGTGACGGCGTTTAACAGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAgttgt                                                              <  2:2237966/148‑1 (MQ=255)
                                                                                      aCGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTtc                                                            >  1:339690/1‑148 (MQ=255)
                                                                                       cGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTtct                                                           >  1:951694/1‑148 (MQ=255)
                                                                                       cGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTtct                                                           <  1:244851/148‑1 (MQ=255)
                                                                                        ggAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCATTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTtctc                                                          >  1:766823/1‑148 (MQ=255)
                                                                                        ggAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTtctc                                                          <  1:1579645/148‑1 (MQ=255)
                                                                                            cAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGt                                                      <  1:26177/148‑1 (MQ=255)
                                                                                             aGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTg                                                     <  1:1314693/148‑1 (MQ=255)
                                                                                             aGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTg                                                     <  1:1314622/148‑1 (MQ=255)
                                                                                             aGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTg                                                     <  1:513852/148‑1 (MQ=255)
                                                                                                ggAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGc                                                  <  1:2332858/148‑1 (MQ=255)
                                                                                                ggAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGc                                                  <  1:2341229/148‑1 (MQ=255)
                                                                                                 gAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCt                                                 >  1:946515/1‑148 (MQ=255)
                                                                                                   gCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTgg                                               >  1:376994/1‑148 (MQ=255)
                                                                                                    cAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGa                                              >  1:2945456/1‑148 (MQ=255)
                                                                                                      ggTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACg                                            >  2:2024302/1‑148 (MQ=255)
                                                                                                      ggTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGGTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACg                                            <  1:1964950/148‑1 (MQ=255)
                                                                                                       gTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGa                                           >  2:3501/1‑148 (MQ=255)
                                                                                                              cGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGaacaa                                    <  2:2610503/148‑1 (MQ=255)
                                                                                                                tttAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACt                                  >  2:3109392/1‑148 (MQ=255)
                                                                                                                 ttAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTg                                 <  2:1489494/148‑1 (MQ=255)
                                                                                                                 ttAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTg                                 <  1:1675511/148‑1 (MQ=255)
                                                                                                                    aCGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACg                              >  2:470890/1‑148 (MQ=255)
                                                                                                                     cGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGc                             >  1:3032705/1‑148 (MQ=255)
                                                                                                                       gCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCaa                           <  2:3656705/148‑1 (MQ=255)
                                                                                                                       gCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCaa                           <  2:3803683/148‑1 (MQ=255)
                                                                                                                        ccccGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCaaa                          <  1:259129/148‑1 (MQ=255)
                                                                                                                         cccGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAAt                         <  1:440719/148‑1 (MQ=255)
                                                                                                                                 gaTGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTg                 >  1:2023555/1‑148 (MQ=255)
                                                                                                                                   tggtggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTgcg               >  2:882088/1‑148 (MQ=255)
                                                                                                                                      tggCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAg            <  1:1579907/148‑1 (MQ=255)
                                                                                                                                           cAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGc       <  2:3606332/148‑1 (MQ=255)
                                                                                                                                           cAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGc       <  2:23003/148‑1 (MQ=255)
                                                                                                                                              tGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGa    <  2:3662306/148‑1 (MQ=255)
                                                                                                                                              tGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGa    <  2:141305/148‑1 (MQ=255)
                                                                                                                                                cGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAAt  <  2:544643/148‑1 (MQ=255)
                                                                                                                                                cGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAAt  <  2:2829536/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAACGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAAT  >  CP009273/1057081‑1057372

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