Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,532,275 T→A 6.8% intergenic (+51/‑134) cysZ → / → cysK sulfate transporter, sulfite inhibited/cysteine synthase A, O‑acetylserine sulfhydrolase A subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,532,2750TA 6.8% 259.5 / 8.6 74intergenic (+51/‑134)cysZ/cysKsulfate transporter, sulfite inhibited/cysteine synthase A, O‑acetylserine sulfhydrolase A subunit
Reads supporting (aligned to +/- strand):  ref base T (30/39);  new base A (3/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

TCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTT  >  NC_000913/2532127‑2532423
                                                                                                                                                    |                                                                                                                                                    
tCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTAt                                                                                                                                                      <  1:762629/149‑1 (MQ=255)
 cccGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATg                                                                                                                                                     <  2:466766/149‑1 (MQ=255)
  ccGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgc                                                                                                                                                    <  4:613426/149‑1 (MQ=255)
   cGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTCAAACAGGGGTCGCTTATgcc                                                                                                                                                   <  2:502613/149‑1 (MQ=255)
       gctTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGccc                                                                                                                                               <  2:448596/149‑1 (MQ=255)
        ctTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGCGGCATAAGCCAcccc                                                                                                                                              <  1:28174/149‑1 (MQ=255)
         ttAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCt                                                                                                                                             <  3:395736/149‑1 (MQ=255)
         ttAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCt                                                                                                                                             <  1:642703/149‑1 (MQ=255)
          tAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACAGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTGTGCCGCCCCtt                                                                                                                                            <  3:526453/149‑1 (MQ=255)
            aTCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTAt                                                                                                                                          >  2:57510/1‑149 (MQ=255)
            aTCTGTTCATCATGCCCGTTGCCGTTTGTGGAGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCCGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTAt                                                                                                                                          >  4:490491/1‑149 (MQ=255)
              cTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgccgc                                                                                                                                                 >  2:103122/1‑140 (MQ=255)
              cTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATgccgc                                                                                                                                                 <  1:103122/140‑1 (MQ=255)
              cTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTc                                                                                                                                        <  1:393421/149‑1 (MQ=255)
               tGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTcc                                                                                                                                       >  4:556089/1‑149 (MQ=255)
                gTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCa                                                                                                                                      >  3:548848/1‑149 (MQ=255)
                  tcatcaTGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATc                                                                                                                                    <  2:113127/149‑1 (MQ=255)
                      catGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGc                                                                                                                                >  2:369753/1‑149 (MQ=255)
                       atGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCa                                                                                                                               <  1:755135/149‑1 (MQ=255)
                           ccGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTc                                                                                                                           >  4:737335/1‑149 (MQ=255)
                              ttGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCatt                                                                                                                        <  1:378867/149‑1 (MQ=255)
                               tGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCatta                                                                                                                       <  3:233927/149‑1 (MQ=255)
                                   gTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTc                                                                                                                   <  3:618435/149‑1 (MQ=255)
                                       gtgGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCtt                                                                                                               <  3:363278/149‑1 (MQ=255)
                                        tgGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTc                                                                                                              <  2:126769/149‑1 (MQ=255)
                                                cGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGtatata                                                                                                      <  2:566807/149‑1 (MQ=255)
                                                cGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGtatata                                                                                                      <  4:548848/149‑1 (MQ=255)
                                                   cGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGat                                                                                                   <  4:219432/149‑1 (MQ=255)
                                                      tgtgGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATg                                                                                                >  2:584986/1‑149 (MQ=255)
                                                        tgGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCt                                                                                              <  1:735201/149‑1 (MQ=255)
                                                               aTTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCt                                                                                       <  1:427980/149‑1 (MQ=255)
                                                               aTTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCt                                                                                       >  1:85561/1‑149 (MQ=255)
                                                                      tCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTTTTCCATCTGGCATGTCATTATTTCCCTTCTGTATATAGATATTCTAAATCCTTACTTcc                                                                                >  2:110300/1‑149 (MQ=255)
                                                                               acacGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATtct                                                                       <  3:737335/149‑1 (MQ=255)
                                                                                cacGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATtctc                                                                      >  1:623772/1‑149 (MQ=255)
                                                                                       tgtgGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATTGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCgg                                                               <  2:85561/149‑1 (MQ=255)
                                                                                               tAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTa                                                       >  1:481998/1‑149 (MQ=255)
                                                                                                aaCAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTAc                                                      >  1:94372/1‑149 (MQ=255)
                                                                                                     tCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCtgtt                                                 <  2:240120/149‑1 (MQ=255)
                                                                                                     tCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCtgtt                                                 <  2:778648/149‑1 (MQ=255)
                                                                                                     tCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCtgtt                                                 <  3:243542/149‑1 (MQ=255)
                                                                                                       tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGa                                                                   >  4:321891/1‑129 (MQ=255)
                                                                                                       tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGa                                                                   <  3:321891/129‑1 (MQ=255)
                                                                                                       tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTAc                                                      <  2:238352/142‑1 (MQ=255)
                                                                                                       tACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTAc                                                      >  1:238352/1‑142 (MQ=255)
                                                                                                          cGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCt                                                                                                             <  3:261507/84‑1 (MQ=255)
                                                                                                          cGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCt                                                                                                             >  4:261507/1‑84 (MQ=255)
                                                                                                          cGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGt                                                             <  3:110595/132‑1 (MQ=255)
                                                                                                          cGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGt                                                             >  4:110595/1‑132 (MQ=255)
                                                                                                          cGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATc                                            <  3:490491/149‑1 (MQ=255)
                                                                                                          cGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATc                                            <  1:57510/149‑1 (MQ=255)
                                                                                                            gTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTctct                                                                     >  3:27300/1‑122 (MQ=255)
                                                                                                            gTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTctct                                                                     <  4:27300/122‑1 (MQ=255)
                                                                                                            gTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCtgtt                                                 >  1:504394/1‑142 (MQ=255)
                                                                                                            gTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCtgtt                                                 <  2:504394/142‑1 (MQ=255)
                                                                                                             ttATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCa                                         >  1:564462/1‑149 (MQ=255)
                                                                                                                    gTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCAt                                  >  4:737125/1‑149 (MQ=255)
                                                                                                                      aaaCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATAc                                >  2:347951/1‑149 (MQ=255)
                                                                                                                      aaaCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATAc                                >  2:730477/1‑149 (MQ=255)
                                                                                                                       aaCCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTGGTATCCCAATTTCATACa                               >  1:625265/1‑149 (MQ=255)
                                                                                                                            tttGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTAt                                                           >  3:445713/1‑116 (MQ=255)
                                                                                                                            tttGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTAt                                                           <  4:445713/116‑1 (MQ=255)
                                                                                                                            tttGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTaa                          <  2:623772/149‑1 (MQ=255)
                                                                                                                                  tgAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAg                    >  4:430659/1‑149 (MQ=255)
                                                                                                                                   gAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAgg                   <  3:713524/149‑1 (MQ=255)
                                                                                                                                     aaCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGcc                 >  4:328813/1‑149 (MQ=255)
                                                                                                                                     aaCAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGcc                 >  4:415584/1‑149 (MQ=255)
                                                                                                                                       cAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCAt               <  2:564462/149‑1 (MQ=255)
                                                                                                                                        aGGGGCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATg              >  1:563317/1‑149 (MQ=255)
                                                                                                                                         ggggCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATc                                            >  3:354051/1‑118 (MQ=255)
                                                                                                                                         ggggCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATc                                            <  4:354051/118‑1 (MQ=255)
                                                                                                                                          gggCGGCATAAGCCACCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCa                                   >  3:360681/1‑126 (MQ=255)
                                                                                                                                             tGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCACACAGTTAAGGACAGGCCATGAGTaa         <  2:416696/149‑1 (MQ=255)
                                                                                                                                                 ttATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAtt     >  1:778500/1‑149 (MQ=255)
                                                                                                                                                 ttATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAtt     >  2:163901/1‑149 (MQ=255)
                                                                                                                                                    tGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGAttttt  <  1:163901/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCCCGCTGCTTAATCTGTTCATCATGCCCGTTGCCGTTTGTGGCGCGACGGCGATGTGGGTCGATTGCTATCGCGATAAACACGCGATGTGGCGGTAACAATCTACCGGTTATTTTGTAAACCGTTTGTGTGAAACAGGGGTGGCTTATGCCGCCCCTTATTCCATCTTGCATGTCATTATTTCCCTTCTGTATATAGATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTT  >  NC_000913/2532127‑2532423

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