Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 163,417 A→T 5.1% K438I (AAA→ATA) ‡ hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913163,4170AT 5.1% 294.1 / 4.5 97K438I (AAA→ATA) ‡hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (53/39);  new base T (3/2);  total (56/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

AAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGCTTA  >  NC_000913/163272‑163566
                                                                                                                                                 |                                                                                                                                                     
aaaCGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAAtt                                                                                                                                                   <  2:160657/150‑1 (MQ=255)
  aCGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTgc                                                                                                                                                 <  1:232045/150‑1 (MQ=255)
    ggcTGATACAAATGCTGGCGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTgccg                                                                                                                                               <  2:1231803/148‑1 (MQ=255)
    gTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTgccg                                                                                                                                               >  1:1319142/1‑150 (MQ=255)
     tCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTgccgc                                                                                                                                              >  1:1605867/1‑150 (MQ=255)
        gTTACAACTGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCCGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAATTTTTgccgcggt                                                                                                                                           <  1:1089218/150‑4 (MQ=255)
        gTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAATTTTTgccgcggt                                                                                                                                           <  1:1277816/150‑4 (MQ=255)
        gTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCAt                                                                                                                                           >  1:78990/1‑150 (MQ=255)
            cAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATtctc                                                                                                                                       >  2:204337/1‑150 (MQ=255)
            cAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATtctc                                                                                                                                       >  2:361526/1‑150 (MQ=255)
             aaaTGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCg                                                                                                                                      <  2:814896/150‑1 (MQ=255)
              aaTGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGCCGACTAACCTGCTACCGCGGCAAAAATTGCCGCCAttaacga                                                                                                                                     >  1:1121922/1‑145 (MQ=255)
                 gCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAaga                                                                                                                                  >  2:1671860/1‑150 (MQ=255)
                  cTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAagag                                                                                                                                 <  1:379791/150‑1 (MQ=255)
                     ggggCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAgccg                                                                                                                              <  2:1668174/150‑1 (MQ=255)
                      gggCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAgccgc                                                                                                                             >  1:1123191/1‑150 (MQ=255)
                         cACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACg                                                                                                                         <  1:1159468/150‑1 (MQ=255)
                             ggAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGAt                                                                                                                      >  1:164835/1‑150 (MQ=255)
                             ggAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGAt                                                                                                                      >  1:1683625/1‑150 (MQ=255)
                              gAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATg                                                                                                                     <  2:1365649/150‑1 (MQ=255)
                                gggTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACTATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATggg                                                                                                                   <  2:708806/150‑1 (MQ=255)
                                 ggTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCCGCGCTGGGTAACGAGCCGCGGTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCGGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGc                                                                                                                  <  2:275230/150‑1 (MQ=255)
                                 ggTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGc                                                                                                                  <  2:1576672/150‑1 (MQ=255)
                                  gTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCa                                                                                                                 >  1:1397902/1‑150 (MQ=255)
                                  gTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCa                                                                                                                 >  1:1168218/1‑150 (MQ=255)
                                  gTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAACAGCCGCCACGGATGGGCa                                                                                                                 >  2:855745/1‑150 (MQ=255)
                                  gTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTATCGAAGAGCCGCCACGGATGGGCa                                                                                                                 >  1:910685/1‑150 (MQ=255)
                                  gTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCAGTCTCGAAGAGCCGCCACGGATGGGCa                                                                                                                 >  1:647476/1‑150 (MQ=255)
                                    gatcGGCGTAGTGCTCATGGCCACAAAGTGGCAGCGCTGGGCAACGCTGCGCGTGTCTCGGCAATGCTGGGTATCGTGAAGAACGACGACGAAGCGGCTACCGCGGCTAAAATTGTCGCCATTCTCGAAGAGCCGCCACGGATGGGCAAt                                                                                                               <  1:282130/147‑1 (MQ=255)
                                      aCGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAg                                                                                                             <  1:1411510/150‑1 (MQ=255)
                                           ttAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGCACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGAcc                                                                                                        >  1:1247009/1‑150 (MQ=255)
                                           ttAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGAcc                                                                                                        <  1:832132/150‑1 (MQ=255)
                                           ttAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGAcc                                                                                                        <  1:1422443/150‑1 (MQ=255)
                                            tAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCt                                                                                                       >  1:712202/1‑150 (MQ=255)
                                             aGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTg                                                                                                      >  1:1370941/1‑150 (MQ=255)
                                              gTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTgg                                                                                                     >  1:1469333/1‑150 (MQ=255)
                                              gTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTgg                                                                                                     >  2:1319421/1‑150 (MQ=255)
                                               ggcgcAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTggg                                                                                                    <  1:1434706/149‑1 (MQ=255)
                                                      gggCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggcgtggcg                                                                                             <  2:78990/150‑1 (MQ=255)
                                                      gggCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggcgtggcg                                                                                             <  1:1324760/150‑1 (MQ=255)
                                                      gggCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggcgtggcg                                                                                             <  1:1232862/150‑1 (MQ=255)
                                                      gggCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGGTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggcgtggcg                                                                                             <  1:321702/150‑1 (MQ=255)
                                                       ggCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggcgtggcgt                                                                                            >  1:1069772/1‑150 (MQ=255)
                                                        gCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCTCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTTCCTCAATTCTCGAAGAGACGCCACGGATGGGCAATAGTGACCTGGGAGTGGCGtt                                                                                           >  2:153545/1‑150 (MQ=255)
                                                        gCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGtt                                                                                           >  1:663918/1‑150 (MQ=255)
                                                         cAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGCAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGCTGGGCAATAGTGACCTGGGCGTGGCGgtt                                                                                          >  1:86120/1‑150 (MQ=255)
                                                         cAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGCGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGttt                                                                                          >  1:1377416/1‑150 (MQ=255)
                                                          aaaaaaTGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGtttt                                                                                         >  1:359640/1‑150 (MQ=255)
                                                           aaaaaTGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTc                                                                                        <  1:469925/150‑1 (MQ=255)
                                                            aaaaTGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTcg                                                                                       >  2:968119/1‑150 (MQ=255)
                                                            aaaaGGGCAGCGCGGGGGAACGAGCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTcg                                                                                       <  1:1642388/150‑1 (MQ=255)
                                                             aaaTGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTcgc                                                                                      >  1:1489067/1‑150 (MQ=255)
                                                              actGGCCGCGCTGGGTCACGATCCGCGGTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGATAGCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTcgcg                                                                                     <  1:1482764/148‑1 (MQ=255)
                                                              aaTGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTcgcg                                                                                     <  2:372300/150‑1 (MQ=255)
                                                               aTGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCgcgc                                                                                    <  2:209083/150‑1 (MQ=255)
                                                                  gCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAAt                                                                                 >  1:347129/1‑150 (MQ=255)
                                                                      cgcTGGGTAACGATCCGCGTGTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAAc                                                                             <  2:116271/150‑1 (MQ=255)
                                                                       gcTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAAcc                                                                            <  2:164835/150‑1 (MQ=255)
                                                                         tGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAg                                                                          >  2:451180/1‑150 (MQ=255)
                                                                         ggggTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAg                                                                          <  1:235992/149‑1 (MQ=255)
                                                                          gggTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACAGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGc                                                                         >  2:855805/1‑150 (MQ=255)
                                                                           ggTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGcc                                                                        >  2:844164/1‑150 (MQ=255)
                                                                             tAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTg                                                                      >  1:1537569/1‑150 (MQ=255)
                                                                              aaCGATCCGCGTTTAGCGGCAATGCTGTGTAGCGCGAAGAACTACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTgg                                                                     <  1:766875/150‑1 (MQ=255)
                                                                               aCGATCAGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGCCGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGgc                                                                    <  2:277663/150‑1 (MQ=255)
                                                                                     cgcgTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAAc                                                              >  2:773395/1‑150 (MQ=255)
                                                                                       cgTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGt                                                            >  2:939440/1‑150 (MQ=255)
                                                                                         tttAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAg                                                          >  1:844188/1‑150 (MQ=255)
                                                                                           tAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTc                                                        <  1:1121483/150‑1 (MQ=255)
                                                                                           tAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTc                                                        <  2:1714855/150‑1 (MQ=255)
                                                                                           tAGCGGCAATGCTGGTTAGAGCGAAGAACGACGCCGAAGCTGCTAACGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCACCGTAGTc                                                        >  2:36914/1‑150 (MQ=255)
                                                                                            aGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCa                                                       <  2:249193/150‑1 (MQ=255)
                                                                                              cGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTTGCGTTTACGCGCAATCAACCAGCATGGCAGCAACGTAGTCAGc                                                     >  2:1064736/1‑150 (MQ=255)
                                                                                              cGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCACGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGc                                                     >  2:799544/1‑150 (MQ=255)
                                                                                                   aTGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTg                                                <  1:204337/150‑1 (MQ=255)
                                                                                                      cTGGTTAGCGCGAAGAACGACGACTAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTa                                             <  2:1247009/150‑1 (MQ=255)
                                                                                                       tGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTaa                                            >  1:1388863/1‑150 (MQ=255)
                                                                                                        ggTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGAGTGGCGTTTTCGCGCAATAAACCAGCCTGGCAGCAACGTAGTCAGCAACTCTTaaa                                           >  2:1589455/1‑150 (MQ=255)
                                                                                                           tAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACg                                        <  1:451180/150‑1 (MQ=255)
                                                                                                            aGCGCGAAGAACGATGACGAAGCTGCTACCGCGACAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGc                                       <  1:1156571/150‑1 (MQ=255)
                                                                                                             gcgcGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGCCCTGGGCGTGGCGTTTTCGCGCCAGCCACCAGCCTGGCAGCCACGTAGTCAGCACCTGTTAAAACGCt                                      >  1:1602454/1‑150 (MQ=255)
                                                                                                                cgAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTaa                                   >  1:1373399/1‑150 (MQ=255)
                                                                                                                    gaaCGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGACGCGCCGCCACGGCTGGGCAGTAGTGCCCTGGGCGTCGCGGGGTCGCGCAATCAAGCAGGCTGGCAGCAACGTAGTCCGCCACTGGTAAACCGCTGACacgt                               >  1:329286/1‑150 (MQ=255)
                                                                                                                         acgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTg                          >  1:1662648/1‑150 (MQ=255)
                                                                                                                         acgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTg                          >  1:1607505/1‑150 (MQ=255)
                                                                                                                           gacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGc                        >  1:372253/1‑150 (MQ=255)
                                                                                                                                aGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTCCGTGGCGGTGa                   >  1:246547/1‑150 (MQ=255)
                                                                                                                                aGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGGTTTCGCGCAATCAACCAGCCGGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTgg                   >  2:764239/1‑149 (MQ=255)
                                                                                                                                     gaatggcgGCAATTTTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAACCGCTTAAACGTACGTGGCGGTGAGGCAg              >  1:1333706/6‑150 (MQ=255)
                                                                                                                                     ctACCGCGGCAAAATTTGCCGCCATTCTCGCAGCGCCGCCACGGACGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCCACCAGCCGGGCAGCAACGTAGTCAGCAACTGTTAACACGCTTAAACTTACGTGGCGGTGAGGCAg              <  1:1489562/150‑1 (MQ=255)
                                                                                                                                       atggcgGCAATTTTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGAc            >  1:165639/4‑150 (MQ=255)
                                                                                                                                       atggcgGCAATTTTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGAAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACCTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGAc            >  1:1025904/4‑150 (MQ=255)
                                                                                                                                       aCCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCGCCCAGCCTGGCCGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGAc            >  1:135040/1‑150 (MQ=255)
                                                                                                                                           cgGCAAAAATTGCCGCCATTCTCGAAGAGCCGCGACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCGGGCAGCAACCTAGTCAGCAACTGTTAAAACGCTTAAACGAACGTGGCGGTGAGGCAGACAGtt        >  1:735144/1‑150 (MQ=255)
                                                                                                                                           cgGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGCCCGGGGCGTGGCGTGTGCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACGGTTAAACCGCTTAAACGGACGGGGCGGTGAGGCAGACGGtt        >  2:1082097/1‑150 (MQ=255)
                                                                                                                                            ggCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTc       <  2:1373399/150‑1 (MQ=255)
                                                                                                                                             gCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCg      >  2:484060/1‑150 (MQ=255)
                                                                                                                                                 aaaTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGCTTa  >  2:1183249/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGCTTA  >  NC_000913/163272‑163566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: