Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,127,252 A→G 7.6% intergenic (+27/+34) rpmE → / ← yiiX 50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,127,2520AG 7.6% 287.5 / 10.9 106intergenic (+27/+34)rpmE/yiiX50S ribosomal subunit protein L31/putative lipid binding hydrolase, DUF830 family protein
Reads supporting (aligned to +/- strand):  ref base A (45/52);  new base G (4/4);  total (49/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

Only 100 of 107 total aligned reads displayed.

GACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTC  >  NC_000913/4127103‑4127399
                                                                                                                                                      |                                                                                                                                                   
gACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCa                                                                                                                                                     <  2:1618091/150‑1 (MQ=255)
 aCCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCAc                                                                                                                                                    >  2:199629/1‑150 (MQ=255)
   cTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACgg                                                                                                                                                  <  2:1454957/150‑1 (MQ=255)
    ggAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGgc                                                                                                                                                 <  2:1512775/149‑1 (MQ=255)
     gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGgcg                                                                                                                                                <  1:449962/150‑1 (MQ=255)
     gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGgcg                                                                                                                                                <  1:107878/150‑1 (MQ=255)
     gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGgcg                                                                                                                                                <  2:33548/150‑1 (MQ=255)
     gAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGgcg                                                                                                                                                <  1:1578793/150‑1 (MQ=255)
        ccTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCtt                                                                                                                                             >  1:1632129/1‑150 (MQ=255)
        ccTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCtt                                                                                                                                             >  1:1607294/1‑150 (MQ=255)
        ccTCGACGTGTGCAGCAAGTGCCACCCGTTCATCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTACAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGgcgcat                                                                                                                                             >  2:1036579/1‑148 (MQ=255)
         cTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCttt                                                                                                                                            <  1:904765/150‑1 (MQ=255)
              cGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTtg                                                                                                                                       >  2:1266727/1‑150 (MQ=255)
               gtgtGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGGTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCGAAAAAAAGCGCCGTGCGGCGCttttttcg                                                                                                                                       <  1:1079539/150‑3 (MQ=255)
               gtgtGCAGCAAGGGCCACCCGTGCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCA‑AAAAAAGCGCCGTGCGGCGCttttttcgg                                                                                                                                      <  1:1211105/150‑4 (MQ=255)
                  tgCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGcc                                                                                                                                   >  1:1387443/1‑150 (MQ=255)
                   gcagcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGGGGCCGGGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGAATTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCt                                                                                                                                  <  1:103885/150‑1 (MQ=255)
                       caAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAt                                                                                                                              >  1:636457/1‑150 (MQ=255)
                         aGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAttt                                                                                                                            >  2:471817/1‑150 (MQ=255)
                          gTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAtttt                                                                                                                           <  2:365006/150‑1 (MQ=255)
                            gCCACCCGTTCTTCACTGGCATACAGCGTGATGTTGCTACCGGTGGACGTGTTGACCGCTTCAACAAGCGTTTCATCATCCCGGGCAGCAAATAAGACCGGATTTCCG‑ATAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTgg                                                                                                                         <  1:277828/150‑1 (MQ=255)
                                  cGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCGTTTATTGTGCCTGAATTTTGGTCGTAt                                                                                                                   >  1:956529/1‑150 (MQ=255)
                                   gTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTAtt                                                                                                                  <  1:456963/150‑1 (MQ=255)
                                    ttcttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTa                                                                                                                 <  1:508228/150‑1 (MQ=255)
                                      cttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTAcc                                                                                                               <  1:186484/150‑1 (MQ=255)
                                       ttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCa                                                                                                              <  2:636457/150‑1 (MQ=255)
                                       ttcACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCa                                                                                                              >  1:1102728/1‑150 (MQ=255)
                                           cTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGCCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGaa                                                                                                          >  2:620587/1‑150 (MQ=255)
                                             ggCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGaaaa                                                                                                        >  2:794373/1‑150 (MQ=255)
                                              gCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAAc                                                                                                       <  2:599947/150‑1 (MQ=255)
                                               cAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAAAAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATATCCT‑AAAAAAGCGCCCCACGGCGCTTTTTTAGTGACAGAATTTTGGTCGTATTACCACGaaaaaa                                                                                                      >  1:589007/1‑148 (MQ=255)
                                                   cAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCgg                                                                                                  <  1:994621/150‑1 (MQ=255)
                                                    aGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGc                                                                                                 <  2:520844/150‑1 (MQ=255)
                                                     gCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGcc                                                                                                >  2:665327/1‑150 (MQ=255)
                                                             ttGCTACCGGGGGCCGTGTTGACCGCTTCAACTAGCGATTCAACATCCCGGGCAGCAAATAAGACCGGAGTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTg                                                                                        <  2:1004426/150‑1 (MQ=255)
                                                              tGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTgg                                                                                       <  2:37133/150‑1 (MQ=255)
                                                              tGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTgg                                                                                       <  1:1288232/150‑1 (MQ=255)
                                                              tGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTgg                                                                                       >  2:1140541/1‑150 (MQ=255)
                                                               gCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGc                                                                                      <  1:1140541/150‑1 (MQ=255)
                                                                cTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCg                                                                                     >  1:1603537/1‑150 (MQ=255)
                                                                cTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCg                                                                                     >  1:622363/1‑150 (MQ=255)
                                                                   ccGGTGGCCGTGTTGACCGCTTCTACAAGCGTTTCAACATCCCGGGCAGCAAATAAGTCCGGATTTCCG‑AAAATAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGAcc                                                                                  <  2:1669236/150‑1 (MQ=255)
                                                                       tGGCCGTTTTGACCGCTTCAACAAGCGTTTCATCATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGTCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtgg                                                                              <  1:1011646/150‑1 (MQ=255)
                                                                       tGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtgg                                                                              <  1:199629/150‑1 (MQ=255)
                                                                         gCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtggtg                                                                            <  2:1239749/150‑1 (MQ=255)
                                                                         gCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCgtggtg                                                                            <  1:1422998/150‑1 (MQ=255)
                                                                         gCCGTGTTGACCGCTTCAACAAGCGTTGCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCatggtg                                                                            <  2:6520/150‑1 (MQ=255)
                                                                           cGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAg                                                                          <  2:1207605/150‑1 (MQ=255)
                                                                            gtgtTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGt                                                                         >  2:262005/1‑150 (MQ=255)
                                                                            gtgtTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGt                                                                         <  1:1036579/150‑1 (MQ=255)
                                                                                tGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCAGTCTTTGGCGACCGTGGTGAGTtgtg                                                                     >  1:1193276/1‑150 (MQ=255)
                                                                                 gACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTgg                                                                    >  2:1545773/1‑150 (MQ=255)
                                                                                 gACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTgg                                                                    >  2:175392/1‑150 (MQ=255)
                                                                                  aCCGCTTCAACAAGCGTTTCATCATGCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGgc                                                                   <  2:853182/150‑1 (MQ=255)
                                                                                   ccGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGgcg                                                                  <  2:225486/150‑1 (MQ=255)
                                                                                     gCTGCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCa                                                                <  2:1210923/150‑1 (MQ=255)
                                                                                      cTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCaa                                                               >  1:480617/1‑149 (MQ=255)
                                                                                       ttCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATc                                                              <  2:768524/150‑1 (MQ=255)
                                                                                        tCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCg                                                             <  2:773292/150‑1 (MQ=255)
                                                                                        tCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTACG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTTCCACGAAAACATCGTCCATTCTTTGGCGACCTTGGTGAGTGGTCGCGAATCg                                                             >  2:324843/1‑150 (MQ=255)
                                                                                            caaGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGaaaa                                                         >  1:732724/1‑150 (MQ=255)
                                                                                              aGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAAGTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACg                                                       >  1:89494/1‑150 (MQ=255)
                                                                                                cGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGc                                                     >  1:939273/1‑150 (MQ=255)
                                                                                                  tttCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGGATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCtt                                                   >  2:151343/1‑150 (MQ=255)
                                                                                                      aaCATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCTCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTgg                                               <  1:1546395/150‑1 (MQ=255)
                                                                                                      aaCATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTgg                                               <  1:809847/150‑1 (MQ=255)
                                                                                                        cATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCg                                             >  1:946725/1‑150 (MQ=255)
                                                                                                         aTCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCga                                            <  1:638003/150‑1 (MQ=255)
                                                                                                            ccGGTCAGCAAATAAGACCTGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGac                                         <  1:400951/150‑1 (MQ=255)
                                                                                                            ccGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGac                                         >  1:1476653/1‑150 (MQ=255)
                                                                                                                  agcaAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTTGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGt                                   <  2:1168721/150‑1 (MQ=255)
                                                                                                                     aaaTAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTc                                <  1:175392/150‑1 (MQ=255)
                                                                                                                      aaTAAGACCTGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGCCCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTcc                               <  2:249943/150‑1 (MQ=255)
                                                                                                                      aaTAAGAACGGATTTCCG‑AAAAAAGCGCCGCGCGGGGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGCCAACAGCGGCGACTGTTTGGCGCCCGTGGGGCGGAGTGGCGCATCGAAACCGCCTTGTGGCGAGACGCCCGTTTcc                               >  1:1067434/1‑150 (MQ=255)
                                                                                                                       aTAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTGTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAACACGGCTTGTGGCGAGACCCCCGTTTccg                              >  1:1511496/1‑149 (MQ=255)
                                                                                                                        tAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTc                             >  1:1540784/1‑150 (MQ=255)
                                                                                                                           gACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGGCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAg                          <  1:858648/150‑1 (MQ=255)
                                                                                                                           gACAGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAg                          >  2:1117195/1‑150 (MQ=255)
                                                                                                                            aCCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGc                         >  2:664717/1‑150 (MQ=255)
                                                                                                                             ccGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCg                        <  2:797301/150‑1 (MQ=255)
                                                                                                                              cGGATTTCCG‑AAAACAGCGCCGCACGGCGCTTTTTTTGTGCCTTAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCgg                       <  1:1027830/150‑1 (MQ=255)
                                                                                                                               ggATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGa                      >  1:1660742/1‑150 (MQ=255)
                                                                                                                               ggATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGa                      <  2:1387443/150‑1 (MQ=255)
                                                                                                                                 aTTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAat                    >  1:476106/1‑150 (MQ=255)
                                                                                                                                      cG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCt               <  2:326142/150‑1 (MQ=255)
                                                                                                                                      cG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCt               >  1:1356134/1‑150 (MQ=255)
                                                                                                                                      cG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCt               >  1:1610668/1‑150 (MQ=255)
                                                                                                                                      acaaaaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTc                >  1:1164593/4‑150 (MQ=255)
                                                                                                                                         aaaaaaGCGCCGTGCGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTg             >  2:1626102/1‑150 (MQ=255)
                                                                                                                                          aaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGc            <  1:404911/150‑1 (MQ=255)
                                                                                                                                          aaaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGc            >  1:1045223/1‑150 (MQ=255)
                                                                                                                                           aaaaGCGCCGCTCGGCGCTTTTTTGGTGCCTTAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAATCGGCTTGTGGCGAGACCACCGTTTCCTCAATCGGAATATTCTTGcc           <  2:589007/150‑1 (MQ=255)
                                                                                                                                           aaaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGcc           >  1:770497/1‑150 (MQ=255)
                                                                                                                                            aaaGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCg          <  1:1545773/150‑1 (MQ=255)
                                                                                                                                              aGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTa        >  1:595827/1‑150 (MQ=255)
                                                                                                                                               gcgcCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACAGTTTCCTCAAGCGGAATATTCTTGCCGTaa       <  1:471817/150‑1 (MQ=255)
                                                                                                                                                cgccgcACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAAc      <  1:1700649/150‑1 (MQ=255)
                                                                                                                                                 gccgcACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACg     >  2:795485/1‑150 (MQ=255)
                                                                                                                                                    gcACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTc  >  1:1270864/1‑150 (MQ=255)
                                                                                                                                                    gcACGGCGCTTTTTTTGAGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCAATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGTCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTc  <  1:261243/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGATTTCCG‑AAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGAATTTTGGTCGTATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTC  >  NC_000913/4127103‑4127399

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