Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,270,146 T→C 6.1% intergenic (+19/‑92) aphA → / → yjbQ acid phosphatase/phosphotransferase/UPF0047 protein YjbQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,270,1460TC 6.1% 310.9 / 7.2 110intergenic (+19/‑92)aphA/yjbQacid phosphatase/phosphotransferase/UPF0047 protein YjbQ
Reads supporting (aligned to +/- strand):  ref base T (53/49);  new base C (3/4);  total (57/53)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-02

Only 100 of 112 total aligned reads displayed.

TACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTC  >  NC_000913/4270004‑4270286
                                                                                                                                              |                                                                                                                                             
tACCGCCGCACTCGATGTCGGCGCTCGTGGTATCCGCCTTCTGCGCGCCTCCAACTCTCCCTACAAACCCTAGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGGCAATTCAGAATACTGACAGAGCGGGAGAGCGTGAt                                                                                                                                               <  1:668836/143‑1 (MQ=255)
tACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGAt                                                                                                                                               <  1:741595/143‑1 (MQ=255)
 acggccgcACGCGATGTCGGAGCCCGTGGTAGCCGCATGCTGCGCGCCTCCAGCTCTACCTACAAACCCGTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATg                                                                                                                                              <  1:364555/140‑1 (MQ=255)
  ccgccgAACGCGATGTCGGCGGTCGGGGTGTCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGACGAGGAGACCGTCAATTCAGAATACTGACAGAGCGGGCGCGCGTGATGc                                                                                                                                             >  2:316939/1‑143 (MQ=255)
        caggcgATGTCGGCGCTCGTGGTATCCGCCTTCTGCGCGCCGCCAACGCTACCTACAAACCCTTGCCACACGCGGGTGCGTTTGGTGAAGAGGGGATCGTCGATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTccc                                                                                                                                       <  2:333387/140‑1 (MQ=255)
          cgcgATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGa                                                                                                                                     <  1:811548/143‑1 (MQ=255)
           gcgATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgct                                                                                                                                    <  1:866928/143‑3 (MQ=255)
           gccaTGCCGGCGCTCGCGGTATCCGCATTCTGCGCGCCTCCAACTCTCCCTACAAACACTTGCCACAAGCGGGTGCGTTTTGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGaa                                                                                                                                    <  1:166336/140‑1 (MQ=255)
            cgATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgctc                                                                                                                                   <  1:503349/143‑4 (MQ=255)
            cgATGTCGGCGCTCGCGGTATGCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAAt                                                                                                                                   <  1:122573/143‑1 (MQ=255)
             gATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATg                                                                                                                                  <  1:788407/143‑1 (MQ=255)
                gTCAGCGCTCGTGGGATCCGCATTCTGCGCGCCTCCAACTCTACCTACAGACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTg                                                                                                                               <  2:341844/143‑1 (MQ=255)
                 tCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGt                                                                                                                              >  2:822315/1‑143 (MQ=255)
                 tCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGGGATGCTCTCCCGAATGCTGt                                                                                                                              >  1:827649/1‑143 (MQ=255)
                    gcgcTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttt                                                                                                                           <  1:829984/143‑1 (MQ=255)
                       cTCGTGGTAGCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttttttt                                                                                                                        >  2:850924/1‑143 (MQ=255)
                        tCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACGAGCGGGGGCGTTTGGTGCAGAGGTGATCGTCAATTCAGACTGCTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttttttc                                                                                                                       >  1:758173/1‑142 (MQ=255)
                        tCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGACTACTGACAGAGCGGGGGAGCGTGGTGCGCTCCCGACGGCTGtttttgtc                                                                                                                       >  2:234009/1‑140 (MQ=255)
                          gTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAAt                                                                                                                     <  2:862643/143‑1 (MQ=255)
                            ggTATCCGGATGCGGCGCGCCTCCAACTCTACCGACAAACCCTGGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATca                                                                                                                   <  2:61341/143‑1 (MQ=255)
                            ggTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATca                                                                                                                   >  1:781180/1‑143 (MQ=255)
                               aTCCGGATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCacac                                                                                                                <  2:614546/143‑1 (MQ=255)
                                tCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACAcc                                                                                                               <  2:745300/143‑1 (MQ=255)
                                  cGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGATGCGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGACAGCTGTTTTTGTAATccagcctc                                                                                                             >  2:726332/1‑136 (MQ=255)
                                  cGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACAccgt                                                                                                             >  1:570003/1‑141 (MQ=255)
                                  cGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCtt                                                                                                             >  2:806636/1‑143 (MQ=255)
                                   gCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTACTCACACCttt                                                                                                            >  1:38806/1‑143 (MQ=255)
                                   gCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCttt                                                                                                            <  1:522239/143‑1 (MQ=255)
                                      ttCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTtctc                                                                                                         >  2:606013/1‑143 (MQ=255)
                                         cggggggcTCCAGCTCTGCCTCCAATCCCTGGCCACGAGCGGGTGCGTGTGGTGAAGAGGTGATGGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCtt                                                                                                      <  2:708759/136‑1 (MQ=255)
                                           cgcgcCTATAACTCTACCTAGAAACACTCGCCACAAGCGAGAGCGTTGGGTGAAGAGGTGATCGCCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTc                                                                                                    <  1:63732/143‑1 (MQ=255)
                                              gcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCt                                                                                                 <  1:228789/143‑1 (MQ=255)
                                              gcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGCTTGGTGCACAGGTGATCGTCAATTCAGAATACTGACCGAGCGGGAGAGCGAGCAGAGCCCTCGACTGCTGTTCTCTTACTCACACCTTCATCCTTCCGCt                                                                                                 >  2:861891/1‑143 (MQ=255)
                                                 tCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGACTACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTACTCACACCTTTATCCTTTCGCTGTc                                                                                              >  1:302593/1‑143 (MQ=255)
                                                 tCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTc                                                                                              >  2:217816/1‑143 (MQ=255)
                                                 tCCAACTCTAACTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCATACACGAATGCTGTTTTGTTAATCGCAGCTTACTCATTTCGCtcgc                                                                                              >  2:43496/1‑140 (MQ=255)
                                                  ccAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTGTCGCTGTCt                                                                                             >  1:27058/1‑143 (MQ=255)
                                                    gacgctacctacAAACCCTTGCCACAAGCGGGTGCGTTGGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtg                                                                                           <  1:874451/139‑1 (MQ=255)
                                                    aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTGTTTTAATCACACCTTTATCCTTTCGcgggcgtg                                                                                           >  1:148274/1‑136 (MQ=255)
                                                     aCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgc                                                                                          >  1:778321/1‑143 (MQ=255)
                                                        ctacctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCCGTTTTTTTAATCACACCTTTATCCTTCCGCTCTCTtgctgc                                                                                       >  1:602152/1‑143 (MQ=255)
                                                        ctacctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAACGCTGTTTTTTTAATCACACCGTTATCCTTTCGCTGGCTtgctgc                                                                                       >  2:374611/1‑143 (MQ=255)
                                                         tacctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCCTTCGCTGTCTtgctgcc                                                                                      >  2:11959/1‑142 (MQ=255)
                                                            ctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGCTGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTCTCGCTGTTTTGCGGCAAAc                                                                                   >  2:768917/1‑143 (MQ=255)
                                                            ctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTACTCACACCTTTATCCTTTCGCTGTCTTGCTGCAAAc                                                                                   >  1:303035/1‑143 (MQ=255)
                                                             tgctggcccTTGGCAGACGCGGGTGCGTGTGGTGGAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACt                                                                                  <  2:401910/137‑1 (MQ=255)
                                                                gtgcccTTGCCACAAGCGGGTGCGTGTGGTGAAGCGGTGATCGTCAACTCAGAATACTTACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTACCCTTTCGCTGTCTTGCTGCAAACTGAt                                                                               <  1:308121/140‑1 (MQ=255)
                                                                aaaCCCATGCCACAAGCGGGTGCGTTTGGGGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAt                                                                               <  1:245033/143‑1 (MQ=255)
                                                                   cccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTTTCTCGCTGCAAACTGATTac                                                                            >  2:658616/1‑142 (MQ=255)
                                                                    ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGCTTAag                                                                           >  1:94777/1‑143 (MQ=255)
                                                                     cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAaga                                                                          <  2:559795/143‑1 (MQ=255)
                                                                       tGCCACACGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagaga                                                                        <  1:627050/143‑1 (MQ=255)
                                                                         ccACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGt                                                                      <  1:53305/143‑1 (MQ=255)
                                                                          cacaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGtt                                                                     >  1:612115/1‑143 (MQ=255)
                                                                           acaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGttt                                                                    <  1:286103/143‑1 (MQ=255)
                                                                               gCGGGTTCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACTCCTTTATCCTTTCGCTGTCTTCCTTCCAACTGCTTAAGAGAGTTTTATc                                                                >  2:590786/1‑143 (MQ=255)
                                                                               gCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATc                                                                <  2:794200/143‑1 (MQ=255)
                                                                                cGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAATAAGAGAGTTTTATCa                                                               >  1:743980/1‑143 (MQ=255)
                                                                                 gggTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCaa                                                              >  1:595410/1‑143 (MQ=255)
                                                                                 gggTGCGTTGGGAGAAGAGGTGATCGGCAATACGGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCaa                                                              <  1:814806/143‑1 (MQ=255)
                                                                                      cGTGTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagc                                                         <  2:566873/143‑1 (MQ=255)
                                                                                       gTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGCCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagca                                                        <  1:653829/143‑1 (MQ=255)
                                                                                        tttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCGTTTCGCTGTCTTGTTGCAAACTGACGAATATAGTGTTATCcagcagcag                                                       >  2:576727/1‑143 (MQ=255)
                                                                                        tttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTCAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagcag                                                       >  1:659276/1‑143 (MQ=255)
                                                                                         ttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagcagc                                                      >  1:709839/1‑143 (MQ=255)
                                                                                            gTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                                                                    >  1:855470/1‑110 (MQ=255)
                                                                                            gTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                                                                    <  2:855470/110‑1 (MQ=255)
                                                                                                   ggTGATCGTCAATTCAGAATACAAACAGGGGGGGAGAGCGTGCTCTGCGTCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTa                                            <  1:790808/143‑1 (MQ=255)
                                                                                                    gTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTAt                                           >  1:418508/1‑143 (MQ=255)
                                                                                                    gTGATCGTCAATTCAGAATACGGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCATGGAGCAGCACATGTGGTAt                                           <  1:519048/143‑1 (MQ=255)
                                                                                                     tGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCCTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGCGCTTTATCAAGGCGCAGCACATGTGGTATc                                          >  2:105141/1‑143 (MQ=255)
                                                                                                        tCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCGGCAGATGTGGGATCgaa                                       >  1:278382/1‑143 (MQ=255)
                                                                                                         cgtgacttCAGAATACTGACAGTGCGGGAGAGCGTGATGATCGCGCCGATGCTGTTGTGTTAATCACACCTTCATCCGTTCGCTGTCTTGCTTCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCaaaa                                      <  2:624701/137‑1 (MQ=255)
                                                                                                          gttaaTTTAGAATACTGACAGAGCGGGGGAGCGTGATCCACTACCGTATTCTTTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAg                                     <  2:816321/140‑1 (MQ=255)
                                                                                                          gTCAATTCAGGATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAg                                     <  2:620993/143‑1 (MQ=255)
                                                                                                          gTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACGTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGAGGTATCAAAAg                                     >  1:379683/1‑143 (MQ=255)
                                                                                                          gTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAg                                     >  2:532520/1‑143 (MQ=255)
                                                                                                           tCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGa                                    <  1:529879/143‑1 (MQ=255)
                                                                                                           tCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATACTTTCACTATCTTCCTGCAAACTGATTAAGCGAGATTTATCAAGGCGCAGCAACTGGGCTATAaaaaca                                    >  2:732369/1‑141 (MQ=255)
                                                                                                           tCAATGCAGAATACTGGCAGAGCGGGAGCGCGTGCTGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTCTCCGTGCGCTTTCGTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGa                                    <  2:351965/143‑1 (MQ=255)
                                                                                                            cAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGAc                                   >  1:179859/1‑143 (MQ=255)
                                                                                                             aaTTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACg                                  <  2:167154/143‑1 (MQ=255)
                                                                                                               ttCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCt                                >  1:110086/1‑143 (MQ=255)
                                                                                                               ttCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTCATCAAGGAGCCGCACATGTGGTATTAAAAGACGCt                                >  2:115366/1‑143 (MQ=255)
                                                                                                               ttCAGAATACTGACAGAGCGGGAGAGCGTGATGCGCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCt                                <  1:376058/143‑1 (MQ=255)
                                                                                                                 cAGCATACTGACAGAGCGGGAGAGCGTCATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTGTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCa                              <  2:327559/143‑1 (MQ=255)
                                                                                                                     aTACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAAGGCTGTGTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCt                         <  2:605247/143‑1 (MQ=255)
                                                                                                                       aCTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTa                       <  1:86038/143‑1 (MQ=255)
                                                                                                                       aCTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTATTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCACACTGATTTAGAGAGTTTTATCAAGGAGCCGACCATGTGGTAGCAACAGACACTCGCGCttt                       >  2:709873/1‑142 (MQ=255)
                                                                                                                             agagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCa                 >  1:568834/1‑143 (MQ=255)
                                                                                                                              attcGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTCTCAAGAAGCAGACCAGCGGGTCTCAAAACACGCCCACAACTAAGCGCCa                 >  2:408108/4‑143 (MQ=255)
                                                                                                                               ttcGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCaaa               >  1:150266/3‑143 (MQ=255)
                                                                                                                               agCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCTTTTCGTTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGTAGCACATGTAGTATCAAAATACGCTCACGCTTAGCGCCaat               >  1:293245/1‑142 (MQ=255)
                                                                                                                               agCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGCTTTATCAAGGAGCAGCACATGTGGTATCATAACAATGTCACGCTTAGCGCCaaa               >  2:811688/1‑143 (MQ=255)
                                                                                                                                  gggAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAACCTCATGAAGCGAGTTTTCTCGAGGCGCAGAACATTACGTAGCCAAAGACGCGCACACTGAGCCCCGAAtct            >  2:111781/1‑143 (MQ=255)
                                                                                                                                      gagCGTGATGCTCTCCTGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGATGCCAACTGATTAAGAGAGTTTTATCAAGGCGCAGCACATGTTGTATCAAACGTCGCTCACGCTTAGCGCCAAATCTCGTg        <  2:465688/143‑1 (MQ=255)
                                                                                                                                      gagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTg        <  1:217816/143‑1 (MQ=255)
                                                                                                                                      gagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTg        >  1:693507/1‑143 (MQ=255)
                                                                                                                                         cGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAGATCTCGGGGGt     >  1:309570/1‑143 (MQ=255)
                                                                                                                                            gATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGTAAACTGATTATTAGAGTTTTATCAAGTAGCAGCACATATGGTATCAAAAAACGTTTACGCTTCACGCCAAATCTCGTGGCtttt  >  2:780682/1‑142 (MQ=255)
                                                                                                                                              |                                                                                                                                             
TACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTC  >  NC_000913/4270004‑4270286

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