Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,060,946 A→G E823G (GAG→GGG)  torA → trimethylamine N‑oxide reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,060,9460AG98.5% 429.1 / ‑2.7 136E823G (GAG→GGG) torAtrimethylamine N‑oxide reductase subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (60/74);  total (61/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

Only 100 of 145 total aligned reads displayed.

TGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGAC  >  NC_000913/1060856‑1061035
                                                                                          |                                                                                         
tGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGc                                                                                         <  2:2497498/93‑1 (MQ=255)
 gACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCa                                                                                        >  1:2087154/1‑93 (MQ=255)
   ccATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAgg                                                                                      <  1:1859080/93‑1 (MQ=255)
    cATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACGCTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGt                                                                                     <  2:238863/93‑1 (MQ=255)
     aTCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTg                                                                                    <  2:4323367/93‑1 (MQ=255)
      tCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGGGa                                                                                   >  1:3503305/1‑93 (MQ=255)
       cGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGAc                                                                                  <  1:929365/93‑1 (MQ=255)
         aCATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACgg                                                                                >  2:310375/1‑93 (MQ=255)
          cATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGc                                                                               >  1:1315717/1‑93 (MQ=255)
            tCGGTACATCGCAGCTGGCGCAGGCGACCCGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGt                                                                             >  2:623361/1‑93 (MQ=255)
            tCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACCCTACGCTGGTGGAAATTGAGAAGTACAACGGCACAGTGGGGCAGGGGACGGGGt                                                                             >  2:4536863/1‑93 (MQ=255)
            tCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGt                                                                             >  1:4652756/1‑93 (MQ=255)
            tCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGt                                                                             >  1:1530279/1‑93 (MQ=255)
              ggTGCATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTGCGCTGGTGGAAATTGAGAAGTACACCGGACCAGTGGGGCAGGTGACGGCGttt                                                                           <  1:3598706/93‑1 (MQ=255)
              ggTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGttt                                                                           <  1:3509853/93‑1 (MQ=255)
                tACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGCGGGGCAGGTGACGGGGTTTaa                                                                         >  2:3863826/1‑93 (MQ=255)
                   aTCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAGCAGTGGGGCAGGTGACGGCGTTTAACgg                                                                      <  1:4159678/93‑1 (MQ=255)
                    tCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGc                                                                     <  2:4448545/93‑1 (MQ=255)
                     cGCAGCTGAAGCAGGAGACCAGTGCGCACACAACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGAAGGTGACGGCGTTTAACGGcc                                                                    <  2:4350314/93‑1 (MQ=255)
                         gCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGt                                                                <  2:1866989/93‑1 (MQ=255)
                          cTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTTACGGCCCCGTg                                                               >  2:1936104/1‑93 (MQ=255)
                              cgcAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGgaga                                                           <  1:669127/93‑1 (MQ=255)
                              cgcAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGgaga                                                           <  1:1546526/93‑1 (MQ=255)
                               gcAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGAt                                                          >  2:1741468/1‑93 (MQ=255)
                                cAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGAtg                                                         <  2:502315/93‑1 (MQ=255)
                                cAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGAtg                                                         >  1:920588/1‑93 (MQ=255)
                                cAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGAtg                                                         <  2:2970599/93‑1 (MQ=255)
                                  ggCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGAtggt                                                       <  1:2034034/93‑1 (MQ=255)
                                   gCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGAtggtg                                                      >  1:4187816/1‑93 (MQ=255)
                                    cGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGAtggtgg                                                     <  1:2613588/93‑1 (MQ=255)
                                     gACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGgc                                                    <  1:3055939/93‑1 (MQ=255)
                                       ccAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTCACGGCCCCGTGGAGATGGTGgcgc                                                  >  2:2743227/1‑93 (MQ=255)
                                         aGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAg                                                >  2:4227448/1‑93 (MQ=255)
                                          gTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAgg                                               >  2:3425337/1‑92 (MQ=255)
                                          gTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGt                                               >  2:2861857/1‑93 (MQ=255)
                                          gTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGt                                               >  2:670934/1‑93 (MQ=255)
                                          gTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACCGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGt                                               <  1:4552821/93‑1 (MQ=255)
                                           tGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTg                                              >  1:3627764/1‑93 (MQ=255)
                                            gcgcACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGc                                             <  1:1698130/93‑1 (MQ=255)
                                            gcgcACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGc                                             >  2:113303/1‑93 (MQ=255)
                                             cgcACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCg                                            <  2:1995156/93‑1 (MQ=255)
                                             cgcACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCg                                            <  2:1995035/93‑1 (MQ=255)
                                             cgcACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCg                                            >  1:1710975/1‑93 (MQ=255)
                                              gcACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGa                                           <  2:3204152/93‑1 (MQ=255)
                                              ccacaCTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGa                                           <  2:3394073/92‑1 (MQ=255)
                                               cacaCTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCGGTGCgga                                          >  2:1198649/1‑91 (MQ=255)
                                               cacaCTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGaa                                          <  1:1259708/93‑1 (MQ=255)
                                               cacaCTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGaa                                          >  1:2603679/1‑93 (MQ=255)
                                               cacaCTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCG‑TTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAat                                         <  1:1129260/93‑1 (MQ=255)
                                                acacTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAat                                         >  2:3197099/1‑93 (MQ=255)
                                                acacTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAat                                         >  2:1676286/1‑93 (MQ=255)
                                                acacTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAat                                         <  2:3344010/93‑1 (MQ=255)
                                                acacTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAat                                         <  2:1778040/93‑1 (MQ=255)
                                                 cacTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAata                                        <  1:2963834/93‑1 (MQ=255)
                                                 cacTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAata                                        <  1:977967/93‑1 (MQ=255)
                                                 cacTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAata                                        <  1:2963833/93‑1 (MQ=255)
                                                 cacTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAata                                        <  2:2461454/93‑1 (MQ=255)
                                                   cTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATg                                      >  2:506328/1‑93 (MQ=255)
                                                    tACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGt                                     >  1:4021970/1‑93 (MQ=255)
                                                    tACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGt                                     >  2:718066/1‑93 (MQ=255)
                                                      cGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTc                                   >  2:2554600/1‑93 (MQ=255)
                                                      cGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTc                                   >  2:1008468/1‑93 (MQ=255)
                                                       gCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTcc                                  <  2:917522/93‑1 (MQ=255)
                                                       gCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTcc                                  <  2:917523/93‑1 (MQ=255)
                                                        cTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTccc                                 <  2:2197775/93‑1 (MQ=255)
                                                        cTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTccc                                 <  1:2871587/93‑1 (MQ=255)
                                                        cTGGTGGAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTccc                                 <  1:2875432/93‑1 (MQ=255)
                                                            tggAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGt                             <  1:1897578/93‑1 (MQ=255)
                                                            tggAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGt                             <  1:3091276/93‑1 (MQ=255)
                                                             ggAAATTGAGAAGTACAACGGAACCGTGGGGCAGGTGACGGCGTTTAACGGCCCCGGGGCGATGGAGGCGCAGTGCGAATATGTTCCCGCGTc                            >  1:4141774/1‑93 (MQ=255)
                                                             ggAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTc                            <  2:2391924/93‑1 (MQ=255)
                                                             ggAAATTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTc                            >  1:4141773/1‑93 (MQ=255)
                                                               aaaTTTAGAAAAACAACAGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCCTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGc                          <  2:1479312/93‑1 (MQ=255)
                                                               aaaGTGAGAAGTACAACGGCACAGTGGGGCAGGTTACGGCGTTTGACGGCCCCGGGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGgcgc                          <  1:1090608/93‑1 (MQ=255)
                                                                aaTTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCa                         <  2:2830181/93‑1 (MQ=255)
                                                                aaTTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCa                         <  1:2562880/93‑1 (MQ=255)
                                                                 aTTGAGAAGTACCACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCCTGGAGGTGGTGGCGCACTGCGTATATGTTCCCGCGCCGCAg                        >  2:1899651/1‑93 (MQ=255)
                                                                 aTTGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAg                        <  1:2293250/93‑1 (MQ=255)
                                                                  ttGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTTGCGCAGCGCGAATATGTTCCCGCGTCGCAgg                       >  2:2216261/1‑93 (MQ=255)
                                                                  atGAGAAGTACAACGGAACAGTGGGGCAGGTGACGTCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAgg                       <  2:3933554/92‑1 (MQ=255)
                                                                   tGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGt                      <  2:4263103/93‑1 (MQ=255)
                                                                   tGAGAAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGACCCGTGGAGATGGTGGCGCAGGGCGAATATGTTCCCGCGTCGCAggg                      >  2:1595516/1‑92 (MQ=255)
                                                                    gagaAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTg                     <  1:4539655/93‑1 (MQ=255)
                                                                    gagaAGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTg                     >  1:3669318/1‑93 (MQ=255)
                                                                       aaGTACAACAGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGaaa                  <  1:4346192/93‑1 (MQ=255)
                                                                        aGTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAAt                 >  2:4478919/1‑93 (MQ=255)
                                                                         gTACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATc                >  2:1432537/1‑93 (MQ=255)
                                                                          tACAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCa               <  1:282678/93‑1 (MQ=255)
                                                                           aCAACGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCAt              >  2:2863890/1‑93 (MQ=255)
                                                                             aaCGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGa            >  2:1473742/1‑93 (MQ=255)
                                                                             aaCGGAACAGTGGGGCAGGTCACGGCGTTTAACGGCCCCGGGGAGATGGGGGCGCAGGGCGAATATGTTCCCGCGCCGCAGGGGAAATcggga            >  1:1116092/1‑89 (MQ=255)
                                                                              aCGGAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGAc           <  2:1698013/93‑1 (MQ=255)
                                                                                ggAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCa         <  1:1338227/93‑1 (MQ=255)
                                                                                 gAACAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCAc        <  2:722619/93‑1 (MQ=255)
                                                                                    cAGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCt     <  1:1293588/93‑1 (MQ=255)
                                                                                     aGTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTg    >  1:674155/1‑93 (MQ=255)
                                                                                      gTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGa   <  1:3007047/93‑1 (MQ=255)
                                                                                      gTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGa   <  1:4436381/93‑1 (MQ=255)
                                                                                      gTGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCCGa   >  2:1729563/1‑93 (MQ=255)
                                                                                       tGGGGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGAc  <  1:1670409/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGAC  >  NC_000913/1060856‑1061035

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