Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,499 A→C 16.2% intergenic (+95/‑87) ydcD → / → yncI putative immunity protein for RhsE/pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,530,4990AC16.2% 411.9 / 76.6 167intergenic (+95/‑87)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base A (74/66);  new base C (12/15);  total (86/81)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

Only 100 of 175 total aligned reads displayed.

CATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAA  >  NC_000913/1530400‑1530598
                                                                                                   |                                                                                                   
cATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTaa                                                                                                     >  1:546262/1‑100 (MQ=255)
 aTAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTaaa                                                                                                    >  1:4068206/1‑100 (MQ=255)
 aTAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTaaa                                                                                                    >  1:835941/1‑100 (MQ=255)
 aTAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTaaa                                                                                                    <  2:1519098/100‑1 (MQ=255)
  tAATGAACGACACCAGAAAGAGAAAAATATGAAAAATATACAGAAAGACATGAGAGAAAAGTACAGGGCAGGTCACGATCATGAACATTGTGCTCTAAAt                                                                                                   >  2:1074041/1‑100 (MQ=255)
  tAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAt                                                                                                   >  1:3584458/1‑100 (MQ=255)
      gAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACg                                                                                               >  2:1795612/1‑100 (MQ=255)
      gAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACg                                                                                               >  2:2114738/1‑100 (MQ=255)
       aacaacACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGa                                                                                              >  1:2334835/1‑100 (MQ=255)
       aacaacACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGa                                                                                              <  2:232837/100‑1 (MQ=255)
        acaacaacaGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAAt                                                                                             <  2:220461/100‑1 (MQ=38)
        acaacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAt                                                                                             <  1:1825973/100‑1 (MQ=255)
          aacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAttt                                                                                           <  2:931159/100‑1 (MQ=255)
           acacCAGAAAGAGAAAAATATAAACAACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                          <  2:399078/100‑1 (MQ=255)
           acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttc                                                                                          >  2:2714955/1‑99 (MQ=255)
            cacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTc                                                                                         >  2:3532369/1‑100 (MQ=255)
            cacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTc                                                                                         >  1:1648317/1‑100 (MQ=255)
             acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCa                                                                                        >  1:2076520/1‑100 (MQ=255)
             acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCa                                                                                        <  2:600224/100‑1 (MQ=255)
              ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                       <  1:3843824/100‑1 (MQ=38)
              ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                       <  1:4101392/100‑1 (MQ=38)
              ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGAGCAAGATCATAAAAATTGTGATGTAAGTCACGATTTTAAt                                                                                       >  2:1749998/1‑100 (MQ=255)
               cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATc                                                                                      <  1:2667558/100‑1 (MQ=255)
               cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                      <  1:1582944/100‑2 (MQ=37)
               cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                      <  1:1403488/100‑2 (MQ=38)
                aGAAAGATAAAAATATAAAATACATACATAAAGACATGACATCAATGTACGGGTCAGTGCAATATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                       <  2:2069388/98‑1 (MQ=16)
                aGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                      <  2:1060540/99‑2 (MQ=37)
                 gAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCtt                                                                                    <  1:3292702/100‑1 (MQ=255)
                 gAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtgat                                                                                    <  2:1658387/100‑4 (MQ=37)
                 gAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTCATTTACATCACAATTTTCAtgat                                                                                    <  1:3553475/100‑4 (MQ=37)
                   aaGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTg                                                                                  >  2:3805412/1‑100 (MQ=255)
                     gagaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCt                                                                                <  1:3363654/100‑1 (MQ=255)
                      agaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCtt                                                                               >  1:106422/1‑100 (MQ=255)
                      agaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCAGCTTTGCt                                                                                >  1:1311121/1‑99 (MQ=255)
                        aacaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTa                                                                             <  1:3490828/97‑1 (MQ=255)
                        aaaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTa                                                                             <  2:407686/100‑1 (MQ=255)
                        aaaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTa                                                                             <  1:4008063/100‑1 (MQ=255)
                         aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTaa                                                                            >  2:2457678/1‑100 (MQ=255)
                         aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTaa                                                                            <  1:778417/100‑1 (MQ=255)
                         aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTaa                                                                            >  2:1744065/1‑100 (MQ=255)
                          aaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAAc                                                                           <  1:695718/100‑1 (MQ=255)
                             tataAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGcc                                                                        >  1:1868204/1‑100 (MQ=255)
                               taAAATACATACAGAAAGACATGAGAGCAAGGTACGGGTCAGTGCAAGAGCATTAAAATTGTGATGTAAATCACGATGTTCATCTTTGCTGTAACGCCTa                                                                      >  2:3500509/1‑100 (MQ=255)
                                  cgatcttacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAg                                                                   <  1:1066330/94‑1 (MQ=255)
                                   atacatacAGAAAGGCATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAAGTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                  <  1:996744/100‑1 (MQ=255)
                                   atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                  <  1:1635929/100‑1 (MQ=255)
                                    tacatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                                 >  1:28749/1‑100 (MQ=255)
                                    tacatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGt                                                                 <  1:2338271/100‑1 (MQ=255)
                                     acatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTg                                                                <  2:1406236/100‑1 (MQ=255)
                                      catacaGAACGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGCTGa                                                               >  2:193530/1‑100 (MQ=255)
                                      catacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGa                                                               >  2:3108831/1‑100 (MQ=255)
                                      catacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGa                                                               >  1:4195624/1‑100 (MQ=255)
                                      catacaGAAAGACATGAGAGCAAGGGACGGAGCGGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGACTACAGGTGa                                                               >  2:2763783/1‑100 (MQ=255)
                                       atacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGAt                                                              <  2:3915235/100‑1 (MQ=255)
                                          caGAAAGACATGAGAGCAAGGGACGGGGCAGAGCAAGATCATAAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaa                                                           >  2:3892305/1‑100 (MQ=255)
                                             aaaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAt                                                        <  2:1041975/100‑1 (MQ=255)
                                              aaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAtt                                                       >  1:1446497/1‑100 (MQ=255)
                                               aGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAttt                                                      >  2:1088810/1‑100 (MQ=255)
                                                gACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTa                                                     <  1:3766018/100‑1 (MQ=255)
                                                  cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                   >  1:2733917/1‑100 (MQ=255)
                                                     gagaGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCtt                                                <  2:4068179/100‑1 (MQ=255)
                                                     gagaGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCtt                                                >  2:1239558/1‑100 (MQ=255)
                                                       gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                              >  2:3390815/1‑100 (MQ=255)
                                                       gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                              >  1:1114692/1‑100 (MQ=255)
                                                         gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGa                                            >  1:2092006/1‑100 (MQ=255)
                                                         gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGa                                            <  1:3862686/100‑1 (MQ=255)
                                                          cAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaa                                           >  1:1899115/1‑100 (MQ=255)
                                                           aaGGGACGGGGCAGGGCAAGATCATTAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAAt                                          <  2:3549605/100‑1 (MQ=255)
                                                              ggACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAAc                                       <  1:1701766/100‑1 (MQ=255)
                                                               gACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAGTCACGATTTTCGTCTTTGCTTTAACGACTACAGGTGATCAAAATTTAACCTTTGCAATAAc                                       >  2:2234005/1‑99 (MQ=255)
                                                                aCGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCAAGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAAc                                       <  2:906106/98‑1 (MQ=255)
                                                                  ggggCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaa                                   >  2:2978443/1‑100 (MQ=255)
                                                                   gggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTTGAATAACTaaaa                                  <  1:2763790/100‑1 (MQ=255)
                                                                   gggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaa                                  <  1:1088810/100‑1 (MQ=255)
                                                                    ggcagagcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                                 <  2:3172220/100‑1 (MQ=255)
                                                                         ggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGTTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaa                            >  2:3288136/1‑100 (MQ=255)
                                                                         ggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaa                            <  1:618851/100‑1 (MQ=255)
                                                                           caAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                          >  2:1332084/1‑100 (MQ=255)
                                                                             caaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaac                        >  2:1844162/6‑99 (MQ=25)
                                                                             caaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa                        >  1:383376/6‑100 (MQ=25)
                                                                             caaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa                        >  1:2171985/6‑100 (MQ=25)
                                                                             aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa                        >  1:2978603/1‑100 (MQ=255)
                                                                              gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg                       <  1:2975251/100‑1 (MQ=255)
                                                                                 gaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc                    >  1:280951/2‑100 (MQ=37)
                                                                                 cATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc                    >  1:1950021/1‑100 (MQ=255)
                                                                                  aTGCAAATTGTGATGTAAATCACGCTTTTCATCTTTGCTTTCACGCCTAAAGTAGATCAAGAATTACCCCACGGCATAACCAACAAGCTCCACAAGGAc                    >  2:3892745/1‑99 (MQ=255)
                                                                                   tGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGc                  >  1:3505606/1‑100 (MQ=255)
                                                                                   tGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGc                  >  2:657544/1‑100 (MQ=255)
                                                                                     aaaaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCa                >  1:3656502/1‑100 (MQ=255)
                                                                                         ttGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTg            <  2:2844409/100‑1 (MQ=255)
                                                                                           gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGa           >  1:2510117/1‑99 (MQ=255)
                                                                                           gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg          <  2:1071462/100‑1 (MQ=255)
                                                                                           gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg          <  2:4107338/100‑1 (MQ=255)
                                                                                             gATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAgg         >  2:2543126/1‑98 (MQ=255)
                                                                                             gATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTa        >  1:1870635/1‑100 (MQ=255)
                                                                                             gATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTa        <  2:4189402/100‑1 (MQ=255)
                                                                                              aTGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTAt       <  1:3785225/100‑1 (MQ=255)
                                                                                                  aaaTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCaa   >  1:3662299/1‑100 (MQ=255)
                                                                                                   aaTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCaan  <  2:229732/100‑2 (MQ=255)
                                                                                                   aaTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCaaa  <  1:2287256/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAA  >  NC_000913/1530400‑1530598

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