Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,766,247 A→C K1355N (AAA→AAC rhsA → Rhs protein with putative toxin 55 domain; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,766,2470AC98.8% 591.3 / ‑4.8 171K1355N (AAA→AACrhsARhs protein with putative toxin 55 domain; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (99/68);  total (99/72)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

GACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAA  >  NC_000913/3766152‑3766343
                                                                                               |                                                                                                 
gACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATg                                                                                               <  1:2803133/100‑1 (MQ=255)
  cATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGCTATTAGTGCTAACGATGCa                                                                                             >  2:2204248/1‑100 (MQ=255)
     aataaATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCaaaa                                                                                          >  2:3079773/1‑100 (MQ=255)
     aataaATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCaaaa                                                                                          <  1:4065692/100‑1 (MQ=255)
        aaaTGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGc                                                                                       <  2:2159683/100‑1 (MQ=255)
          aTGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGCGCTAACGATGCAAAAAGcac                                                                                     <  1:3179166/100‑1 (MQ=255)
            gACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGcacac                                                                                   <  1:1061414/100‑1 (MQ=255)
            gACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGcacac                                                                                   >  2:3643706/1‑100 (MQ=255)
            gACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGcacac                                                                                   <  1:2450758/100‑1 (MQ=255)
            gACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGcacac                                                                                   <  1:415162/100‑1 (MQ=255)
            gACAGATTCTGAAGTAATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGcacac                                                                                   <  2:1078515/100‑1 (MQ=255)
             aCAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGcacaca                                                                                  >  1:2203956/1‑100 (MQ=255)
             aCAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGcacaca                                                                                  >  1:3162731/1‑100 (MQ=255)
              cAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAg                                                                                 >  1:597567/1‑100 (MQ=255)
               aGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGa                                                                                >  1:3444038/1‑100 (MQ=255)
                gATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGaa                                                                               >  2:1889782/1‑100 (MQ=255)
                      ggAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTg                                                                         >  1:1788180/1‑100 (MQ=255)
                      gaaGTGATGGGAATAAGCCCGACAGATGTGATGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAAATAAAAAAAGCACACAGAAAGCCTg                                                                         <  2:3140304/98‑1 (MQ=255)
                         gTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGaa                                                                      >  1:2625152/1‑100 (MQ=255)
                         gTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGaa                                                                      >  1:998383/1‑100 (MQ=255)
                          tGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAAt                                                                     >  1:4060074/1‑100 (MQ=255)
                          tGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAAt                                                                     <  1:756068/100‑1 (MQ=255)
                           gATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAAtt                                                                    >  1:1377972/1‑100 (MQ=255)
                            aTGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTg                                                                   <  1:1000975/100‑1 (MQ=255)
                             tGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgt                                                                  >  1:1057672/1‑100 (MQ=255)
                              gggAAAAAACCCGACAGATGTGGCGTACTCCAGGCGCGAAGTGATTGCGCAGATATTAGTGCTAACGCTGCAAAAAGCACACAGAAAGCCTGGAATTggc                                                                 >  2:1741344/1‑98 (MQ=255)
                              gggAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgtc                                                                 >  2:1012389/1‑100 (MQ=255)
                              gggAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgtc                                                                 >  2:1883464/1‑100 (MQ=255)
                              gggAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgtc                                                                 >  2:2307827/1‑100 (MQ=255)
                              gggAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgtc                                                                 >  1:1149066/1‑100 (MQ=255)
                               ggAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgtcg                                                                <  2:1070218/100‑1 (MQ=255)
                               ggAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgtcg                                                                <  1:1758238/100‑1 (MQ=255)
                               ggAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgtcg                                                                >  1:3892893/1‑100 (MQ=255)
                               ggAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgtcg                                                                >  1:3524732/1‑100 (MQ=255)
                               agAAAAAACCCAACAGATGTGACGTACCGCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgtcg                                                                <  2:1580069/99‑1 (MQ=255)
                                 aaaaaaCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgtcgtc                                                              <  1:2801182/100‑1 (MQ=255)
                                 aaaaaaCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTgtcgtc                                                              <  1:3176315/100‑1 (MQ=255)
                                    aaaCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAGGCCTGGAATTGTCGTCACt                                                           >  1:3902960/1‑100 (MQ=255)
                                      aCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTcc                                                         <  2:3599783/100‑1 (MQ=255)
                                      aCCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTcc                                                         <  2:79905/100‑1 (MQ=255)
                                       cccGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCa                                                        <  2:522948/100‑1 (MQ=255)
                                       cccGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCa                                                        >  2:1614455/1‑100 (MQ=255)
                                        ccGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAg                                                       <  1:3692014/100‑1 (MQ=255)
                                        ccGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAg                                                       <  1:1036911/100‑1 (MQ=255)
                                         cGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGa                                                      >  1:689753/1‑100 (MQ=255)
                                         cGACAGATGTGACGTACTTCAGGAACGAATTGATTGTGGTGATATTATTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGa                                                      >  2:3802479/1‑100 (MQ=255)
                                          gACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGAc                                                     <  2:3860816/100‑1 (MQ=255)
                                              gATGTGACGTACTTCAAGAACTAATTGGTTGGGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTc                                                 >  2:696006/1‑100 (MQ=255)
                                                tgtgACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCaa                                               >  2:1218633/1‑100 (MQ=255)
                                                tgtgACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCaa                                               >  1:2556265/1‑100 (MQ=255)
                                                 gtgACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCaaa                                              >  1:1611696/1‑100 (MQ=255)
                                                     cGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGAt                                          >  1:2374103/1‑100 (MQ=255)
                                                      gTACTTTAGAAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGTAAGCCTGGAATTGTCGCCACTCCAGACAGTCAAACGATa                                        <  2:3412938/100‑1 (MQ=255)
                                                       tACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATaa                                       >  1:960376/1‑100 (MQ=255)
                                                       tACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATaa                                       >  2:3950869/1‑100 (MQ=255)
                                                         cTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATaaaa                                     >  2:2584924/1‑100 (MQ=255)
                                                          ttCAGGAACTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATaaaaa                                    <  2:2517068/100‑1 (MQ=255)
                                                             attaaCTAATTGATTGTGGTGATATTAGTTCTAACGATGTAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCTAGACAGTCAAACTATAAAAAAAg                                 <  2:2849371/97‑1 (MQ=255)
                                                                aaCTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATa                              <  1:3706302/100‑1 (MQ=255)
                                                                aaCTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATa                              <  2:2034326/100‑1 (MQ=255)
                                                                 aCTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAg                             <  2:1100693/100‑1 (MQ=255)
                                                                  cTAATTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGc                            <  1:2537766/100‑1 (MQ=255)
                                                                    aaTTGATTGTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGccc                          >  1:2911459/1‑100 (MQ=255)
                                                                     attgattgTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCt                         >  1:2709729/1‑100 (MQ=255)
                                                                       tgattgTGGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTtg                       <  1:1843861/100‑1 (MQ=255)
                                                                           tgtgGTGATATTAGTGCTAATGATGCAAAAAGCATACAGAAAGCTTGGAATTGTCGTTACTCCAGACAGTAAAACGATAAAAAAAGATAGCCCTTGTGGa                   <  2:215845/100‑1 (MQ=255)
                                                                           tgtgGTGATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGa                   <  2:539390/100‑1 (MQ=255)
                                                                            gtggtgATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAg                  >  2:3631291/1‑100 (MQ=255)
                                                                            gtggtgATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAg                  <  1:1302618/100‑1 (MQ=255)
                                                                            gtggtgATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAg                  >  2:1641164/1‑100 (MQ=255)
                                                                             tggtgATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAgg                 >  2:1174037/1‑100 (MQ=255)
                                                                              ggtgATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGt                <  2:3358841/100‑1 (MQ=255)
                                                                              ggtgATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGt                >  1:1535073/1‑100 (MQ=255)
                                                                               gtgATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGtt               >  2:585696/1‑100 (MQ=255)
                                                                                tgATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTc              >  1:1789478/1‑100 (MQ=255)
                                                                                 tatatTAGTGCAAACGATGCAAAAATCACACAGAAAGACTTGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTcc             <  2:3275513/99‑1 (MQ=37)
                                                                                 gATATTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTcc             <  1:2948633/100‑1 (MQ=255)
                                                                                   tatTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTg           <  1:1483531/100‑1 (MQ=255)
                                                                                   tatTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTg           <  1:1943676/100‑1 (MQ=255)
                                                                                   tatTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTg           >  1:1434881/1‑100 (MQ=255)
                                                                                   tatTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTg           >  1:1468308/1‑100 (MQ=255)
                                                                                    atTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCGCTCCAGACAGTAAAAAAAAAAAAAAAAAATAGCCCTTGTGGAGGTTCCTg           <  2:1230150/100‑1 (MQ=255)
                                                                                    atTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGc          <  2:1443020/100‑1 (MQ=255)
                                                                                    atTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGc          <  2:1442744/100‑1 (MQ=255)
                                                                                    atTAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGc          >  1:1604846/1‑100 (MQ=255)
                                                                                     ttAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCa         >  1:3830796/1‑100 (MQ=255)
                                                                                     ttAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCa         >  1:393248/1‑100 (MQ=255)
                                                                                      tAGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCaa        <  2:1323049/100‑1 (MQ=255)
                                                                                       aGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAAt       <  1:552185/100‑1 (MQ=255)
                                                                                       aGTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAAt       >  1:4030097/1‑100 (MQ=255)
                                                                                        gTGCTAACGATGCGAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATg      <  1:1557175/100‑1 (MQ=255)
                                                                                        gTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATg      >  2:497647/1‑100 (MQ=255)
                                                                                        gTGCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATg      >  2:498174/1‑100 (MQ=255)
                                                                                          gCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTc    >  1:1692883/1‑100 (MQ=255)
                                                                                          gCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTc    >  1:1066431/1‑100 (MQ=255)
                                                                                          gCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTc    <  1:2047415/100‑1 (MQ=255)
                                                                                          gCTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTc    <  1:20896/100‑1 (MQ=255)
                                                                                           cTAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCa   <  1:2813275/100‑1 (MQ=255)
                                                                                            tAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCaa  <  1:1942284/100‑1 (MQ=255)
                                                                                            tAACGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCaa  >  2:1806568/1‑100 (MQ=255)
                                                                                               |                                                                                                 
GACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAATTGATTGTGGTGATATTAGTGCTAAAGATGCAAAAAGCACACAGAAAGCCTGGAATTGTCGTCACTCCAGACAGTCAAACGATAAAAAAAGATAGCCCTTGTGGAGGTTCCTGCAATGTCAA  >  NC_000913/3766152‑3766343

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