Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,194,577 2 bp→GC coding (97‑98/774 nt) yjaD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,194,5770CG100.0% 437.8 / NA 141R33A (CGG→GCG) yjaDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (49/92);  total (49/92)
*NC_0009134,194,5780GC99.3% 438.3 / NA 141R33A (CGG→GCG) yjaDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (48/92);  minor base A (1/0);  total (49/92)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

ATCGTATAATTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAACGGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGC  >  NC_000913/4194485‑4194664
                                                                                            ||                                                                                      
     ataATTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCa                                                                                    <  2:433693/93‑1 (MQ=255)
     ataATTGAAAAATTAGATCAAGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGAGAAGCGGCa                                                                                    >  2:3813475/1‑93 (MQ=255)
       aaTTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCaaa                                                                                  >  2:1855449/1‑93 (MQ=255)
       aaTTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCaaa                                                                                  <  1:3764469/93‑1 (MQ=255)
         ttGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAtt                                                                                >  2:92546/1‑93 (MQ=255)
         ttGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAtt                                                                                <  2:4057585/93‑1 (MQ=255)
         ttGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAtt                                                                                >  1:2429339/1‑93 (MQ=255)
          tGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAttt                                                                               <  1:2612951/93‑1 (MQ=255)
           gAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTc                                                                              <  2:3481822/93‑1 (MQ=255)
             aaaaTTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGa                                                                            <  2:3854566/93‑1 (MQ=255)
             aaaaTTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGa                                                                            <  2:3908079/93‑1 (MQ=255)
               aaTTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATc                                                                          >  1:98863/1‑93 (MQ=255)
                aTTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCt                                                                         >  2:1743810/1‑93 (MQ=255)
                 ttAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCtt                                                                        >  2:2679923/1‑93 (MQ=255)
                   aGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgt                                                                      <  2:1314800/93‑1 (MQ=255)
                   aGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgt                                                                      <  2:772037/93‑1 (MQ=255)
                    gATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgtg                                                                     <  1:1411341/93‑1 (MQ=255)
                    gATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgtg                                                                     <  1:114590/93‑1 (MQ=255)
                    gATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgtg                                                                     <  1:334606/93‑1 (MQ=255)
                     aTCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTgg                                                                    <  2:2900522/93‑1 (MQ=255)
                     aTCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTgg                                                                    <  2:1672911/93‑1 (MQ=255)
                     aTCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTCGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTgg                                                                    <  2:959190/93‑1 (MQ=255)
                      tCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATTGCGAAGCGGCAAATTTCGATCTTGTggg                                                                   >  2:53135/1‑93 (MQ=255)
                      tCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTggg                                                                   <  1:699774/93‑1 (MQ=255)
                        aCGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTc                                                                 >  2:2210093/1‑93 (MQ=255)
                         cGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCa                                                                <  1:1780704/93‑1 (MQ=255)
                         cGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCa                                                                <  2:42157/93‑1 (MQ=255)
                         cGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGAGGCAAATTTCGATCTTGTGGGTCa                                                                >  1:860750/1‑93 (MQ=255)
                          ggCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAg                                                               <  2:2866625/93‑1 (MQ=255)
                           gCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgc                                                              >  1:2078836/1‑93 (MQ=255)
                            cTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcg                                                             <  2:4613603/93‑1 (MQ=255)
                            cTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcg                                                             <  1:3929055/93‑1 (MQ=255)
                             tggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgc                                                            <  2:3021958/93‑1 (MQ=255)
                             tggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgc                                                            >  1:3693969/1‑93 (MQ=255)
                              ggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgcg                                                           >  1:738846/1‑93 (MQ=255)
                              ggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgcg                                                           >  1:4364474/1‑93 (MQ=255)
                                  ggTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACt                                                       <  2:3083436/93‑1 (MQ=255)
                                   gtcgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTa                                                      >  2:976607/1‑93 (MQ=255)
                                    tcgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTAc                                                     <  2:2710497/93‑1 (MQ=255)
                                     cgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACa                                                    >  1:640577/1‑93 (MQ=255)
                                      gtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACaa                                                   <  2:2609760/93‑1 (MQ=255)
                                       tcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACaaa                                                  >  2:3531928/1‑93 (MQ=255)
                                        cAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAAt                                                 <  2:1606343/93‑1 (MQ=255)
                                        cAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAAt                                                 <  1:2041165/93‑1 (MQ=255)
                                        cAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAAt                                                 <  1:909643/93‑1 (MQ=255)
                                         aGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATc                                                <  2:4038140/93‑1 (MQ=255)
                                          gCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCg                                               <  2:2474178/93‑1 (MQ=255)
                                            cATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGc                                             <  1:3285334/93‑1 (MQ=255)
                                            cATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGc                                             >  2:565224/1‑93 (MQ=255)
                                             aTGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCg                                            >  2:2638515/1‑93 (MQ=255)
                                             aTGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCg                                            >  2:123450/1‑93 (MQ=255)
                                              tGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGa                                           >  1:961941/1‑93 (MQ=255)
                                              tGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGa                                           >  1:1607706/1‑93 (MQ=255)
                                               gAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGaa                                          <  1:1648644/93‑1 (MQ=255)
                                               gAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGaa                                          <  2:2261581/93‑1 (MQ=255)
                                                aaCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAAt                                         <  1:2388617/93‑1 (MQ=255)
                                                aaCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAAt                                         <  2:3273056/93‑1 (MQ=255)
                                                aaCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAAt                                         <  2:3816989/93‑1 (MQ=255)
                                                 aCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATg                                        >  1:3017400/1‑93 (MQ=255)
                                                  cAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                                       <  1:10498/93‑1 (MQ=255)
                                                  cAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                                       <  1:2595415/93‑1 (MQ=255)
                                                  cAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                                       <  1:99379/93‑1 (MQ=255)
                                                   aaaaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGc                                      >  1:2167604/1‑93 (MQ=255)
                                                     aaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAg                                    <  2:1418993/93‑1 (MQ=255)
                                                     aaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAg                                    <  1:3644280/93‑1 (MQ=255)
                                                     aaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAg                                    <  1:1627385/93‑1 (MQ=255)
                                                      aaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgg                                   <  2:1499009/93‑1 (MQ=255)
                                                      aaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgg                                   <  1:3516925/93‑1 (MQ=255)
                                                       aTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAggg                                  <  1:280498/93‑1 (MQ=255)
                                                       aTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAggg                                  <  1:2223319/93‑1 (MQ=255)
                                                        ttATGGTTGCCGAAGGTAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgggg                                 >  2:759310/1‑93 (MQ=255)
                                                        ttATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgggg                                 <  1:1179028/93‑1 (MQ=255)
                                                        ttATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgggg                                 <  2:3667145/93‑1 (MQ=255)
                                                            ggTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAAc                             >  1:589120/1‑93 (MQ=255)
                                                             gTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAAcc                            <  1:4291253/93‑1 (MQ=255)
                                                             gTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAAcc                            <  2:3708807/93‑1 (MQ=255)
                                                              ttGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCt                           <  1:3025232/93‑1 (MQ=255)
                                                               tGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTg                          <  1:2757385/93‑1 (MQ=255)
                                                                gCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGt                         >  1:2000692/1‑93 (MQ=255)
                                                                gCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGt                         <  2:596183/93‑1 (MQ=255)
                                                                 ccGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGtt                        <  1:148007/93‑1 (MQ=255)
                                                                  cGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGttt                       <  1:877315/93‑1 (MQ=255)
                                                                  cGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGttt                       >  1:105660/1‑93 (MQ=255)
                                                                   gAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTg                      >  2:2479849/1‑93 (MQ=255)
                                                                   gAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTg                      >  1:2975440/1‑93 (MQ=255)
                                                                    aaGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTgg                     <  1:344213/93‑1 (MQ=255)
                                                                      gggAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGtt                   >  2:2141689/1‑93 (MQ=255)
                                                                      gggAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGAGCGCACTACAAATCGGCGAATGCCAGGAGGAACCTGTTTGGtt                   >  2:3367346/1‑93 (MQ=255)
                                                                       ggAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTa                  <  1:3063418/93‑1 (MQ=255)
                                                                       ggAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTa                  <  1:4142050/93‑1 (MQ=255)
                                                                             ttGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTacaa            <  2:738560/93‑1 (MQ=255)
                                                                                  atatGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCCCTACAAATCGGCGAATGGCAGGGGGGAACTGTTTGGTTAGTacaaccaca       >  2:4473731/1‑89 (MQ=255)
                                                                                  atatGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcagca       <  1:3351308/93‑1 (MQ=255)
                                                                                  atatGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcagca       <  1:449277/93‑1 (MQ=255)
                                                                                  atatGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcagca       <  2:2210353/93‑1 (MQ=255)
                                                                                   tatGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcagcag      <  1:101065/93‑1 (MQ=255)
                                                                                   tatGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcagcag      <  2:3357450/93‑1 (MQ=255)
                                                                                   tatGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcagcag      <  1:4174785/93‑1 (MQ=255)
                                                                                     tGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAAAGCAgcg    >  2:1211217/1‑93 (MQ=255)
                                                                                       gCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAgcggc  <  2:2106275/93‑1 (MQ=255)
                                                                                            ||                                                                                      
ATCGTATAATTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAACGGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGC  >  NC_000913/4194485‑4194664

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