Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,490 T→C 15.6% intergenic (+86/‑96) 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,4900TC15.6% 399.8 / 61.5 167intergenic (+86/‑96)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base T (74/67);  new base C (12/14);  total (86/81)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-01

Only 100 of 177 total aligned reads displayed.

TAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTG  >  NC_000913/1530392‑1530588
                                                                                                  |                                                                                                  
tAAGTGGTCATAATGCCAAACAACAGAAAGAGAAAAATATCAAATACATACAGAAAGACATGAGAGCAAGAGACGGGGCAGGGCAAGATCATGAAAATtg                                                                                                   >  2:1137837/1‑100 (MQ=255)
tAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATtg                                                                                                   <  1:2888275/100‑1 (MQ=255)
  aaTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATtgtg                                                                                                 <  1:2608595/100‑1 (MQ=255)
  aaTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATtgtg                                                                                                 <  1:203923/100‑1 (MQ=255)
  aaTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATAAAGAAAGACATGAGAGAAAGGGACGAGGCAGGGCAAGATCAAGAAAAAtg                                                                                                   >  2:1471431/1‑98 (MQ=255)
   aTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGa                                                                                                <  1:2204251/100‑1 (MQ=255)
   aTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGa                                                                                                <  1:1666828/100‑1 (MQ=255)
    tGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGAt                                                                                               <  1:1618182/100‑1 (MQ=255)
     ggTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATg                                                                                              <  2:3235714/100‑1 (MQ=255)
      gTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGt                                                                                             <  1:1057334/100‑1 (MQ=255)
      gTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGt                                                                                             >  1:445272/1‑100 (MQ=255)
      gTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGt                                                                                             >  1:1271032/1‑100 (MQ=255)
       tCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTa                                                                                            <  1:674379/100‑1 (MQ=255)
       ctaTAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTa                                                                                            <  1:2353741/98‑1 (MQ=255)
        cATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTaa                                                                                           >  1:546262/1‑100 (MQ=255)
         aTAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTaag                                                                                          >  2:2353736/1‑99 (MQ=255)
         aTAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTaaa                                                                                          >  1:835941/1‑100 (MQ=255)
          tAATGAACGACACCAGAAAGAGAAAAATATGAAAAATATACAGAAAGACATGAGAGAAAAGTACAGGGCAGGTCACGATCATGAACATTGTGCTCTAAAt                                                                                         >  2:1074041/1‑100 (MQ=255)
          tAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAt                                                                                         >  1:3584494/1‑100 (MQ=255)
          tAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAt                                                                                         >  1:1071676/1‑100 (MQ=255)
            aTGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATTTACATCa                                                                                       <  2:2192762/100‑1 (MQ=255)
             tGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCAc                                                                                      >  1:1438691/1‑100 (MQ=255)
               aacaacACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGa                                                                                    >  1:2334835/1‑100 (MQ=255)
                acaacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAt                                                                                   <  1:1825973/100‑1 (MQ=255)
                  aacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAttt                                                                                 <  2:931159/100‑1 (MQ=255)
                   acacCAGAAAGAGAAAAATATAAACAACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                <  2:399078/100‑1 (MQ=255)
                   acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                >  2:3477279/1‑100 (MQ=255)
                   acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttc                                                                                >  2:2714955/1‑99 (MQ=255)
                    cacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTc                                                                               >  1:1648317/1‑100 (MQ=255)
                    cacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTc                                                                               >  2:3532369/1‑100 (MQ=255)
                     acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCa                                                                              <  2:600224/100‑1 (MQ=255)
                     acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCa                                                                              >  1:2076520/1‑100 (MQ=255)
                      ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCAt                                                                             >  2:1363858/1‑100 (MQ=255)
                      ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                             <  1:4101392/100‑1 (MQ=38)
                      ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                             <  1:3843824/100‑1 (MQ=38)
                      ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGAGCAAGATCATAAAAATTGTGATGTAAGTCACGATTTTAAt                                                                             >  2:1749998/1‑100 (MQ=255)
                      ccAGAAAGAGAAAAAGATAAAATACATACAGAAAGGCATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAAGTGTGATGTAAATCACGATGTTCa                                                                              >  2:839161/1‑99 (MQ=255)
                       cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCGAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                            <  1:730755/100‑2 (MQ=37)
                       cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATc                                                                            >  2:2260246/1‑100 (MQ=255)
                       cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATc                                                                            <  1:2667558/100‑1 (MQ=255)
                        aGAAAGATAAAAATATAAAATACATACATAAAGACATGACATCAATGTACGGGTCAGTGCAATATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                             <  2:2069388/98‑1 (MQ=16)
                        aGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                            <  2:1060540/99‑2 (MQ=37)
                         gAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCtt                                                                          <  1:3292702/100‑1 (MQ=255)
                          aaaGAGAAAAATATAAAATACATACAGAAAGCCATGACAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCttt                                                                         <  2:2822030/100‑1 (MQ=255)
                           aaGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTg                                                                        <  2:4105645/100‑1 (MQ=255)
                            agagAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGc                                                                       <  2:3901409/100‑1 (MQ=255)
                             gagaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCt                                                                      <  1:3363654/100‑1 (MQ=255)
                              agaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCtt                                                                     >  1:106422/1‑100 (MQ=255)
                              agaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCAGCTTTGCt                                                                      >  1:1311121/1‑99 (MQ=255)
                               gaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCttt                                                                    <  1:371271/100‑1 (MQ=255)
                                aacaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTa                                                                   <  1:3490828/97‑1 (MQ=255)
                                aaaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTa                                                                   <  2:407686/100‑1 (MQ=255)
                                 aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTaa                                                                  <  1:1111409/100‑1 (MQ=255)
                                 aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTaa                                                                  <  1:778417/100‑1 (MQ=255)
                                 aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTaa                                                                  >  2:1725961/1‑100 (MQ=255)
                                 aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTaa                                                                  >  2:1744065/1‑100 (MQ=255)
                                    atatAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGc                                                               <  2:3538292/100‑1 (MQ=255)
                                     tataAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGcc                                                              >  1:1868204/1‑100 (MQ=255)
                                     tataAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGcc                                                              <  1:2906226/100‑1 (MQ=255)
                                       taAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTa                                                            <  2:1045742/100‑1 (MQ=255)
                                         aaaTACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACa                                                          <  2:4006694/100‑1 (MQ=255)
                                           atacatacAGAAAGGCATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAAGTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                        <  1:996744/100‑1 (MQ=255)
                                           atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                        <  2:1799755/100‑1 (MQ=255)
                                              catacaGAACGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGCTGa                                                     >  2:193530/1‑100 (MQ=255)
                                              catacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGa                                                     >  1:4195624/1‑100 (MQ=255)
                                              catacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGa                                                     >  2:518049/1‑100 (MQ=255)
                                                 acaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                   <  2:3616071/99‑1 (MQ=255)
                                                  caGAAAGACATGAGAGCAAGGGACGGGGCAGAGCAAGATCATAAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaa                                                 >  2:3892305/1‑100 (MQ=255)
                                                   aGAAAGACATGCGTGCAAGGGACGGGTCAGGTCAAGATCATGAAAATTGTGATGTAAATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCa                                                  <  2:3079991/98‑1 (MQ=255)
                                                        gACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTa                                           <  1:3766018/100‑1 (MQ=255)
                                                             gagaGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCtt                                      >  2:1239558/1‑100 (MQ=255)
                                                             gagaGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCtt                                      <  1:3805436/100‑1 (MQ=255)
                                                               gagCAAGGGACGGGTCAGGGCAAGATCAGGAAAATTGTGATGTAAAGCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                    <  1:2681329/100‑1 (MQ=255)
                                                                 gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGa                                  <  1:3862686/100‑1 (MQ=255)
                                                                  cAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaa                                 >  1:1899115/1‑100 (MQ=255)
                                                                  cAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGa                                  >  1:1601130/1‑99 (MQ=255)
                                                                      ggACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAAc                             <  1:1701766/100‑1 (MQ=255)
                                                                      ggACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAAc                             >  1:4060410/1‑100 (MQ=255)
                                                                       gACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAGTCACGATTTTCGTCTTTGCTTTAACGACTACAGGTGATCAAAATTTAACCTTTGCAATAAc                             >  2:2234005/1‑99 (MQ=255)
                                                                       gACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAAc                             <  2:908215/99‑1 (MQ=255)
                                                                           gggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTTGAATAACTaaaa                        <  1:2763790/100‑1 (MQ=255)
                                                                            ggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                       <  1:1506512/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)
                                                                                 ggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaa                  >  1:3955930/1‑100 (MQ=255)
                                                                                 aaagcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaa                  >  2:1247096/10‑100 (MQ=21)
                                                                                  aagcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTACATTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATcaa                 >  2:1932978/9‑100 (MQ=18)
                                                                                   agcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                >  2:1549234/8‑100 (MQ=25)
                                                                                   agcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                >  1:3576321/8‑100 (MQ=25)
                                                                                     caaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa              >  1:3513781/6‑100 (MQ=25)
                                                                                     aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa              >  1:2978603/1‑100 (MQ=255)
                                                                                     aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa              >  1:2178311/1‑100 (MQ=255)
                                                                                      gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg             <  2:337971/100‑1 (MQ=255)
                                                                                      aaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg             >  1:2147203/5‑100 (MQ=37)
                                                                                      aaagaTGAAAATCGTGATTTACATCACAATTTTAATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg             >  2:2510045/5‑100 (MQ=25)
                                                                                         gaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc          >  1:280951/2‑100 (MQ=37)
                                                                                         cATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc          >  1:1950021/1‑100 (MQ=255)
                                                                                                 ttGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTg  <  1:1795616/100‑1 (MQ=255)
                                                                                                 ttGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTg  <  2:2844409/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
TAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTG  >  NC_000913/1530392‑1530588

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: