Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,499 A→C 16.5% 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.5% 345.9 / 60.1 138intergenic (+95/‑87)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base A (52/63);  new base C (9/14);  total (61/77)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 140 total aligned reads displayed.

CATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAA  >  NC_000913/1530400‑1530597
                                                                                                   |                                                                                                  
cATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTaa                                                                                                    >  2:2362165/1‑100 (MQ=255)
  tAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGATCAATGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAt                                                                                                  <  1:1921056/100‑1 (MQ=255)
  tAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAt                                                                                                  >  2:2023785/1‑100 (MQ=255)
  tAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAt                                                                                                  <  2:3041100/100‑1 (MQ=255)
   aaTGATCAACACCAAAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGATCAAGGGACGGTGCAGGGCAAGATCATGAAAATTGTGATGTAAATc                                                                                                 <  1:2633830/100‑1 (MQ=255)
         caacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtt                                                                                           >  2:2816216/1‑100 (MQ=255)
           acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                         >  1:3731149/1‑100 (MQ=255)
           acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                         <  1:2539721/100‑1 (MQ=255)
           acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAAtttt                                                                                         <  2:591824/100‑1 (MQ=38)
           acacCAGAAAGAGAAAAATATAAAAGACATACAGAAAGACATGAAAGCATGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                         <  1:3610815/100‑1 (MQ=255)
           acaacaGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                         <  1:1200495/100‑1 (MQ=255)
             acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                                       <  2:3600043/100‑1 (MQ=38)
             acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                                       <  2:1632746/100‑1 (MQ=38)
             acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                                       <  2:774146/100‑1 (MQ=38)
             acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                                       <  2:3704836/100‑1 (MQ=38)
             acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                                       <  2:308728/100‑1 (MQ=38)
              ccAGAAAGAGAAAAATATAAAATACATGCAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                      <  2:3081755/100‑1 (MQ=37)
              ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCAt                                                                                      >  2:3407229/1‑100 (MQ=255)
              ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                      <  1:427668/100‑1 (MQ=38)
              ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                      <  2:1176875/100‑1 (MQ=38)
               cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                     <  1:2364434/100‑2 (MQ=37)
               cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                     <  2:2024741/100‑2 (MQ=37)
                aGAAAGAGTAAAATATAAAATACATACAGAAAGACATGAGAGCAAGAGACGGGGCAGGGCAAGATCATGAAATTCGTGATTTACATCACAATTTTCAt                                                                                      <  2:1377457/98‑1 (MQ=37)
                  aaaGAGAAGATTATAAAACACATATACAAAGACATGAGAGCAAGATACAGATCAGGTTAAGCTCAACAAAATAGTAATTTACAACTCAATTTTCAtg                                                                                     <  1:1792911/97‑2 (MQ=255)
                  aaaGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGTACGGGGCAGGGCAAGATTATAAAAATCGTGATTTACATCACAATTTTCAtgat                                                                                   <  1:1158111/99‑4 (MQ=25)
                   aaGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTg                                                                                 <  2:1496519/100‑1 (MQ=255)
                     gagaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCt                                                                               >  1:2096847/1‑100 (MQ=255)
                      agaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCtt                                                                              <  2:3538954/100‑1 (MQ=255)
                      agaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCtt                                                                              <  2:2763712/100‑1 (MQ=255)
                       gaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCttt                                                                             >  2:2376982/1‑100 (MQ=255)
                       gaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGAAGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATATTTGCttt                                                                             >  1:2278334/1‑100 (MQ=255)
                          aaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAAc                                                                          <  1:460506/100‑1 (MQ=255)
                           aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                                         <  1:820567/100‑1 (MQ=255)
                           aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                                         <  1:819521/100‑1 (MQ=255)
                           aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                                         <  1:800251/100‑1 (MQ=255)
                                   atacatacAGAAAGACATGAGAGCAAGGGCCGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                 <  2:759769/100‑1 (MQ=255)
                                   atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgn                                                                 <  2:1461873/100‑2 (MQ=255)
                                   atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                 <  2:145528/100‑1 (MQ=255)
                                   atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                 >  2:3271242/1‑100 (MQ=255)
                                   atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                 <  2:2635560/100‑1 (MQ=255)
                                      catacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGa                                                              >  2:2519413/1‑100 (MQ=255)
                                       atacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAAGTGTGATGTAAATCACGATTTTCATCTTTGCTTTAAAGCCTAAAg                                                                  >  1:37431/1‑95 (MQ=255)
                                        tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                            <  2:1311441/100‑1 (MQ=255)
                                        tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                            >  1:1475638/1‑100 (MQ=255)
                                        tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                            >  2:742730/1‑100 (MQ=255)
                                        tacaGAAAGACATGAGAGCAAGGGACGGGGAAGGGCAAGATCATGAAAATTATGATGTAAATCACGATTTTCATCTTTGCTTTAACGCATACAGGTg                                                               >  1:123532/1‑97 (MQ=255)
                                          caGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaa                                                          >  2:3280529/1‑100 (MQ=255)
                                           aGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaa                                                         >  1:1378671/1‑100 (MQ=255)
                                               aGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAttt                                                     <  1:3040660/100‑1 (MQ=255)
                                                gACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTa                                                    <  2:3788447/100‑1 (MQ=255)
                                                  cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                  <  1:1075710/100‑1 (MQ=255)
                                                  cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                  <  2:1746661/100‑1 (MQ=255)
                                                  cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                  >  2:3595869/1‑100 (MQ=255)
                                                   aTGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAcc                                                 <  1:3184157/100‑1 (MQ=255)
                                                    tGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                                >  1:3362433/1‑100 (MQ=255)
                                                      agagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCttt                                              >  1:3642527/1‑100 (MQ=255)
                                                       gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                             <  2:3380170/100‑1 (MQ=255)
                                                       gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                             >  2:392389/1‑100 (MQ=255)
                                                       gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                             >  2:2181803/1‑100 (MQ=255)
                                                        agCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTgg                                            <  2:105579/100‑1 (MQ=255)
                                                         gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGa                                           >  2:1731911/1‑100 (MQ=255)
                                                          cAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaa                                          <  2:2640824/100‑1 (MQ=255)
                                                          cAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaa                                          >  1:1925510/1‑100 (MQ=255)
                                                           aaGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAAt                                         <  1:3385376/100‑1 (MQ=255)
                                                           aaGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAAt                                         >  1:1836759/1‑100 (MQ=255)
                                                            aGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATa                                        <  1:1447479/100‑1 (MQ=255)
                                                              ggACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAGTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAAc                                      >  1:1108671/1‑100 (MQ=255)
                                                               gACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACt                                     >  1:1598609/1‑100 (MQ=255)
                                                               gACGCGGCAGGGCAAGATAATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGTTGATCAAAATTTAACCTTTGGATTAACt                                     >  1:37990/1‑100 (MQ=255)
                                                                   gggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaa                                 >  2:2500822/1‑100 (MQ=255)
                                                                    ggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                                >  2:2170580/1‑100 (MQ=255)
                                                                        nggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                            >  1:3007427/2‑100 (MQ=255)
                                                                        gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                            <  2:1549271/100‑1 (MQ=255)
                                                                        gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                            <  2:1088404/100‑1 (MQ=255)
                                                                         ggcaAGATGATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGATTTAACGCCTACCGCTAATCAAAATTTAACCTTTGGAATAACTAAAAAGATaa                           >  1:278311/1‑100 (MQ=255)
                                                                         aaagcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaa                           >  1:3027808/10‑100 (MQ=21)
                                                                           ngcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                         >  2:1944319/8‑100 (MQ=25)
                                                                           caAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                         <  1:3651172/100‑1 (MQ=255)
                                                                              gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg                      <  1:742722/100‑1 (MQ=255)
                                                                              aaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg                      >  1:3220820/5‑100 (MQ=25)
                                                                              aaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg                      >  2:118858/5‑100 (MQ=37)
                                                                                 gaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc                   >  2:1887526/2‑100 (MQ=37)
                                                                                 gaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc                   >  1:706802/2‑100 (MQ=37)
                                                                                  aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAgg                     >  1:2926395/1‑97 (MQ=255)
                                                                                  aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACg                  <  2:1066364/100‑1 (MQ=255)
                                                                                   tGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGc                 <  1:2362154/100‑1 (MQ=255)
                                                                                   tGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGc                 >  2:3800986/1‑100 (MQ=38)
                                                                                     aaaaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCa               <  1:1463338/100‑1 (MQ=255)
                                                                                     aaaaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCa               <  1:309014/100‑1 (MQ=255)
                                                                                       aaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAgg             <  2:937305/100‑1 (MQ=255)
                                                                                       aaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCCTTGGAATAACTAAAAAGATAAAAAAGGACGCCAgg             <  1:2102369/100‑1 (MQ=255)
                                                                                           gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg         <  2:2453678/100‑1 (MQ=255)
                                                                                           gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg         <  2:810123/100‑1 (MQ=255)
                                                                                             gATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTa       >  1:282269/1‑100 (MQ=255)
                                                                                              aTGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACATTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTAt      >  1:2447229/1‑100 (MQ=255)
                                                                                               tGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTAtt     <  2:2626560/100‑1 (MQ=255)
                                                                                                 tAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCa   >  2:2017882/1‑100 (MQ=255)
                                                                                                 tAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCa   >  1:2990837/1‑100 (MQ=255)
                                                                                                 tAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCa   >  1:1745494/1‑100 (MQ=255)
                                                                                                  aaaTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCaa  >  2:1086357/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAA  >  NC_000913/1530400‑1530597

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: