Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,453,330 A→G N254D (AAC→GAC)  ytfT → putative sugar ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,453,3300AG96.3% 375.2 / 3.7 109N254D (AAC→GAC) ytfTputative sugar ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (56/49);  total (57/52)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.59e-02

Only 100 of 112 total aligned reads displayed.

ATCGAAGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCAACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCTTTCGG  >  NC_000913/4453183‑4453477
                                                                                                                                                   |                                                                                                                                                   
  cGAAGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCgacaac                                                                                                                                                <  1:1379934/151‑1 (MQ=255)
   gAAGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACg                                                                                                                                               >  1:1400444/1‑151 (MQ=255)
    aaGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGc                                                                                                                                              >  2:113261/1‑151 (MQ=255)
     aGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCt                                                                                                                                             <  1:482773/151‑1 (MQ=255)
       ccGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTg                                                                                                                                            <  2:746161/150‑1 (MQ=255)
       ccGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTg                                                                                                                                            <  2:1525127/150‑1 (MQ=255)
         gTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTggg                                                                                                                                          >  1:1180888/1‑150 (MQ=255)
         gTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTggg                                                                                                                                          >  1:893933/1‑150 (MQ=255)
         gTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGt                                                                                                                                         >  2:1393670/1‑151 (MQ=255)
         gTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGt                                                                                                                                         >  1:743186/1‑151 (MQ=255)
          ttGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGtt                                                                                                                                        >  1:285399/1‑151 (MQ=255)
             gTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATg                                                                                                                                     <  2:1078016/151‑1 (MQ=255)
                 cAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctg                                                                                                                                 >  2:133316/1‑151 (MQ=255)
                  aaCATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctgg                                                                                                                                <  2:1701057/151‑1 (MQ=255)
                   aCATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctgga                                                                                                                               <  1:79883/151‑1 (MQ=255)
                     aTTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctggagc                                                                                                                             <  1:379208/151‑1 (MQ=255)
                     aTTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGAAAACGCTGGGTTATGgctggag                                                                                                                              >  1:1189655/1‑150 (MQ=255)
                         gggCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctggagctgga                                                                                                                         <  1:1608808/151‑1 (MQ=255)
                         gggCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGgctggagctgg                                                                                                                          >  1:1634596/1‑150 (MQ=255)
                                 aaaaaTGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATtc                                                                                                                 >  1:66796/1‑151 (MQ=255)
                                 aaaaaTGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATtc                                                                                                                 <  1:96049/151‑1 (MQ=255)
                                  aaaaTGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATtct                                                                                                                <  2:615858/151‑1 (MQ=255)
                                   aaaTGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATtctc                                                                                                               <  2:256151/151‑1 (MQ=255)
                                   aaaTGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATtctc                                                                                                               >  1:1069576/1‑151 (MQ=255)
                                     aTGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGAGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCAtttt                                                                                                                >  1:152516/1‑146 (MQ=255)
                                       gCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTcgcg                                                                                                            >  2:501693/1‑150 (MQ=255)
                                         cGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCgggg                                                                                                         >  1:1155392/1‑149 (MQ=255)
                                          ggggTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCggtgg                                                                                                        >  2:1075403/1‑151 (MQ=255)
                                           gggTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCggtggt                                                                                                       <  2:265326/151‑1 (MQ=255)
                                             gTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGa                                                                                                     <  1:415334/151‑1 (MQ=255)
                                               aaaCACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGAt                                                                                                    <  2:894013/150‑1 (MQ=255)
                                                 acacGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTgg                                                                                                 >  2:428849/1‑151 (MQ=255)
                                                 acacGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCAACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTgg                                                                                                 <  1:457949/151‑1 (MQ=255)
                                                    cgcgAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTggcg                                                                                               <  1:201130/150‑1 (MQ=255)
                                                         atcatcGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGGGCTGGACGCCATTCTCGCGGTGGTGATTGGcggcgg                                                                                           >  1:1125061/1‑149 (MQ=255)
                                                         atcatcGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGAt                                                                                         <  2:575082/151‑1 (MQ=255)
                                                            atcGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCg                                                                                       >  2:1128776/1‑150 (MQ=255)
                                                            atcGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCAACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGc                                                                                      >  2:1090486/1‑151 (MQ=255)
                                                             tcGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGc                                                                                      <  2:325202/150‑1 (MQ=255)
                                                                    gCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATgggc                                                                              <  1:1681863/151‑1 (MQ=255)
                                                                    gCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATggg                                                                               >  1:1702610/1‑150 (MQ=255)
                                                                    gCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATggg                                                                               <  2:387992/150‑1 (MQ=255)
                                                                      acACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGCCAACGCTGGGGTATGGCTGGAGCTTGAAGCCCATCTTGCGGTTGCGATTGGCGGCGGTTGGCGGGTgggcgg                                                                            >  1:1115193/2‑151 (MQ=255)
                                                                       cACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATgggcggg                                                                           <  2:211310/151‑1 (MQ=255)
                                                                        aCTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATgggcgggc                                                                          <  2:874465/151‑1 (MQ=255)
                                                                            aTGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCtt                                                                      >  2:1659723/1‑151 (MQ=255)
                                                                             tgtgTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCttt                                                                     >  1:1686873/1‑151 (MQ=255)
                                                                              gtgtTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCttt                                                                     <  1:1069865/150‑1 (MQ=255)
                                                                               tgtTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGGGATTGGCGGCGGGTCGCTGGTGGGCGGGCGCttt                                                                     >  1:106214/1‑149 (MQ=255)
                                                                                gtTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAAc                                                                  >  1:1618792/1‑151 (MQ=255)
                                                                                 ttGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGAAAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAAcc                                                                 >  2:1655377/1‑151 (MQ=255)
                                                                                  tGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAAcc                                                                 >  1:113224/1‑150 (MQ=255)
                                                                                  tGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCt                                                                >  2:86906/1‑151 (MQ=255)
                                                                                  tGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGGGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAAcc                                                                 >  1:1661842/1‑150 (MQ=255)
                                                                                   gAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCt                                                                <  1:71877/150‑1 (MQ=255)
                                                                                   gAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTg                                                               <  2:489796/151‑1 (MQ=255)
                                                                                      cGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCt                                                             >  1:1106979/1‑150 (MQ=255)
                                                                                        ggCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACt                                                          <  1:387467/151‑1 (MQ=255)
                                                                                            gtgtgCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCg                                                      <  2:1542942/151‑1 (MQ=255)
                                                                                            gtgtgCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCg                                                      <  2:732193/151‑1 (MQ=255)
                                                                                            gtgtgCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCg                                                      <  2:1306940/151‑1 (MQ=255)
                                                                                             tgtgCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCgg                                                     >  2:916090/1‑151 (MQ=255)
                                                                                              gtgCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCggt                                                    <  2:438576/151‑1 (MQ=255)
                                                                                                 cggcgATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCggggg                                                  >  1:365492/1‑150 (MQ=255)
                                                                                                  ggcgATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCggtggta                                                <  1:256931/151‑2 (MQ=255)
                                                                                                    cgATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTggg                                              <  2:256153/151‑1 (MQ=255)
                                                                                                    cgATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTggg                                              >  2:266163/1‑151 (MQ=255)
                                                                                                     gATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGGGCTGGACGCCATTCTCGCGGTGGTGATTGGCGACGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGgggggg                                             >  1:1145679/1‑151 (MQ=255)
                                                                                                        tGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGGGGATCGCTGATGGGCGGGCGCTTTAACGTGCTACTTTtgggggtggggggg                                          >  1:772987/1‑149 (MQ=255)
                                                                                                        tGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGgcg                                          >  2:1704172/1‑151 (MQ=255)
                                                                                                        tGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGgcg                                          >  2:1124730/1‑151 (MQ=255)
                                                                                                          cGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCt                                        <  2:1049269/151‑1 (MQ=255)
                                                                                                          cGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCt                                        >  2:1255513/1‑151 (MQ=255)
                                                                                                           gggCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTg                                       <  2:557212/151‑1 (MQ=255)
                                                                                                            ggCATTATCGTGGCGGCGGATATTCGCGGTGCCTATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTg                                       <  1:470185/150‑1 (MQ=255)
                                                                                                             gCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGat                                     <  2:952979/151‑1 (MQ=255)
                                                                                                               aTTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGatta                                   <  1:939551/151‑1 (MQ=255)
                                                                                                               aTTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGatt                                    >  2:1096014/1‑150 (MQ=255)
                                                                                                                  aTCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTc                                >  2:1376816/1‑151 (MQ=255)
                                                                                                                  aTCGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTc                                <  2:1235164/151‑1 (MQ=255)
                                                                                                                    cGTGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCa                               <  2:372719/150‑1 (MQ=255)
                                                                                                                     ttGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGGGCTGATTATTCAgg                             >  1:558601/2‑151 (MQ=255)
                                                                                                                      tGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAggg                            >  2:353876/1‑151 (MQ=255)
                                                                                                                      tGGCGGCGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAggg                            >  1:832105/1‑151 (MQ=255)
                                                                                                                       ggcggcGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAgggg                           <  1:113252/151‑1 (MQ=255)
                                                                                                                       ggcggcGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAgggg                           <  1:428814/151‑1 (MQ=255)
                                                                                                                       ggcggcGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAgggg                           <  2:1370500/151‑1 (MQ=255)
                                                                                                                       ggcggcGGATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAggg                            <  2:1400536/150‑1 (MQ=255)
                                                                                                                         cggcggATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGAt                         <  1:1322171/151‑1 (MQ=255)
                                                                                                                          ggcggATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATg                        <  1:1504896/151‑1 (MQ=255)
                                                                                                                            cggATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGaa                      <  1:934501/151‑1 (MQ=255)
                                                                                                                            cggATATTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGACCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGGGGGGGCGCTGATTATTCAGGGGATGa                       >  1:35412/1‑150 (MQ=255)
                                                                                                                            cggATATTCGCGGTGCCGATGCCAACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGaa                      <  1:1090401/151‑1 (MQ=255)
                                                                                                                                 atTCGCGGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCg                 >  2:1058112/1‑151 (MQ=255)
                                                                                                                                     gcgGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAAt             >  1:560195/1‑151 (MQ=255)
                                                                                                                                     gcgGTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAAt             >  1:1310388/1‑151 (MQ=255)
                                                                                                                                       ggTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAAttt           >  2:1444461/1‑151 (MQ=255)
                                                                                                                                        gTGCCGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAAtttt          <  2:521868/151‑1 (MQ=255)
                                                                                                                                           ccGATGCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCt       <  2:1634699/151‑1 (MQ=255)
                                                                                                                                                gCCGACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCTTTCgg  >  2:969216/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCGAAGCCGTTGGTATCAACATTCGGGCGGCAAAAAATGCCGGGGTAAACACGCGAATCATCGTCATGCTCACTTATGTGTTGAGCGGGCTGTGTGCGGCGATTGCGGGCATTATCGTGGCGGCGGATATTCGCGGTGCCGATGCCAACAACGCTGGGTTATGGCTGGAGCTGGACGCCATTCTCGCGGTGGTGATTGGCGGCGGATCGCTGATGGGCGGGCGCTTTAACCTGCTACTTTCGGTGGTGGGGGCGCTGATTATTCAGGGGATGAACACCGGAATTTTGCTTTCGG  >  NC_000913/4453183‑4453477

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