Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,485,103 G→C 5.9% G681G (GGG→GGC)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,485,1030GC 5.9% 370.1 / 12.8 101G681G (GGG→GGC) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (57/38);  new base C (3/3);  total (60/41)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 115 total aligned reads displayed.

CCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGATTGTT  >  NC_000913/1484956‑1485250
                                                                                                                                                   |                                                                                                                                                   
ccATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGgcg                                                                                                                                                   <  2:1042156/150‑1 (MQ=255)
ccATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGgcg                                                                                                                                                   <  2:1502074/150‑1 (MQ=255)
     tCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATtg                                                                                                                                              >  1:1211697/1‑150 (MQ=255)
     tCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATtg                                                                                                                                              >  1:120664/1‑150 (MQ=255)
        gCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATtgctg                                                                                                                                           >  2:212736/1‑150 (MQ=255)
           tGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTgcgc                                                                                                                                        >  1:1275679/1‑150 (MQ=255)
           tGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTgcgc                                                                                                                                        >  2:162594/1‑150 (MQ=255)
             aaaGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGt                                                                                                                                      <  1:510953/150‑1 (MQ=255)
                  tGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGa                                                                                                                                 <  2:264442/150‑1 (MQ=255)
                   gCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGAc                                                                                                                                >  2:244019/1‑150 (MQ=255)
                    ccGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGAcc                                                                                                                               >  2:1288816/1‑150 (MQ=255)
                     cGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGAccc                                                                                                                              <  1:60459/150‑1 (MQ=255)
                      gATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCg                                                                                                                             >  1:769931/1‑150 (MQ=255)
                          aaCAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAAt                                                                                                                         >  1:1009617/1‑150 (MQ=255)
                           aCAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATg                                                                                                                        <  2:204123/150‑1 (MQ=255)
                             aaGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATggg                                                                                                                      >  2:869876/1‑150 (MQ=255)
                             aaGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATggg                                                                                                                      >  1:993430/1‑150 (MQ=255)
                             aaGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATggg                                                                                                                      >  1:585989/1‑150 (MQ=255)
                                   ataCCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGc                                                                                                                <  1:981590/150‑1 (MQ=255)
                                         gcgcACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTg                                                                                                          <  2:1498394/150‑1 (MQ=255)
                                            cACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTc                                                                                                       <  2:642679/150‑1 (MQ=255)
                                            cACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTc                                                                                                       <  2:1093914/150‑1 (MQ=255)
                                            cACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGCCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTc                                                                                                       >  2:1364092/1‑150 (MQ=255)
                                            cACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGCCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTc                                                                                                       >  2:1364090/1‑150 (MQ=255)
                                            cACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGCCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTc                                                                                                       >  2:158191/1‑150 (MQ=255)
                                             aCGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGACTGGGTGGAGCCAGTTGCTCa                                                                                                      >  1:1697769/1‑150 (MQ=255)
                                             aCGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCa                                                                                                      >  1:346081/1‑150 (MQ=255)
                                               gTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGc                                                                                                    <  2:1050372/150‑1 (MQ=255)
                                               gTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGc                                                                                                    <  2:1480788/150‑1 (MQ=255)
                                                     cgcgTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGa                                                                                              >  1:1261634/1‑150 (MQ=255)
                                                         tttCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTaa                                                                                          >  1:402992/1‑150 (MQ=255)
                                                            ctcCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACg                                                                                       >  2:857308/1‑150 (MQ=255)
                                                            ctcCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACg                                                                                       >  2:733367/1‑150 (MQ=255)
                                                               cATCTTCCCCGGTTCTGGTTTATTCAATAAACCGCCTTAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCAc                                                                                    <  2:1711744/150‑1 (MQ=255)
                                                               cATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCAc                                                                                    <  2:1295661/150‑1 (MQ=255)
                                                                 tCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCt                                                                                  >  1:727328/1‑150 (MQ=255)
                                                                   ttCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTAc                                                                                >  2:822523/1‑150 (MQ=255)
                                                                      cccGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCAGTTGATTAAACGCACCTACAGc                                                                             >  2:602302/1‑150 (MQ=255)
                                                                       ccGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCg                                                                            >  2:1445524/1‑150 (MQ=255)
                                                                       ccGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCg                                                                            >  2:1445507/1‑150 (MQ=255)
                                                                       ccGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCg                                                                            >  1:1073303/1‑150 (MQ=255)
                                                                         ggTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACAAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGaa                                                                          >  1:810909/1‑150 (MQ=255)
                                                                           ttCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAAcc                                                                        <  2:335756/150‑1 (MQ=255)
                                                                            tCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCg                                                                       >  1:156333/1‑150 (MQ=255)
                                                                            tCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATGTGATTAAACGCACCTACAGCGAACCg                                                                       >  2:260752/1‑150 (MQ=255)
                                                                                 tttATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGACTGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTg                                                                  >  2:1271187/1‑150 (MQ=255)
                                                                                 tttATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTg                                                                  >  2:1411897/1‑150 (MQ=255)
                                                                                 tttATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTg                                                                  >  2:1551597/1‑150 (MQ=255)
                                                                                 tttATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTg                                                                  <  1:22440/150‑1 (MQ=255)
                                                                                     ttCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGaa                                                              >  1:125419/1‑150 (MQ=255)
                                                                                     ttCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGaa                                                              >  1:181274/1‑150 (MQ=255)
                                                                                      tCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAAc                                                             >  1:313854/1‑150 (MQ=255)
                                                                                        aaaaaaCCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACgg                                                           >  1:40909/1‑150 (MQ=255)
                                                                                         aaaaaCCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACTGCGAACCGCACTGGGAACggg                                                          <  1:209539/150‑1 (MQ=255)
                                                                                         aaaaaCCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACggg                                                          <  1:468384/150‑1 (MQ=255)
                                                                                         aaaaaCCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACggg                                                          <  1:468284/150‑1 (MQ=255)
                                                                                            aaCCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGgcgc                                                       >  2:1379164/1‑150 (MQ=255)
                                                                                                gCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAggg                                                   <  1:481256/150‑1 (MQ=255)
                                                                                                  cTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGgcg                                                 >  2:224374/1‑150 (MQ=255)
                                                                                                      ctGGGTAATGGTGGCGGAACTGGTAGAACCCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGAGTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGt                                             <  2:283041/149‑1 (MQ=255)
                                                                                                      ctGGGGAATGGTGGCGGAACTGGTAGAACCCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGt                                             <  2:283037/149‑1 (MQ=255)
                                                                                                      aTGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGt                                             >  1:1562051/1‑150 (MQ=255)
                                                                                                      aTGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGt                                             <  2:756415/150‑1 (MQ=255)
                                                                                                       tGGGTAATGGTGGTGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTg                                            >  1:789367/1‑150 (MQ=255)
                                                                                                         ggTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCTGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGAt                                          <  2:324065/150‑1 (MQ=255)
                                                                                                               ggtggCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAAc                                    >  1:1506850/1‑150 (MQ=255)
                                                                                                                gtggCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACg                                   >  1:111738/1‑150 (MQ=255)
                                                                                                                  ggCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGGATTTGATTTAACGCACCTACAGCGAACCGCTCTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGa                                 >  2:167873/1‑150 (MQ=255)
                                                                                                                  ggCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGa                                 >  2:180457/1‑150 (MQ=255)
                                                                                                                  ggCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGa                                 >  2:180388/1‑150 (MQ=255)
                                                                                                                      gAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGaaaaa                             <  1:1156402/150‑1 (MQ=255)
                                                                                                                      gAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGaaaaa                             <  1:635420/150‑1 (MQ=255)
                                                                                                                       aaCTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAg                            <  1:212736/150‑1 (MQ=255)
                                                                                                                          tGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCACCGGAAAAAGTCa                         >  1:1106472/1‑150 (MQ=255)
                                                                                                                          tGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCa                         >  2:1615496/1‑150 (MQ=255)
                                                                                                                          tGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCa                         >  2:1608136/1‑150 (MQ=255)
                                                                                                                           ggTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCAc                        <  2:165881/150‑1 (MQ=255)
                                                                                                                                 aaCCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTa                  <  1:333429/150‑1 (MQ=255)
                                                                                                                                 aaCCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTa                  <  2:51864/150‑1 (MQ=255)
                                                                                                                                  aCCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTAt                 <  1:1288816/150‑1 (MQ=255)
                                                                                                                                   ccAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATg                >  1:569032/1‑150 (MQ=255)
                                                                                                                                   ccAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATg                >  2:805121/1‑150 (MQ=255)
                                                                                                                                   ccAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATg                >  1:653811/1‑150 (MQ=255)
                                                                                                                                          ccTGTGGGGCCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGcc         <  1:158191/150‑1 (MQ=255)
                                                                                                                                          ccTGTGGGGCCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGcc         <  1:1364090/150‑1 (MQ=255)
                                                                                                                                            tgtgGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGa       <  1:822523/150‑1 (MQ=255)
                                                                                                                                            tgtgGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGa       >  2:1140480/1‑150 (MQ=255)
                                                                                                                                            tgtgGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGa       >  1:1052478/1‑150 (MQ=255)
                                                                                                                                              tgGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGACCCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAtt     >  2:1102555/1‑150 (MQ=255)
                                                                                                                                              tgGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAtt     >  2:1102516/1‑150 (MQ=255)
                                                                                                                                                ggggCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAttgt   >  2:313961/1‑150 (MQ=255)
                                                                                                                                                ggggCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAttgt   >  2:10669/1‑150 (MQ=255)
                                                                                                                                                ggggCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAttgt   >  2:659485/1‑150 (MQ=255)
                                                                                                                                                ggggCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAttgt   >  2:659569/1‑150 (MQ=255)
                                                                                                                                                ggggCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGAACCTACAGCGAACCGCACTGGGAACGGTCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAttgt   >  2:502969/1‑150 (MQ=255)
                                                                                                                                                 gggCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAttgtt  <  1:995776/150‑1 (MQ=255)
                                                                                                                                                 gggCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAttgtt  <  1:1602138/150‑1 (MQ=255)
                                                                                                                                                 gggCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAttgtt  <  1:1054003/150‑1 (MQ=255)
                                                                                                                                                 gggCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAttgtt  <  1:1602244/150‑1 (MQ=255)
                                                                                                                                                 gggCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGAttgtt  >  2:1235853/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACGCACCTACAGCGAACCGCACTGGGAACGGGCGCAGGGCGCGGTGATGGCAACGGAAAAAGTCACTGTTTATGGTTTGCCGATTGTT  >  NC_000913/1484956‑1485250

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