Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,656,923 T→C P277P (CCT→CCC b1587 → putative oxidoreductase major subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,656,9230TC99.4% 562.6 / ‑5.1 177P277P (CCT→CCCb1587putative oxidoreductase major subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (96/80);  total (97/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 186 total aligned reads displayed.

AAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAA  >  NC_000913/1656834‑1657012
                                                                                         |                                                                                         
aaaGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCAtt                                                                                        <  2:4661314/93‑1 (MQ=255)
  aGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTcc                                                                                      >  2:3905763/1‑92 (MQ=255)
   gcgcgcGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAAGGGCTCACCATTCGc                                                                                     >  1:2105389/1‑93 (MQ=255)
   gcgcgcGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGGCGTGAAGACGAATGGCTCCCCATTCGc                                                                                     >  1:1271968/1‑93 (MQ=255)
   gcgcgcGAAAAATCGAACGCCAAAATGATTGTTACCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGc                                                                                     <  1:2513438/93‑1 (MQ=255)
    cgcgcgAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCTCATTCGcc                                                                                    <  2:6134463/93‑1 (MQ=255)
        cgAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCggg                                                                                >  1:795230/1‑93 (MQ=255)
          agaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCggggc                                                                              >  1:3337960/1‑91 (MQ=255)
          agaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGcc                                                                              >  2:4791952/1‑92 (MQ=255)
          agaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCa                                                                              <  1:5586032/93‑1 (MQ=255)
          agaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCa                                                                              >  2:2671644/1‑93 (MQ=255)
           gaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCAc                                                                             <  2:115928/93‑1 (MQ=255)
             aaTCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCg                                                                           >  2:3260550/1‑93 (MQ=255)
              aTCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGa                                                                          <  2:1736836/93‑1 (MQ=255)
              aTCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGa                                                                          >  1:6276349/1‑93 (MQ=255)
               tCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGAt                                                                         <  1:5471706/93‑1 (MQ=255)
               tCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGAt                                                                         <  1:1337915/93‑1 (MQ=255)
                 gAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGc                                                                       >  2:4410431/1‑93 (MQ=255)
                 gAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGc                                                                       <  1:5927161/93‑1 (MQ=255)
                   aCGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCcg                                                                     <  1:1543248/93‑1 (MQ=255)
                    cGCAAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCcgc                                                                    <  2:216161/93‑1 (MQ=255)
                       cAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGGCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCCCt                                                                 >  2:1440066/1‑93 (MQ=255)
                        aaaaTGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTg                                                                >  1:1105901/1‑93 (MQ=255)
                        aaaaTGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTg                                                                >  2:2543634/1‑93 (MQ=255)
                        aaaaTGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTg                                                                >  2:3851559/1‑93 (MQ=255)
                         aaaTGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTgg                                                               >  1:936505/1‑93 (MQ=255)
                         aaaTGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTgg                                                               >  1:1679048/1‑93 (MQ=255)
                           aTGATTGTTATCGATCCGCGATATACGGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCTATGCCGCGCTggcg                                                             >  2:458216/1‑91 (MQ=255)
                            tGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCCTTCGCCCGGGCACCGATGCCGCGCTGGTTg                                                            >  1:1170254/1‑93 (MQ=255)
                            tGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTg                                                            >  1:3613545/1‑93 (MQ=255)
                            tGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTg                                                            >  1:5133615/1‑93 (MQ=255)
                             gATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGc                                                           >  1:4919873/1‑93 (MQ=255)
                             gATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGCGAAGACGAATGGCTACCCATTCGCCCGGGCACCGATGCCGCGCCGGCTGc                                                           >  2:2415584/1‑93 (MQ=255)
                              aTTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCg                                                          >  1:4706288/1‑93 (MQ=255)
                              aTTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCg                                                          <  1:5111301/93‑1 (MQ=255)
                               ttGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCgg                                                         >  2:4747634/1‑93 (MQ=255)
                                tGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCggg                                                        <  2:663013/93‑1 (MQ=255)
                                tGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCggg                                                        >  1:2416776/1‑93 (MQ=255)
                                 gTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCgggg                                                       >  2:2791223/1‑92 (MQ=255)
                                  ttATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTa                                                      >  2:1597211/1‑93 (MQ=255)
                                  ttATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTa                                                      <  2:163551/93‑1 (MQ=255)
                                   tATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTAt                                                     >  1:3415073/1‑93 (MQ=255)
                                   tATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTAt                                                     >  1:4490513/1‑93 (MQ=255)
                                    aTCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGGGGGTAtt                                                    >  2:565937/1‑93 (MQ=255)
                                        aTCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCt                                                >  2:4916535/1‑93 (MQ=255)
                                        aTCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCt                                                <  1:2763280/93‑1 (MQ=255)
                                        aTCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCt                                                <  2:4651446/93‑1 (MQ=255)
                                        aTCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCt                                                >  2:425003/1‑93 (MQ=255)
                                         tCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTg                                               >  2:1304387/1‑93 (MQ=255)
                                         tCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTg                                               <  2:1908734/93‑1 (MQ=255)
                                            gcgATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGt                                            >  2:5444647/1‑93 (MQ=255)
                                             cgATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTa                                           <  1:3222634/93‑1 (MQ=255)
                                             cgATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTa                                           <  1:2498391/93‑1 (MQ=255)
                                              gATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTAt                                          >  2:4161307/1‑93 (MQ=255)
                                              gATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTAt                                          <  2:2289628/93‑1 (MQ=255)
                                                 ataCCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGGATTGa                                       >  1:1784723/1‑93 (MQ=255)
                                                  taCCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGAt                                      >  1:3693514/1‑93 (MQ=255)
                                                  taCCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGAt                                      >  1:2326360/1‑93 (MQ=255)
                                                   aCCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGAtt                                     <  1:765273/93‑1 (MQ=255)
                                                   aCCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGAtt                                     <  2:2959578/93‑1 (MQ=255)
                                                      gATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAAc                                  >  2:3165385/1‑93 (MQ=255)
                                                       aTACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACg                                 >  2:3880506/1‑93 (MQ=255)
                                                        tACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGa                                >  2:2613655/1‑93 (MQ=255)
                                                           ggCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGaaaa                             <  1:3885208/93‑1 (MQ=255)
                                                            gCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAt                            >  1:673206/1‑93 (MQ=255)
                                                             cTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtc                           <  1:703660/93‑1 (MQ=255)
                                                             cTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtc                           <  2:4371332/93‑1 (MQ=255)
                                                              tGCGGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtct                          <  1:893956/93‑1 (MQ=255)
                                                              tGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtct                          <  1:2356681/93‑1 (MQ=255)
                                                              tGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtct                          >  2:4088559/1‑93 (MQ=255)
                                                              tGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtct                          <  1:3915794/93‑1 (MQ=255)
                                                               gCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtctc                         <  2:4847940/93‑1 (MQ=255)
                                                                ccGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCg                        <  2:4436967/93‑1 (MQ=255)
                                                                 cGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGt                       >  2:2074946/1‑93 (MQ=255)
                                                                 cGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGt                       >  1:4081771/1‑93 (MQ=255)
                                                                 cGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGt                       <  2:5370037/93‑1 (MQ=255)
                                                                  ggTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGtt                      >  2:337417/1‑93 (MQ=255)
                                                                    tCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGa                    >  2:2015081/1‑93 (MQ=255)
                                                                       tGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCa                 >  2:3492933/1‑93 (MQ=255)
                                                                        gAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCaa                <  1:1859447/93‑1 (MQ=255)
                                                                        gAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCaa                <  1:2613962/93‑1 (MQ=255)
                                                                        gAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCaa                >  1:349954/1‑93 (MQ=255)
                                                                        gAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCaa                <  2:624457/93‑1 (MQ=255)
                                                                         aaGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAAc               >  1:1063473/1‑93 (MQ=255)
                                                                         aaGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTAGTGATTAGCGAAAATCTCGTTGATCAAc               <  2:5885647/93‑1 (MQ=255)
                                                                             cGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCttt           >  1:4938888/1‑93 (MQ=255)
                                                                             cGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCttt           <  1:1582039/93‑1 (MQ=255)
                                                                             cGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCttt           >  1:5236099/1‑93 (MQ=255)
                                                                              gAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCtttt          <  2:3990993/93‑1 (MQ=255)
                                                                              gAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCtttt          >  1:6151926/1‑93 (MQ=255)
                                                                               aaTGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTtc         >  1:4697277/1‑93 (MQ=255)
                                                                               aaTGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTtc         >  2:5107070/1‑93 (MQ=255)
                                                                                aTGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTtct        <  2:2845766/93‑1 (MQ=255)
                                                                                 tGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTtctc       >  2:3414449/1‑93 (MQ=255)
                                                                                  ggCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCg      <  1:3515307/93‑1 (MQ=255)
                                                                                  ggCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCg      <  1:502023/93‑1 (MQ=255)
                                                                                    cTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGAt    <  2:1571391/93‑1 (MQ=255)
                                                                                    cTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGAt    <  1:4746760/93‑1 (MQ=255)
                                                                                     tCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATa   <  2:5607371/93‑1 (MQ=255)
                                                                                      ccccATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATaa  <  2:5477580/93‑1 (MQ=255)
                                                                                         |                                                                                         
AAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAA  >  NC_000913/1656834‑1657012

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