Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,622,796 T→G 100% L1202R (CTG→CGG)  rhsB → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,622,7960TG98.9% 456.9 / ‑2.5 174L1202R (CTG→CGG) rhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (92/79);  total (93/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

Only 100 of 178 total aligned reads displayed.

ATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGAT  >  NC_000913/3622700‑3622891
                                                                                                |                                                                                               
atgatgAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTa                                                                                              <  1:653887/100‑1 (MQ=25)
atgatgAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTa                                                                                              <  1:510886/100‑1 (MQ=25)
atgatgAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTa                                                                                              <  2:657076/100‑1 (MQ=25)
atgatgAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTa                                                                                              >  2:308312/1‑100 (MQ=25)
 tgatgaATGGGGCAACCTGCTGAATGAAGAGAGCCCGGATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAt                                                                                             <  1:856586/100‑1 (MQ=255)
  natgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAt                                                                                             >  1:4081527/2‑99 (MQ=17)
  gatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGTTTATCCGGCTGTCGGGGCAGCAGTATGATGAGTAGTCAGGCCGGTAtt                                                                                            >  1:1864942/1‑100 (MQ=11)
  gatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAtt                                                                                            <  1:1399964/100‑1 (MQ=18)
  gatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAtt                                                                                            >  2:2560932/1‑100 (MQ=18)
     gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                         <  1:289032/100‑1 (MQ=18)
     gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                         >  1:2377365/1‑100 (MQ=18)
      aaTGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACaa                                                                                        >  1:2668752/1‑100 (MQ=17)
        tGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAAcc                                                                                      >  2:3304336/1‑100 (MQ=18)
         ggggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCg                                                                                     <  2:843002/100‑1 (MQ=18)
           ggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGcc                                                                                   <  2:1383318/100‑1 (MQ=18)
            gCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCa                                                                                  <  2:612299/100‑1 (MQ=18)
             cAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCAGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAc                                                                                 >  2:3378435/1‑100 (MQ=255)
              aaCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAcc                                                                                <  2:1109381/100‑1 (MQ=18)
               aCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCg                                                                               <  1:429806/100‑1 (MQ=17)
                ccTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAcc                                                                                <  2:1548189/98‑1 (MQ=17)
                  tgctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCta                                                                            <  1:3181933/100‑1 (MQ=18)
                   gctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtat                                                                           <  1:1920514/100‑1 (MQ=18)
                     tgaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtatta                                                                         >  1:3238817/1‑100 (MQ=17)
                         tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAc                                                                     >  1:588876/1‑100 (MQ=17)
                         tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAc                                                                     <  2:1954333/100‑1 (MQ=18)
                           aaGAGAACCCGCATCAGCTGCAGCCGCTTATCCAGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAccc                                                                   >  1:805947/1‑100 (MQ=255)
                           aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAccc                                                                   >  1:1392971/1‑100 (MQ=18)
                           aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAccc                                                                   >  1:2177128/1‑100 (MQ=17)
                           aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGGATTACAACCGCCACCTCTATTATGAccc                                                                   >  1:4245983/1‑100 (MQ=255)
                               gaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTg                                                               <  2:280957/100‑1 (MQ=17)
                                aaCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGc                                                              <  1:32456/100‑1 (MQ=17)
                                 aCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCa                                                             <  2:2377371/100‑1 (MQ=18)
                                  cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                            >  1:3842797/1‑100 (MQ=18)
                                  cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                            >  1:3513718/1‑100 (MQ=18)
                                    cGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAggg                                                          >  2:3910229/1‑100 (MQ=18)
                                    cGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAggg                                                          >  2:1444498/1‑100 (MQ=17)
                                    cGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAggg                                                          <  2:732210/100‑1 (MQ=18)
                                       aTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCg                                                       <  1:1853224/100‑1 (MQ=17)
                                          aGCTGCAGCAGCTTCTCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACACCCGCCCCCGCTTTTATGACCCGCTCCacg                                                           >  1:4226519/1‑91 (MQ=11)
                                          aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                    <  2:3541496/100‑1 (MQ=17)
                                          aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                    >  2:3878411/1‑100 (MQ=18)
                                          aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                    <  2:454851/100‑1 (MQ=18)
                                          aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGGATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                    <  2:1011409/100‑1 (MQ=12)
                                           gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatat                                                   >  2:1625131/1‑100 (MQ=18)
                                           gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatat                                                   >  1:2820433/1‑100 (MQ=18)
                                            cTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatata                                                  >  1:3230366/1‑100 (MQ=18)
                                             tGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatatat                                                 <  1:2541366/100‑1 (MQ=17)
                                              gcagcaGCTTATCCGGCTGCCTGCGCACCTGTATCATGAGTAGTCCGGCCAGTATTACAACCTCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                <  1:2476771/100‑1 (MQ=11)
                                              gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                <  2:1007241/100‑1 (MQ=17)
                                              gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                >  1:95545/1‑100 (MQ=18)
                                               cagcagCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCa                                               >  2:2684751/1‑100 (MQ=18)
                                               cagcagCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCa                                               >  1:3630154/1‑100 (MQ=18)
                                                agcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCAc                                              <  2:1394891/100‑1 (MQ=18)
                                                 gcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACt                                             >  1:2981834/1‑100 (MQ=17)
                                                    gcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAg                                          >  2:1307222/1‑100 (MQ=18)
                                                     cTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAgg                                         >  1:321182/1‑100 (MQ=18)
                                                     cTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAgg                                         <  2:210670/100‑1 (MQ=18)
                                                      ttATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGa                                        <  2:3860566/100‑1 (MQ=17)
                                                       tATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGAt                                       >  1:1686492/1‑100 (MQ=17)
                                                        aTCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATc                                      >  1:239184/1‑100 (MQ=18)
                                                        aTCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGGATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGAt                                       <  2:1848663/99‑1 (MQ=255)
                                                         tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGGCCCGCTGCAGGGGCGATATATCACTCAGGATcc                                     <  1:2781394/100‑1 (MQ=255)
                                                         tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATcc                                     <  2:3231611/100‑1 (MQ=18)
                                                         tCCGGCTGCCGGAGCAGCAGTATTATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGCCCCGCTGCAGGGGCGATATATCACTCAGGATcc                                     <  1:2499313/100‑1 (MQ=11)
                                                          ccGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATcc                                     <  1:3987345/99‑1 (MQ=17)
                                                          ccGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCg                                    >  2:3326228/1‑100 (MQ=17)
                                                            ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt                                  >  2:3432345/1‑100 (MQ=14)
                                                            ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt                                  <  1:1307218/100‑1 (MQ=14)
                                                                 ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                             >  1:2698290/1‑100 (MQ=14)
                                                                 ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                             >  2:4124824/1‑100 (MQ=14)
                                                                 ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                             <  1:2711422/100‑1 (MQ=14)
                                                                    gggCAGCAGTATGATGAGGAGTCCGGCCGGTATTACCACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGa                          <  2:194647/100‑1 (MQ=12)
                                                                       cagcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgg                       >  2:2595123/1‑100 (MQ=25)
                                                                       cagcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgg                       >  2:4187390/1‑100 (MQ=25)
                                                                       cagcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgg                       >  2:1379254/1‑100 (MQ=25)
                                                                         gcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggg                     >  2:3439534/1‑100 (MQ=25)
                                                                          cagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAggggg                    <  2:3644422/100‑1 (MQ=25)
                                                                           agTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                   <  2:1622689/100‑1 (MQ=255)
                                                                              atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                <  2:2799410/100‑1 (MQ=255)
                                                                               tgatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATgg               <  1:1567366/100‑1 (MQ=255)
                                                                                gatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa              >  1:631784/1‑100 (MQ=25)
                                                                                gatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa              >  2:1736455/1‑100 (MQ=255)
                                                                                 atgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGaa             >  1:3478187/1‑100 (MQ=255)
                                                                                  tgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAc            <  2:3603729/100‑1 (MQ=255)
                                                                                  tgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAc            >  2:4018005/1‑100 (MQ=38)
                                                                                  tgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAc            <  2:3171925/100‑1 (MQ=255)
                                                                                   gaggagTCCGGGCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACGCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc           >  1:2240453/1‑100 (MQ=255)
                                                                                   gaggagTACGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGAGAtcggaag           >  2:2240051/1‑94 (MQ=11)
                                                                                   gaggagTACGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACACAGGATCCGATTGGACTGAAGGGGGGATGGAAcc           >  1:2241239/1‑100 (MQ=255)
                                                                                     ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         >  1:1882920/1‑100 (MQ=255)
                                                                                     ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         >  2:718730/1‑100 (MQ=255)
                                                                                     ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         >  1:105633/1‑100 (MQ=255)
                                                                                       agTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTa       <  1:229247/100‑1 (MQ=255)
                                                                                         tCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATg     >  1:3866285/1‑100 (MQ=255)
                                                                                          ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATgg    >  2:2597565/1‑100 (MQ=255)
                                                                                          ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATgg    >  1:1466487/1‑100 (MQ=255)
                                                                                           cGGCCGGTATTACAACCGCCAGCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa   <  2:752630/100‑1 (MQ=255)
                                                                                           cGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGCTGGAACCTGTATGGa   >  2:299605/1‑100 (MQ=255)
                                                                                           cGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa   >  2:4167143/1‑100 (MQ=255)
                                                                                           cGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa   >  1:1080550/1‑100 (MQ=255)
                                                                                                |                                                                                               
ATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGAT  >  NC_000913/3622700‑3622891

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