Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,057,390 C→T G28G (GGC→GGT torC → trimethylamine N‑oxide reductase, cytochrome c‑type subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,057,3900CT96.6% 590.3 / 3.4 178G28G (GGC→GGTtorCtrimethylamine N‑oxide reductase, cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base T (91/81);  total (93/85)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 182 total aligned reads displayed.

GGGGTTATTATGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTA  >  NC_000913/1057298‑1057479
                                                                                            |                                                                                         
ggggTTATTATGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGCTTGTGATTGGc                                                                                           >  1:2825865/1‑93 (MQ=255)
   gTTATTATGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTAtt                                                                                        <  2:2959175/93‑1 (MQ=255)
   gTTATTATGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTAtt                                                                                        >  2:2402786/1‑93 (MQ=255)
    ttattaTGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTg                                                                                       >  2:739610/1‑93 (MQ=255)
    ttattaTGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTg                                                                                       >  1:4104179/1‑93 (MQ=255)
     tattatGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTgc                                                                                      <  2:2937449/93‑1 (MQ=255)
     tattatGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTgc                                                                                      <  1:1711175/93‑1 (MQ=255)
      attatGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTgcg                                                                                     >  1:4002919/1‑93 (MQ=255)
      attatGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTgcg                                                                                     >  2:42946/1‑93 (MQ=255)
       ttatGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTgcgc                                                                                    >  2:4137538/1‑93 (MQ=255)
       ttatGCGGAAACTCTGGAACGCGCTACGCCGACCCAGCGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTgcgc                                                                                    >  1:3624415/1‑93 (MQ=255)
       ttatGCGGAAACCCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTgcgc                                                                                    <  1:508897/93‑1 (MQ=255)
           gCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGAt                                                                                >  1:382689/1‑93 (MQ=255)
             ggAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTg                                                                              <  2:3916542/93‑1 (MQ=255)
              gAAACTCTGGAACGCGCTACGCCGCCCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGt                                                                             <  2:1275755/93‑1 (MQ=255)
              gAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGt                                                                             <  2:3627986/93‑1 (MQ=255)
                 aCTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTAtt                                                                          <  1:195698/93‑1 (MQ=255)
                 aCTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTAtt                                                                          <  2:1224840/93‑1 (MQ=255)
                   tctGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGc                                                                        >  2:4256383/1‑93 (MQ=255)
                    ctGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGcc                                                                       <  2:2945761/93‑1 (MQ=255)
                     tGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGccc                                                                      >  2:1621589/1‑92 (MQ=255)
                     tGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCca                                                                      <  1:1281482/93‑1 (MQ=255)
                      ggAATGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCcac                                                                     >  2:2595925/1‑93 (MQ=255)
                       gAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCcaca                                                                    <  1:2592412/93‑1 (MQ=255)
                       gAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCcaca                                                                    >  1:4350458/1‑93 (MQ=255)
                       gAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCcaca                                                                    <  1:1811290/93‑1 (MQ=255)
                        aaCGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCacac                                                                   >  2:3934176/1‑93 (MQ=255)
                         aCGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACg                                                                  >  2:3265905/1‑93 (MQ=255)
                           gcgcTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGtt                                                                <  2:3368579/93‑1 (MQ=255)
                           gcgcTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGtt                                                                >  2:2140088/1‑93 (MQ=255)
                           gcgcTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGtt                                                                >  1:3778542/1‑93 (MQ=255)
                            cgcTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTg                                                               <  1:1546317/93‑1 (MQ=255)
                              cTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTGCGCTGATTGTATTGCCACACGTTggg                                                             <  2:3724887/93‑1 (MQ=255)
                                aCGCCGCCCCAGTGCTCGTTGGTCGTGACTGGCGCTGGTCGTAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGAt                                                           <  2:2064199/93‑1 (MQ=255)
                                aCGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGAt                                                           >  2:4472013/1‑93 (MQ=255)
                                aCGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGAt                                                           <  1:3763247/93‑1 (MQ=255)
                                 cGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATc                                                          <  1:892551/93‑1 (MQ=255)
                                 cGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATc                                                          >  2:4185645/1‑93 (MQ=255)
                                   ccGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCaa                                                        >  1:3303528/1‑93 (MQ=255)
                                    cGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCaaa                                                       >  1:1633565/1‑93 (MQ=255)
                                    cGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCaaa                                                       <  1:2078670/93‑1 (MQ=255)
                                     gACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAg                                                      >  1:2547448/1‑93 (MQ=255)
                                      aCCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGt                                                     >  1:103391/1‑93 (MQ=255)
                                      aCCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGt                                                     >  1:3954144/1‑93 (MQ=255)
                                       cccAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTc                                                    >  1:4521960/1‑93 (MQ=255)
                                       cccAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTc                                                    >  1:2042739/1‑93 (MQ=255)
                                         cAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCAc                                                  >  2:4415812/1‑93 (MQ=255)
                                           gTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCa                                                >  1:2849569/1‑93 (MQ=255)
                                           gTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCa                                                >  1:2378023/1‑93 (MQ=255)
                                           gTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCa                                                >  1:1936848/1‑93 (MQ=255)
                                             gCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGc                                              <  2:1380325/93‑1 (MQ=255)
                                              cTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGca                                             <  2:572447/93‑1 (MQ=255)
                                               tCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGcac                                            >  2:2961609/1‑93 (MQ=255)
                                               tCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGcac                                            >  2:92602/1‑93 (MQ=255)
                                               tCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGcac                                            <  2:2337099/93‑1 (MQ=255)
                                               tCGTTGGTCGGTACTGGCGCTGGTCGCAACTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCTCCAGcac                                            >  2:1970272/1‑93 (MQ=255)
                                                 gTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACaa                                          >  2:3572267/1‑93 (MQ=255)
                                                  ttGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAAc                                         >  2:1450859/1‑93 (MQ=255)
                                                  ttGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAAc                                         <  1:584518/93‑1 (MQ=255)
                                                  ttGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAAc                                         >  1:2584317/1‑93 (MQ=255)
                                                      tCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGaa                                     <  2:2216084/93‑1 (MQ=255)
                                                       cGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAAt                                    <  1:3424107/93‑1 (MQ=255)
                                                         gTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAGCCGAAttt                                  <  2:3310558/93‑1 (MQ=255)
                                                          tACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAAtttt                                 <  2:2904884/93‑1 (MQ=255)
                                                          tACTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAAtttt                                 <  1:104185/93‑1 (MQ=255)
                                                            cTGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTtgt                               <  1:2490971/93‑1 (MQ=255)
                                                             tGGCGCTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTtgtg                              >  1:3525496/1‑93 (MQ=255)
                                                               gcgcTGGTGGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTc                            <  2:483291/93‑1 (MQ=255)
                                                               gcgcTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTc                            <  2:2649861/93‑1 (MQ=255)
                                                                cgcTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCa                           <  2:1770917/93‑1 (MQ=255)
                                                                  cTGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGt                         <  1:2727756/93‑1 (MQ=255)
                                                                   tGGTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGtt                        <  1:4704008/93‑1 (MQ=255)
                                                                    ggTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTg                       >  1:4598190/1‑93 (MQ=255)
                                                                     gTCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGc                      >  1:2986203/1‑93 (MQ=255)
                                                                      tCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGcc                     <  1:2136344/93‑1 (MQ=255)
                                                                      tCGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGcc                     <  2:325934/93‑1 (MQ=255)
                                                                       cGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCca                    <  2:2622345/93‑1 (MQ=255)
                                                                       cGCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCca                    <  1:3261254/93‑1 (MQ=255)
                                                                        gCAATTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCcac                   >  2:2499847/1‑93 (MQ=255)
                                                                          aaTTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAg                 >  1:1064064/1‑93 (MQ=255)
                                                                           aTTGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGt                >  1:714987/1‑93 (MQ=255)
                                                                            ttGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTa               <  2:4542226/93‑1 (MQ=255)
                                                                            ttGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTa               >  1:3036229/1‑93 (MQ=255)
                                                                            ttGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTa               <  2:948631/93‑1 (MQ=255)
                                                                            ttGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTa               >  1:3585264/1‑93 (MQ=255)
                                                                            ttGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTa               >  2:1452147/1‑93 (MQ=255)
                                                                             tGGGATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTAt              <  2:3538175/93‑1 (MQ=255)
                                                                              gggATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATg             >  2:3206086/1‑93 (MQ=255)
                                                                              gggATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATg             >  2:1737891/1‑93 (MQ=255)
                                                                               ggATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGc            >  1:3349777/1‑93 (MQ=255)
                                                                                gATTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCa           >  1:49789/1‑93 (MQ=255)
                                                                                 aTTGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCaa          >  1:3379052/1‑93 (MQ=255)
                                                                                  ttGTGATTGGTATTTCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAAc         <  1:1898994/93‑1 (MQ=255)
                                                                                  ttGTGATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAAc         >  2:2274217/1‑93 (MQ=255)
                                                                                   tgtgATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAAcc        <  2:4655004/93‑1 (MQ=255)
                                                                                     tgATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCgg      >  1:1103890/1‑93 (MQ=255)
                                                                                      gATTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGgt     >  1:4278664/1‑93 (MQ=255)
                                                                                       aTTGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGgtg    >  1:1995537/1‑93 (MQ=255)
                                                                                        ttGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGgtgt   >  1:3647951/1‑93 (MQ=255)
                                                                                         tGGTATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTa  <  2:1033722/93‑1 (MQ=255)
                                                                                            |                                                                                         
GGGGTTATTATGCGGAAACTCTGGAACGCGCTACGCCGACCCAGTGCTCGTTGGTCGGTACTGGCGCTGGTCGCAATTGGGATTGTGATTGGCATTGCGCTGATTGTATTGCCACACGTTGGGATCAAAGTCACCAGCACAACCGAATTTTGTGTCAGTTGCCACAGTATGCAACCGGTGTA  >  NC_000913/1057298‑1057479

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