Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,643,776 C→G 5.4% D177H (GAT→CAT)  ydfU ← Qin prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,643,7760CG 5.4% 398.7 / 10.8 112D177H (GAT→CAT) ydfUQin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (54/52);  new base G (3/3);  total (57/55)
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.47e-01

Only 100 of 117 total aligned reads displayed.

GCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTG  >  NC_000913/1643637‑1643923
                                                                                                                                           |                                                                                                                                                   
gCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCt                                                                                                                                             >  3:735623/1‑148 (MQ=255)
   aGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGAt                                                                                                                                          <  4:735623/148‑1 (MQ=255)
                 ttgttTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAg                                                                                                                                     <  4:883587/139‑1 (MQ=255)
                 ttgttTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAg                                                                                                                                     >  3:883587/1‑139 (MQ=255)
                     ttCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAgg                                                                                                                        <  2:79008/148‑1 (MQ=255)
                            tGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGc                                                                                                                                       >  2:972457/1‑126 (MQ=255)
                            tGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGc                                                                                                                                       <  1:972457/126‑1 (MQ=255)
                             gCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATccacca                                                                                                                                                <  2:932248/116‑1 (MQ=255)
                             gCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATccacca                                                                                                                                                >  1:932248/1‑116 (MQ=255)
                                        tCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATaa                                                                                                     >  1:620799/1‑148 (MQ=255)
                                          aGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAAtt                                                                                                   <  2:620799/148‑1 (MQ=255)
                                           gCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTc                                                                                                  >  3:625378/1‑148 (MQ=255)
                                                               cTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTTCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTa                                                                              <  2:762145/148‑1 (MQ=255)
                                                               cTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCt                                                                                          <  3:1043825/136‑1 (MQ=255)
                                                               cTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCt                                                                                          >  4:1043825/1‑136 (MQ=255)
                                                                 ttCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCtt                                                                               >  1:865356/1‑145 (MQ=255)
                                                                 ttCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCtt                                                                               <  2:865356/145‑1 (MQ=255)
                                                                    gcgGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCg                                                                         <  4:625378/148‑1 (MQ=255)
                                                                            ccGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAg                                                                                                               <  1:452624/102‑1 (MQ=255)
                                                                            ccGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAg                                                                                                               >  2:452624/1‑102 (MQ=255)
                                                                            ccGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTcc                                                                                                 >  1:703553/1‑116 (MQ=255)
                                                                              gAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCg                                                                                               >  1:7347/1‑116 (MQ=255)
                                                                              gAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCg                                                                                               <  2:7347/116‑1 (MQ=255)
                                                                               aGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGAc                                                              >  2:785474/1‑148 (MQ=255)
                                                                                   gATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCa                                                          >  2:837924/1‑148 (MQ=255)
                                                                                     ttttttCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACgg                                                                   >  1:824774/1‑137 (MQ=255)
                                                                                     ttttttCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACgg                                                                   <  2:824774/137‑1 (MQ=255)
                                                                                     ttttttCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAAt                                                        >  3:881645/1‑148 (MQ=255)
                                                                                          tCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCaa                                                                                                           >  4:1059801/1‑92 (MQ=255)
                                                                                          tCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCaa                                                                                                           <  3:1059801/92‑1 (MQ=255)
                                                                                            cGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCAt                                                 <  4:881645/148‑1 (MQ=255)
                                                                                                  gattaccAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATCGCGTGACGGATGACGTCAATCACCCATGCTGTc                                           >  3:675079/6‑148 (MQ=255)
                                                                                                   aaTCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAAGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCa                                          >  4:175633/1‑148 (MQ=255)
                                                                                                      cccAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCaa                                                                                                           <  1:168825/80‑1 (MQ=255)
                                                                                                      cccAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCaa                                                                                                           >  2:168825/1‑80 (MQ=255)
                                                                                                        cAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTc                                                           >  4:238402/1‑126 (MQ=255)
                                                                                                        cAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTc                                                           <  3:238402/126‑1 (MQ=255)
                                                                                                        cAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATCGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAttt                                     <  4:675079/148‑1 (MQ=255)
                                                                                                         aGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCa                                                                                                                <  3:727821/72‑1 (MQ=255)
                                                                                                         aGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCa                                                                                                                >  4:727821/1‑72 (MQ=255)
                                                                                                         aGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAtttt                                    >  1:161141/1‑148 (MQ=255)
                                                                                                         aGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAtttt                                    <  1:785474/148‑1 (MQ=255)
                                                                                                         aGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAtttt                                    >  4:656520/1‑148 (MQ=255)
                                                                                                          gCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTg                                                                           <  3:375946/108‑1 (MQ=255)
                                                                                                          gCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTg                                                                           >  4:375946/1‑108 (MQ=255)
                                                                                                           cGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTa                                                                              <  4:1035627/104‑1 (MQ=255)
                                                                                                           cGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgt                                  <  3:175633/148‑1 (MQ=255)
                                                                                                           cGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgt                                  <  3:656520/148‑1 (MQ=255)
                                                                                                               cgacgCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgt                                  >  2:111357/1‑144 (MQ=255)
                                                                                                               cgacgCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgt                                  <  1:111357/144‑1 (MQ=255)
                                                                                                                 acgCGATACAGCCTCAGGTAGTTCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAtt                                      >  4:70107/1‑138 (MQ=255)
                                                                                                                 acgCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGGTGTCAGAtt                                      <  3:70107/138‑1 (MQ=255)
                                                                                                                    cgATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTg                                                                       <  2:928094/102‑1 (MQ=255)
                                                                                                                    cgATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGt                                            >  4:136100/1‑129 (MQ=255)
                                                                                                                       tACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGa                                                                      <  4:528511/100‑1 (MQ=255)
                                                                                                                       tACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGa                                                                      >  3:528511/1‑100 (MQ=255)
                                                                                                                       tACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAt                                       >  1:599881/1‑131 (MQ=255)
                                                                                                                       tACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAt                                       <  2:599881/131‑1 (MQ=255)
                                                                                                                        aCAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCt                     >  3:370887/1‑148 (MQ=255)
                                                                                                                         cAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGa                                                                  >  2:137593/1‑102 (MQ=255)
                                                                                                                         cAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGa                                                                  <  1:137593/102‑1 (MQ=255)
                                                                                                                         cAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGa                                                               >  2:558593/1‑105 (MQ=255)
                                                                                                                         cAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGa                                                               <  1:927295/105‑1 (MQ=255)
                                                                                                                         cAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGa                                                               >  2:927295/1‑105 (MQ=255)
                                                                                                                         cAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGa                                                               <  1:558593/105‑1 (MQ=255)
                                                                                                                           gCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCa                  <  1:555654/148‑1 (MQ=255)
                                                                                                                             cTCAGTTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgttg                                <  1:399490/132‑1 (MQ=255)
                                                                                                                             cTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCCCCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGa                                                                  <  2:222833/98‑1 (MQ=255)
                                                                                                                             cTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgttg                                >  2:399490/1‑132 (MQ=255)
                                                                                                                             atCAGGTGGTGGATGCACTACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGt                >  1:526870/2‑148 (MQ=255)
                                                                                                                              tCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTg               <  1:837924/148‑1 (MQ=255)
                                                                                                                               cAGGTGGTGGATGCACTACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCtt                                                                               >  2:368192/1‑83 (MQ=38)
                                                                                                                               cAGGTGGTGGATGCACTACCTGATTGCAGACCGCCCACCAGGATAATTCAGCAAAAGATAATTCCCGCTCCTGCGTACCGCtt                                                                               <  1:368192/83‑1 (MQ=38)
                                                                                                                               cAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttg                                   >  4:582579/1‑127 (MQ=255)
                                                                                                                               cAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttg                                   <  2:136006/127‑1 (MQ=255)
                                                                                                                               cAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttg                                   <  3:582579/127‑1 (MQ=255)
                                                                                                                               cAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttg                                   >  1:136006/1‑127 (MQ=255)
                                                                                                                               cAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttg                                   >  3:1131826/1‑127 (MQ=255)
                                                                                                                                  gtggTGGATGCACTACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCaa                                                         <  2:643658/99‑1 (MQ=255)
                                                                                                                                  gtggTGGATGCACTACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAtttt                                    <  2:570229/120‑1 (MQ=255)
                                                                                                                                  gtggTGGATGCACTACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGAtttt                                    >  1:570229/4‑123 (MQ=255)
                                                                                                                                  gTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGGCGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTc           <  3:515531/148‑1 (MQ=255)
                                                                                                                                  gTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGGCGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTc           >  4:515531/1‑148 (MQ=255)
                                                                                                                                   tAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCAcc                                                    >  4:1111566/1‑106 (MQ=255)
                                                                                                                                   tAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCAcc                                                    <  3:1111566/106‑1 (MQ=255)
                                                                                                                                   tAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGt                                            <  2:117695/114‑1 (MQ=255)
                                                                                                                                   tAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGt                                            >  1:117695/1‑114 (MQ=255)
                                                                                                                                      tGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTAt                                                                             <  4:1072916/78‑1 (MQ=255)
                                                                                                                                      tGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTAt                                                                             >  3:1072916/1‑78 (MQ=255)
                                                                                                                                        cATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgt                                  <  3:884696/119‑1 (MQ=255)
                                                                                                                                         aTCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgt                                  >  4:61885/1‑118 (MQ=255)
                                                                                                                                         aTCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgt                                  >  4:1048711/1‑118 (MQ=255)
                                                                                                                                         aTCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgt                                  >  2:557860/1‑118 (MQ=255)
                                                                                                                                         aTCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgt                                  <  3:61885/118‑1 (MQ=255)
                                                                                                                                         aTCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTttgt                                  <  3:1048711/118‑1 (MQ=255)
                                                                                                                                          tCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATc                                                       <  1:399533/96‑1 (MQ=255)
                                                                                                                                          tCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATc                                                       >  2:399533/1‑96 (MQ=255)
                                                                                                                                          tCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCAccc                                                   >  4:576257/1‑100 (MQ=255)
                                                                                                                                          tCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCAccc                                                   <  3:576257/100‑1 (MQ=255)
                                                                                                                                           ccaccaCCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGt                                                            >  1:847123/1‑90 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATCACCCATGCTGTCAGATTTTGTTGAGCAAGCAGCTCCAGTGATTCCGATGTCTG  >  NC_000913/1643637‑1643923

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