Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 841,817 A→G G80G (GGT→GGC mcbA ← DUF1471 domain‑containing protein McbA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913841,8170AG97.2% 1125.0 / 12.6 286G80G (GGT→GGCmcbADUF1471 domain‑containing protein McbA
Reads supporting (aligned to +/- strand):  ref base A (6/2);  new base G (145/131);  total (151/135)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 286 total aligned reads displayed.

AGTTCATATTTGATGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTACCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTA  >  NC_000913/841719‑841914
                                                                                                  |                                                                                                 
aGTTCATATTTGATGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTAcc                                                                                                 >  2:1308147/1‑101 (MQ=255)
  ttCATATTTGATGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGa                                                                                               >  1:1308149/1‑101 (MQ=255)
   tCATATTTGATGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGaa                                                                                              <  2:1307964/101‑1 (MQ=255)
     atatTTGATGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTACCGAACa                                                                                            >  1:1308153/1‑101 (MQ=255)
      tatTTGATGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACAt                                                                                           >  2:1308155/1‑101 (MQ=255)
         ttGATGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTg                                                                                        <  2:1307463/101‑1 (MQ=255)
         ttGATGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTg                                                                                        <  2:1307749/101‑1 (MQ=255)
         ttGATGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTg                                                                                        >  2:1308159/1‑101 (MQ=255)
          tGATGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTgg                                                                                       >  2:1308163/1‑101 (MQ=255)
            aTGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTc                                                                                     >  2:1308165/1‑101 (MQ=255)
             tGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCa                                                                                    <  1:1307729/101‑1 (MQ=255)
              gTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTACCGAACATCTGGTCAt                                                                                   >  1:1308166/1‑101 (MQ=255)
                tttttGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCAttt                                                                                 <  1:1307630/101‑1 (MQ=255)
                tttttGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCAttt                                                                                 >  1:1308171/1‑101 (MQ=255)
                tttttGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTACCGAACATCTGGTCAttt                                                                                 <  1:1307410/101‑1 (MQ=255)
                 ttttGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTc                                                                                >  2:1308172/1‑101 (MQ=255)
                 ttttGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTc                                                                                <  1:1307716/101‑1 (MQ=255)
                   ttGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCg                                                                              >  1:1308174/1‑101 (MQ=255)
                    tGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGc                                                                             >  2:1308179/1‑101 (MQ=255)
                      aTTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCt                                                                           <  2:1307514/101‑1 (MQ=255)
                      aTTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCt                                                                           <  2:1307808/101‑1 (MQ=255)
                       ttACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTg                                                                          >  2:1308185/1‑101 (MQ=255)
                        tACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGc                                                                         >  1:1308187/1‑101 (MQ=255)
                         aCAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTgg                                                                                       >  1:1308191/1‑86 (MQ=255)
                          cAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCgg                                                                       <  1:1307720/101‑1 (MQ=255)
                          cAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCgg                                                                       <  1:1307614/101‑1 (MQ=255)
                           aaaaCAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGa                                                                      >  1:1308194/1‑101 (MQ=255)
                           aaaaCAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGa                                                                      <  1:1307580/101‑1 (MQ=255)
                           aaaaCAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGa                                                                      <  1:1307803/101‑1 (MQ=255)
                            aaaCAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGaa                                                                     >  2:1308198/1‑101 (MQ=255)
                            aaaCAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGaa                                                                     >  1:1308200/1‑101 (MQ=255)
                             aaCAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAAt                                                                    >  2:1308202/1‑101 (MQ=255)
                              aCAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAAtt                                                                   <  2:1308051/101‑1 (MQ=255)
                               cAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTa                                                                  <  2:1307734/101‑1 (MQ=255)
                                aGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTaa                                                                 >  2:1308208/1‑101 (MQ=255)
                                aGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTaa                                                                 >  1:1308206/1‑101 (MQ=255)
                                  aCGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATg                                                               <  1:1307687/101‑1 (MQ=255)
                                  aCGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATg                                                               <  1:1307762/101‑1 (MQ=255)
                                  aCGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATg                                                               <  1:1307694/101‑1 (MQ=255)
                                  aCGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTACCGAACATCTGGTCATTTCCGCCTGCGGAATTAATg                                                               >  1:1308213/1‑101 (MQ=255)
                                  aCGAAAACTGCATCCAGAGGAGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATg                                                               >  2:1308212/1‑101 (MQ=255)
                                     aaaaCTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACa                                                            >  1:1308218/1‑101 (MQ=255)
                                     aaaaCTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACa                                                            <  1:1307535/101‑1 (MQ=255)
                                      aaaCTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACAt                                                           >  2:1308219/1‑101 (MQ=255)
                                       aaCTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATa                                                          <  1:1307932/101‑1 (MQ=255)
                                           gCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCt                                                      >  1:1308224/1‑101 (MQ=255)
                                           gCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCt                                                      >  2:1308225/1‑101 (MQ=255)
                                            cATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCtt                                                     >  2:1308230/1‑101 (MQ=255)
                                            cATCCAGAGGCGAATGCAGTTACGCTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCtt                                                     <  2:1307876/101‑1 (MQ=255)
                                             aTCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCttt                                                    <  2:1307656/101‑1 (MQ=255)
                                             aTCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCttt                                                    <  1:1307437/101‑1 (MQ=255)
                                             aTCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCttt                                                    >  1:1308231/1‑101 (MQ=255)
                                                cAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTcgc                                                 <  1:1307756/101‑1 (MQ=255)
                                                 agagGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTcgcg                                                >  1:1308234/1‑101 (MQ=255)
                                                    ggCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCt                                             <  1:1307673/101‑1 (MQ=255)
                                                    ggCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCt                                             <  2:1307671/101‑1 (MQ=255)
                                                    ggCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCt                                             <  2:1307787/101‑1 (MQ=255)
                                                    ggCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCt                                             >  1:1308239/1‑101 (MQ=255)
                                                    ggCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCt                                             <  2:1307785/101‑1 (MQ=255)
                                                     gCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCtt                                            <  2:1307633/101‑1 (MQ=255)
                                                      cGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTg                                           >  2:1308240/1‑101 (MQ=255)
                                                       gAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGt                                          >  1:1308244/1‑101 (MQ=255)
                                                       gAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTACCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGt                                          <  1:1307828/101‑1 (MQ=255)
                                                        aaTGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGtt                                         <  1:1307774/101‑1 (MQ=255)
                                                          tGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTcg                                       <  2:1307759/101‑1 (MQ=255)
                                                          tGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTcg                                       >  1:1308247/1‑101 (MQ=255)
                                                             aGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTcgcgc                                    <  1:1307845/101‑1 (MQ=255)
                                                             aGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTcgcgc                                    >  2:1308248/1‑101 (MQ=255)
                                                               ttACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCgcgcgc                                  >  2:1308253/1‑101 (MQ=255)
                                                                 aCGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCtt                                >  2:1308255/1‑101 (MQ=255)
                                                                 aCGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCtt                                >  1:1308257/1‑101 (MQ=255)
                                                                 aCGTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCtt                                <  1:1308178/101‑1 (MQ=255)
                                                                   gTTGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCtttt                              >  2:1308261/1‑101 (MQ=255)
                                                                    ttGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCttttt                             <  2:1307900/101‑1 (MQ=255)
                                                                    ttGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCttttt                             <  2:1307557/101‑1 (MQ=255)
                                                                    ttGGTGGGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCttttt                             <  1:1307866/101‑1 (MQ=255)
                                                                     tggtggGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTc                            >  1:1308263/1‑101 (MQ=255)
                                                                      ggtggGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCt                           <  2:1307827/101‑1 (MQ=255)
                                                                      ggtggGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCt                           <  1:1307879/101‑1 (MQ=255)
                                                                       gtggGGTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTg                          >  2:1308265/1‑101 (MQ=255)
                                                                         ggggTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGcc                        >  1:1308269/1‑101 (MQ=255)
                                                                          gggTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCa                       <  2:1307769/101‑1 (MQ=255)
                                                                            gTTATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGt                     >  2:1308273/1‑101 (MQ=255)
                                                                             ttATTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGtt                    >  1:1308274/1‑101 (MQ=255)
                                                                              tATTTGTAGATGGGTGCAGTACCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGttt                   >  1:1308279/1‑101 (MQ=255)
                                                                               aTTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTa                  >  2:1308280/1‑101 (MQ=255)
                                                                               aTTTGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTa                  >  1:1308283/1‑101 (MQ=255)
                                                                                tttGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTAt                 >  1:1308286/1‑101 (MQ=255)
                                                                                tttGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTAt                 <  2:1307812/101‑1 (MQ=255)
                                                                                tttGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTAt                 <  1:1307746/101‑1 (MQ=255)
                                                                                 ttGTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTAtc                >  1:1308288/1‑101 (MQ=255)
                                                                                   gTAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTAtcct              >  2:1308293/1‑101 (MQ=255)
                                                                                    tAGATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTAtcctc             <  2:1307846/101‑1 (MQ=255)
                                                                                      gATGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCa           <  2:1307794/101‑1 (MQ=255)
                                                                                        tGGTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGa         <  1:1307859/101‑1 (MQ=255)
                                                                                          gTTGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATc       <  1:1307740/101‑1 (MQ=255)
                                                                                           ttGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCg      <  2:1307707/101‑1 (MQ=255)
                                                                                            tGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGc     >  1:1308296/1‑101 (MQ=255)
                                                                                            tGCAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGc     <  1:1307886/101‑1 (MQ=255)
                                                                                              cAGTGCCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCtt   >  2:1308297/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
AGTTCATATTTGATGTTTTTTGATTACAAAACAGACGAAAACTGCATCCAGAGGCGAATGCAGTTACGTTGGTGGGGTTATTTGTAGATGGTTGCAGTACCGAACATCTGGTCATTTCCGCCTGCGGAATTAATGACATAACCTTTCGCGCCTTGTTCGCGCGCTTTTTCTGCCAGTTTATCCTCCAGATCGCTTA  >  NC_000913/841719‑841914

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