Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,511,139 G→C 9.5% W379C (TGG→TGC ydcR → fused putative DNA‑binding transcriptional regulator/putative aminotransferase YdcR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,511,1390GC 9.5% 222.8 / 11.2 71W379C (TGG→TGCydcRfused putative DNA‑binding transcriptional regulator/putative aminotransferase YdcR
Reads supporting (aligned to +/- strand):  ref base G (39/25);  new base C (2/5);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGTAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTA  >  NC_000913/1510990‑1511287
                                                                                                                                                     |                                                                                                                                                    
cGTAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTgg                                                                                                                                                      >  1:1355327/1‑150 (MQ=255)
 gTAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTggc                                                                                                                                                     >  2:763852/1‑150 (MQ=255)
 gTAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTCCATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCACCTTGCCGAACGTAAACAACCTGCCTggc                                                                                                                                                     >  2:167155/1‑150 (MQ=255)
 gTAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTAACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTggc                                                                                                                                                     >  2:161982/1‑150 (MQ=255)
      aTTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTggcaggca                                                                                                                                                >  1:1088610/1‑150 (MQ=255)
         cAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCActg                                                                                                                                             >  1:983301/1‑150 (MQ=255)
           aCGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGCCAGGCAcgttg                                                                                                                                           <  1:1020349/150‑5 (MQ=255)
           aCGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGCCAGGCAcgttg                                                                                                                                           <  2:928094/150‑5 (MQ=255)
            cGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCActgctg                                                                                                                                          <  2:359637/150‑1 (MQ=255)
            cGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGCCAGGCAcgttgt                                                                                                                                          <  2:183608/150‑6 (MQ=255)
              cTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCg                                                                                                                                        <  1:330584/150‑1 (MQ=255)
              cTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGCCAGGCAcgttgttt                                                                                                                                        <  1:1188513/150‑8 (MQ=255)
              cTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGCCAGGCAcgttgttt                                                                                                                                        <  1:1193466/150‑8 (MQ=255)
                 gCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTa                                                                                                                                     <  1:1019622/150‑1 (MQ=255)
                 gCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTa                                                                                                                                     >  1:1113108/1‑150 (MQ=255)
                   aGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATc                                                                                                                                   >  1:1301606/1‑150 (MQ=255)
                      tgatgaGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATctgc                                                                                                                                >  1:1057858/1‑150 (MQ=255)
                       gatgaGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATctgcc                                                                                                                               >  2:833572/1‑150 (MQ=255)
                              aCGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGaa                                                                                                                        >  1:178191/1‑150 (MQ=255)
                               cGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAg                                                                                                                       <  2:202287/150‑1 (MQ=255)
                               cGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAg                                                                                                                       >  2:357890/1‑150 (MQ=255)
                                gCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGt                                                                                                                      <  2:1355327/150‑1 (MQ=255)
                                       accaGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAAtt                                                                                                               >  1:167501/1‑150 (MQ=255)
                                       accaGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAAtt                                                                                                               >  1:80578/1‑150 (MQ=255)
                                          agatCACCTATTCAACTTGAGATTGTCGATTACATTTCCACGCGCCGATACCACGCACATAATCGTCGCCTGCCTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATtcat                                                                                                            >  2:1268177/4‑150 (MQ=255)
                                            cTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATtcatca                                                                                                          >  2:684434/1‑150 (MQ=255)
                                              cACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATa                                                                                                        >  2:1107552/1‑150 (MQ=255)
                                                    tGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACa                                                                                                  <  2:167501/150‑1 (MQ=255)
                                                     gCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCAAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATATGACTGCGGAAGTGAAAATTCATCATAATGACAg                                                                                                 >  1:516654/1‑150 (MQ=255)
                                                      cAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGt                                                                                                >  1:1044690/1‑150 (MQ=255)
                                                            gcgcTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTAc                                                                                          <  2:637327/150‑1 (MQ=255)
                                                            gcgcTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTAc                                                                                          >  1:887844/1‑150 (MQ=255)
                                                                tggtggATTAGCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTtc                                                                                      <  1:161982/150‑1 (MQ=255)
                                                                tggtggATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCTGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTtc                                                                                      <  1:167155/150‑1 (MQ=255)
                                                                     gATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTTCTTTCTCTgg                                                                                 >  2:1320679/1‑150 (MQ=255)
                                                                          ccTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGa                                                                            <  1:155741/150‑1 (MQ=255)
                                                                                    cgcCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAg                                                                  >  2:117565/1‑150 (MQ=255)
                                                                                    cgcCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAg                                                                  >  2:210745/1‑150 (MQ=255)
                                                                                    cgcCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAg                                                                  <  2:143523/150‑1 (MQ=255)
                                                                                      ccGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGcc                                                                >  2:1312681/1‑150 (MQ=255)
                                                                                        gATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGt                                                              >  1:225934/1‑150 (MQ=255)
                                                                                         aTACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGtt                                                             <  1:654808/150‑1 (MQ=255)
                                                                                          tACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTa                                                            <  1:763852/150‑1 (MQ=255)
                                                                                           acgaccCCCATATTCGTCGACTGCGTCCCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGTCACTGCTGCGTTATCTTCCTGCGGAAGTGAAAATTCATCAAAATGAAAGTGGTTACTATCTCTGGTTGGAGCTCCCCGAGCCGTTAg                                                           >  2:620992/1‑150 (MQ=255)
                                                                                            cgacgCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGa                                                          >  1:244331/1‑150 (MQ=255)
                                                                                                 cccATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACCTGCCTGGCAGGCACTGCTGCGTTATATGCCTGCGGAAGTGAAAATTAATAATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGACCCGTTAGATGCCg                                                     >  1:1017954/1‑150 (MQ=255)
                                                                                                  ccATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCgg                                                    <  1:357890/150‑1 (MQ=255)
                                                                                                   cATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACTGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGc                                                   >  2:589907/1‑150 (MQ=255)
                                                                                                   cATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGc                                                   <  2:957769/150‑1 (MQ=255)
                                                                                                      cTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCgccagg                                                >  2:611251/1‑145 (MQ=255)
                                                                                                        tcgtcgCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAAtt                                              >  1:358702/1‑150 (MQ=255)
                                                                                                            cgCCTCCGTCGCCAGCTAGCGCAACGTAAACAACGTGCCTGCCAGGCACTGCTGCGTTATCTGCCTGCGTAAGTGAAAATTCATCATAATGACACAGGTTACATTCTCTGGTTGGAGCTCCCCGAGCCGTTACATGCCGGCGAATTAAGc                                          >  2:87336/1‑150 (MQ=255)
                                                                                                              ccTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCt                                        <  1:1107552/150‑1 (MQ=255)
                                                                                                                tGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTgg                                      >  1:1276353/1‑150 (MQ=255)
                                                                                                                 gCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTggc                                     >  1:128997/1‑150 (MQ=255)
                                                                                                                      gCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCAc                                >  1:528838/1‑150 (MQ=255)
                                                                                                                       ccAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACt                               <  2:1301606/150‑1 (MQ=255)
                                                                                                                        cAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTg                              <  1:1268177/150‑1 (MQ=255)
                                                                                                                        cAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTg                              >  1:537816/1‑150 (MQ=255)
                                                                                                                          gCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGAc                            >  1:105124/1‑150 (MQ=255)
                                                                                                                               cGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGcatc                       <  2:1057858/150‑1 (MQ=255)
                                                                                                                                ggAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGcatca                      >  1:648192/1‑150 (MQ=255)
                                                                                                                                 gAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGcatcat                     >  1:174924/1‑150 (MQ=255)
                                                                                                                                  aaCGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACTGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCata                    <  1:589907/150‑1 (MQ=255)
                                                                                                                                    cGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATc                  <  2:1088610/150‑1 (MQ=255)
                                                                                                                                    cGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATc                  >  1:570137/1‑150 (MQ=255)
                                                                                                                                    cGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTACTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATc                  <  2:80578/150‑1 (MQ=255)
                                                                                                                                     gTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCa                 <  2:225934/150‑1 (MQ=255)
                                                                                                                                      tAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAg                <  2:178191/150‑1 (MQ=255)
                                                                                                                                         cagcagTGCCTGCCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTAt             >  1:1303032/6‑150 (MQ=255)
                                                                                                                                          caacGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTAtt            >  1:1228753/1‑150 (MQ=255)
                                                                                                                                             cGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTgcg         >  2:5184/1‑150 (MQ=255)
                                                                                                                                               tGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGAt                                                         >  2:900049/1‑100 (MQ=255)
                                                                                                                                               tGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGAt                                                         <  1:900049/100‑1 (MQ=255)
                                                                                                                                                  cTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCggg    >  1:123178/1‑150 (MQ=255)
                                                                                                                                                  cTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCggg    >  1:1114644/1‑150 (MQ=255)
                                                                                                                                                    ggcaggcaCTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTa  >  1:678296/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGTAAAATTCAACGCTTGCAGTTGATGAGTACGCTTTCCACCAGCTCACCGATGCAACTTGCGCTGGTGGATTACCTTTCCACGCGCCGATACGACGCCCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTA  >  NC_000913/1510990‑1511287

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