Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,576 T→C 5.4% intergenic (+21/+220) ybiT → / ← ybiU ABC‑F family putative regulatory ATPase/DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913857,5760TC 5.4% 281.8 / 4.8 92intergenic (+21/+220)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base T (39/48);  new base C (3/2);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

CACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAA  >  NC_000913/857437‑857722
                                                                                                                                           |                                                                                                                                                  
cacaCGAGCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgg                                                                                                                                            <  6:241023/148‑1 (MQ=255)
cacaCGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgg                                                                                                                                            <  3:324015/148‑1 (MQ=255)
 acacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgggg                                                                                                                                          >  7:148597/1‑149 (MQ=255)
 acacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCggg                                                                                                                                           <  4:161753/148‑1 (MQ=255)
  cacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                                                                         <  2:193325/149‑1 (MQ=255)
  cacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                                                                         <  1:152256/149‑1 (MQ=255)
   acGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTg                                                                                                                                        <  6:136640/149‑1 (MQ=255)
    cGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGa                                                                                                                                       <  8:148597/149‑1 (MQ=255)
     gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAc                                                                                                                                      >  8:89566/1‑149 (MQ=255)
     gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAc                                                                                                                                      >  1:92423/1‑149 (MQ=255)
     gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAc                                                                                                                                      <  7:337842/149‑1 (MQ=255)
     gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAc                                                                                                                                      >  7:47155/1‑149 (MQ=255)
     gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAc                                                                                                                                      >  7:279434/1‑149 (MQ=255)
       ccGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcgg                                                                                                                                    >  6:26053/1‑149 (MQ=255)
       ccGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAcgg                                                                                                                                    <  1:290589/149‑1 (MQ=255)
          tGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcggcgg                                                                                                                                 <  2:109881/149‑1 (MQ=255)
               tCGTAATCTTCCTGGCGACCCGCATTCTTGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgat                                                                                                                            <  2:223354/149‑1 (MQ=255)
                 gTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACggcggc                                                                                                                                <  6:8711/143‑1 (MQ=255)
                 gTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACggcggc                                                                                                                                >  5:8711/1‑143 (MQ=255)
                  tAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGAc                                                                                                                         <  2:239327/149‑1 (MQ=255)
                   aaGCTCCCTGGCGCCACGCATTCTGGCAATCACCCCGGAACGCGTGATAGACTTTAGCGGTAATTAAGAAGATTACCTGAGTAGTAAAGGGATCGAGGAAGTTATTATACGTCACCACGATACTGACGGGGTGCCGGCGGAtgatga                                                                                                                          >  5:227911/1‑147 (MQ=255)
                      cTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAAt                                                                                                                     <  6:196590/149‑1 (MQ=255)
                      cTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAAt                                                                                                                     <  4:224562/149‑1 (MQ=255)
                       tCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATg                                                                                                                    >  2:136667/1‑149 (MQ=255)
                         ccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGc                                                                                                                   >  2:121004/1‑148 (MQ=255)
                          cTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCa                                                                                                                  >  2:74802/1‑148 (MQ=255)
                           tGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                                                                                 >  1:153119/1‑148 (MQ=255)
                              cGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATa                                                                                                             <  1:86621/149‑1 (MQ=255)
                                 cccGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGcc                                                                                                          <  1:136667/149‑1 (MQ=255)
                                  ccGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                                                                         >  4:61890/1‑149 (MQ=255)
                                  ccGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGAGTAGTAAAGGTATCGAGTAAGTGAGTATCCGTCACCCCGATACTGTCGGGGTGAATGCGGCTGATGACAAATGCAAAATAGCCt                                                                                                         >  2:141675/1‑149 (MQ=255)
                                       ttCTGGAAATCACCCCGGACCGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGAGGCTGATGACAAATGCAAAATAGCCTGATgc                                                                                                    >  5:74482/1‑149 (MQ=255)
                                          tGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGAtaagcg                                                                                                 >  7:84600/1‑144 (MQ=255)
                                          tGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgcgc                                                                                                  >  7:188811/1‑148 (MQ=255)
                                                  cccccGTACCGCTTGATCGACTTTACCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGAGCGAGTAAGTGATTATCCGTAACCCCGATACTGTCGGTGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTa                                                                                          <  5:209190/147‑1 (MQ=255)
                                                  aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                                                                         <  1:12246/101‑1 (MQ=255)
                                                  aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                                                                         >  2:12246/1‑101 (MQ=255)
                                                  aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                         >  5:337444/1‑149 (MQ=255)
                                                   ccccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTCTCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATc                                                                                        <  6:87444/149‑1 (MQ=255)
                                                    cccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCa                                                                                       <  3:149651/149‑1 (MQ=255)
                                                     ccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAg                                                                                      <  1:74802/149‑1 (MQ=255)
                                                       ggAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                                                                    <  1:121004/149‑1 (MQ=255)
                                                         aaCGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCt                                                                                  <  6:337444/149‑1 (MQ=255)
                                                           cgcgTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                <  6:115794/149‑1 (MQ=255)
                                                              gTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACg                                                                               >  3:79402/1‑147 (MQ=255)
                                                                gATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAg                                                                            >  8:19481/1‑148 (MQ=255)
                                                                  tCGACTTTAGCGGTAATTAAGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCACGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                         >  1:259841/1‑149 (MQ=255)
                                                                    gACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTACGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCt                                                                       <  6:243282/149‑1 (MQ=255)
                                                                     aCTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTg                                                                      >  7:37354/1‑149 (MQ=255)
                                                                      cTTTAGCGTTAATTACGAAGATTACCTGCGTAGTAATGGGATCGAGTATGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGc                                                                     <  8:30075/149‑1 (MQ=255)
                                                                         tAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCaa                                                                   <  5:26053/148‑1 (MQ=255)
                                                                              gTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt                                                             >  1:143714/1‑149 (MQ=255)
                                                                              gTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt                                                             >  5:333702/1‑149 (MQ=255)
                                                                                aaTTACTAATATTCCCTGCGTAGTAACGGTATCGAGTATGTTATTATCCGTCACCCCTATACTGGGGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTg                                                            <  8:47155/148‑1 (MQ=25)
                                                                                  ttACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAAt                                                         <  7:219476/149‑1 (MQ=255)
                                                                                     cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                                                                    <  4:149889/119‑1 (MQ=255)
                                                                                     cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                                                                    >  3:149889/1‑119 (MQ=255)
                                                                                     cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                                                                  <  3:47126/137‑1 (MQ=255)
                                                                                     cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                                                                  >  4:47126/1‑137 (MQ=255)
                                                                                     cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTg                                                      >  5:294330/1‑149 (MQ=255)
                                                                                      gAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCTATACTGTCGGGGTGACGGAGGCTGATGACAAATTCAAAATATCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGc                                                     <  1:141675/149‑1 (MQ=255)
                                                                                       aaGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTg                                                      >  7:219354/1‑147 (MQ=255)
                                                                                       aaGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCa                                                    <  3:210446/149‑1 (MQ=255)
                                                                                       aaGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCa                                                    <  3:61890/149‑1 (MQ=255)
                                                                                          aTTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGa                                                 >  7:115467/1‑149 (MQ=255)
                                                                                             aCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACt                                               <  6:294330/148‑1 (MQ=255)
                                                                                              ccTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACttt                                             <  2:153119/149‑1 (MQ=255)
                                                                                               cTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTg                                            <  6:116114/149‑1 (MQ=255)
                                                                                                  cGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGt                                           <  2:92423/147‑1 (MQ=255)
                                                                                                  cGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTAAATCTGCAGGACTTTGTAg                                         <  8:75241/149‑1 (MQ=255)
                                                                                                      gtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCg                                     <  8:188811/149‑1 (MQ=255)
                                                                                                       taAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCg                                     >  7:150802/1‑148 (MQ=255)
                                                                                                         aaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGa                                   >  2:321172/1‑148 (MQ=255)
                                                                                                         aaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGa                                   <  2:143714/148‑1 (MQ=255)
                                                                                                           gggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGAt                                  >  5:306578/1‑147 (MQ=255)
                                                                                                              aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                             <  7:89566/149‑1 (MQ=255)
                                                                                                              aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                             >  1:319861/1‑149 (MQ=255)
                                                                                                                  aGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGtt                          >  6:301545/1‑148 (MQ=255)
                                                                                                                  aGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGtt                          >  2:3415/1‑148 (MQ=255)
                                                                                                                   gtattttaTTATCCGTCACCCCTATACTGTCGGGGTTACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTa                         <  8:37354/141‑1 (MQ=255)
                                                                                                                   gTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTaa                        >  2:188073/1‑149 (MQ=255)
                                                                                                                                 gTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCa          >  8:195189/1‑149 (MQ=255)
                                                                                                                                  tCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAt         >  8:223404/1‑149 (MQ=255)
                                                                                                                                    aCCCCGATTCTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAtt        <  4:322968/148‑1 (MQ=255)
                                                                                                                                    aCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAtt        >  1:129689/1‑148 (MQ=255)
                                                                                                                                    aCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAtt        <  8:219354/148‑1 (MQ=255)
                                                                                                                                    aCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa       <  2:259841/149‑1 (MQ=255)
                                                                                                                                      cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                                                                    >  4:83724/1‑70 (MQ=255)
                                                                                                                                      cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                                                                    <  3:83724/70‑1 (MQ=255)
                                                                                                                                      cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTaa      >  7:226559/1‑148 (MQ=255)
                                                                                                                                      cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAAc     <  6:306578/149‑1 (MQ=255)
                                                                                                                                         gATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACaaa  <  7:125930/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
CACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAA  >  NC_000913/857437‑857722

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