Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,584 A→G 6.7% D108G (GAC→GGC) ‡ mcbR → DNA‑binding transcriptional dual regulator McbR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,520,5840AG 6.7% 265.8 / 7.1 90D108G (GAC→GGC) ‡mcbRDNA‑binding transcriptional dual regulator McbR
Reads supporting (aligned to +/- strand):  ref base A (36/48);  new base G (4/2);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

GGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGC  >  NC_000913/1520435‑1520730
                                                                                                                                                     |                                                                                                                                                  
ggTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGa                                                                                                                                                    >  1:1545062/1‑150 (MQ=255)
ggTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGa                                                                                                                                                    <  1:1456444/150‑1 (MQ=255)
 gTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGAc                                                                                                                                                   <  2:92433/150‑1 (MQ=255)
          aaTGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGGTCTTGCgggg                                                                                                                                          <  1:834300/150‑3 (MQ=255)
            tGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAACTGGGGAAACGTAAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGAAAGGTCTTGCggggtg                                                                                                                                        <  2:1501406/150‑5 (MQ=255)
               gcTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAACCCTCACCCCGCAAGACCTTGCGGAACTCCa                                                                                                                                     >  2:1146960/1‑150 (MQ=255)
                cTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAg                                                                                                                                    <  2:581424/150‑1 (MQ=255)
                     tGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAAtt                                                                                                                               >  1:1587913/1‑150 (MQ=255)
                     tGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAAtt                                                                                                                               <  2:722282/150‑1 (MQ=255)
                      gTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTg                                                                                                                              >  1:832649/1‑150 (MQ=255)
                           aCCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCga                                                                                                                         <  1:277032/150‑1 (MQ=255)
                               gCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGaaa                                                                                                                     <  2:1507580/150‑1 (MQ=255)
                                 acaAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAAtt                                                                                                                   >  2:1053504/1‑150 (MQ=255)
                                 acaAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAAtt                                                                                                                   <  2:1361752/150‑1 (MQ=255)
                                 acaAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAAtt                                                                                                                   <  2:1587913/150‑1 (MQ=255)
                                  caAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTa                                                                                                                  >  1:199521/1‑150 (MQ=255)
                                    aGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAca                                                                                                                <  1:881138/150‑1 (MQ=255)
                                       aTTTACAGTTCAGGAATTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCTCTACGAACTGGAATTAATGGCAGTTGCTCTTGCTGTTGAAAACCTCACCCATCAAGACCTTGCGGAACTCAATTAATTGCTCGAGAAATTAcagca                                                                                                             <  2:810384/150‑1 (MQ=255)
                                        tttACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAcagcag                                                                                                            >  1:816428/1‑150 (MQ=255)
                                        tttACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAcagcag                                                                                                            >  1:642653/1‑150 (MQ=255)
                                         ttACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAgg                                                                                                           >  1:837454/1‑150 (MQ=255)
                                          tACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGgc                                                                                                          >  1:1631362/1‑150 (MQ=255)
                                          tACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGgc                                                                                                          <  2:642653/150‑1 (MQ=255)
                                                 ccGGAACTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAACCCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGaa                                                                                                   >  2:403174/1‑150 (MQ=255)
                                                       gTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGc                                                                                             <  1:320969/150‑1 (MQ=255)
                                                        tGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCg                                                                                            <  2:1631362/150‑1 (MQ=255)
                                                         ggggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGa                                                                                           >  2:887595/1‑150 (MQ=255)
                                                          gggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGAc                                                                                          >  1:1334175/1‑150 (MQ=255)
                                                           ggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACa                                                                                         >  1:1503576/1‑150 (MQ=255)
                                                             aaaCGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATg                                                                                       <  2:199521/150‑1 (MQ=255)
                                                                cGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGATATTACAGCAGGCGCAAGAAAAGGGCGACATGGaa                                                                                    <  2:422474/150‑1 (MQ=255)
                                                                cGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGaa                                                                                    >  2:71921/1‑150 (MQ=255)
                                                                       tAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCa                                                                             <  2:1102222/150‑1 (MQ=255)
                                                                        aGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCAt                                                                            <  2:409835/150‑1 (MQ=255)
                                                                        aGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCAGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCAt                                                                            >  1:1046409/1‑150 (MQ=255)
                                                                          aCGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTa                                                                          <  1:510086/150‑1 (MQ=255)
                                                                          aCGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTa                                                                          >  2:1369395/1‑150 (MQ=255)
                                                                           cGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTaa                                                                         >  1:928436/1‑150 (MQ=255)
                                                                           cGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTaa                                                                         >  2:1499375/1‑150 (MQ=255)
                                                                           cGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTaa                                                                         <  2:1284464/150‑1 (MQ=255)
                                                                             aaaTCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATg                                                                       <  1:800860/150‑1 (MQ=255)
                                                                              aatcaatcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGt                                                                      >  1:155098/1‑150 (MQ=255)
                                                                               atcaatcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTa                                                                     <  1:94665/150‑1 (MQ=255)
                                                                                  aatcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAAc                                                                  <  2:465473/150‑1 (MQ=255)
                                                                                   atcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACa                                                                 >  2:389014/1‑150 (MQ=255)
                                                                                      ggATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGc                                                              <  2:298330/150‑1 (MQ=255)
                                                                                          ccGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTAtt                                                          >  2:630274/1‑150 (MQ=255)
                                                                                            gCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTc                                                        <  2:155098/150‑1 (MQ=255)
                                                                                            gCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTc                                                        >  1:239587/1‑150 (MQ=255)
                                                                                            gCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTc                                                        >  2:1397860/1‑150 (MQ=255)
                                                                                            gCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTc                                                        <  2:832649/150‑1 (MQ=255)
                                                                                              tACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGc                                                      >  1:4133/1‑150 (MQ=255)
                                                                                                 gAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTa                                                   <  1:666160/150‑1 (MQ=255)
                                                                                                   aCTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGc                                                 <  1:1053504/150‑1 (MQ=255)
                                                                                                      ggAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAt                                              >  1:1247261/1‑150 (MQ=255)
                                                                                                           tAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatc                                         >  1:456025/1‑150 (MQ=255)
                                                                                                              tGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatcatc                                      <  1:1069652/150‑1 (MQ=255)
                                                                                                               ggCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCg                                     <  1:1507673/150‑1 (MQ=255)
                                                                                                                gCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAAGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGc                                    >  2:547461/1‑150 (MQ=255)
                                                                                                                   gTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCa                                 >  2:429871/1‑150 (MQ=255)
                                                                                                                      tctgtGGCTGTTGAAATCCTGACCCCGCAAGAGCTTGCGGAACTGCAGGAATTGCTTGAGAAATTACAGCATGCGCATGAATAGGTCGACATGGAACAATTCATTTATGTATACAGTCTATTTCGGTTAGCGATTTATCATCGCTCaatt                              <  2:376326/146‑3 (MQ=255)
                                                                                                                      gCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGGTATTTCGCTTAGCGATTTATCATCGCTCAAat                              <  2:816428/150‑1 (MQ=255)
                                                                                                                       ctctGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAtt                                             <  2:533768/134‑1 (MQ=255)
                                                                                                                       ctctGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAtt                                             >  1:533768/1‑134 (MQ=255)
                                                                                                                       ctctGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTAGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAata                             <  2:456025/150‑1 (MQ=255)
                                                                                                                         ctGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCa                                 <  1:1342376/144‑1 (MQ=255)
                                                                                                                         ctGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCa                                 >  2:1342376/1‑144 (MQ=255)
                                                                                                                            gCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGccc                        <  2:1503576/150‑1 (MQ=255)
                                                                                                                                ttGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATcc                    <  2:239587/150‑1 (MQ=255)
                                                                                                                                ttGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATcc                    <  2:1611414/150‑1 (MQ=255)
                                                                                                                                  gAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGACCATCCtg                  >  1:749314/1‑150 (MQ=255)
                                                                                                                                  gAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAAACATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCtg                  >  1:1165777/1‑150 (MQ=255)
                                                                                                                                   aaaaCCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCtgt                 >  1:1524023/1‑150 (MQ=255)
                                                                                                                                    aaaCCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCtgtg                <  2:1585269/150‑1 (MQ=255)
                                                                                                                                     aaCCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCtgtgt               <  2:928436/150‑1 (MQ=255)
                                                                                                                                      aCCTCACCCCGCAAGACCTTGCGGAACTCCAGTAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCtgtgtg              <  1:1149664/150‑1 (MQ=255)
                                                                                                                                      aCCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCtgtgtg              <  2:577168/150‑1 (MQ=255)
                                                                                                                                       tggagttccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGAGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTAAAATATGCCCATCCTGTGTga             >  1:1247231/8‑150 (MQ=255)
                                                                                                                                       ccTCACCCCTCAAGACCTTGCGGAACTCCAGGAATTTCTCTAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTATCGATTTATCATCGCTCAAATATGACCATCCTGTGTga             <  2:749314/150‑1 (MQ=255)
                                                                                                                                       ccTCACCCCGCAAGACCTTGCGGAACTCCAGTAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGTCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTga             <  1:440242/150‑1 (MQ=255)
                                                                                                                                       ccTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTga             <  2:521851/150‑1 (MQ=255)
                                                                                                                                       ccTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTga             >  1:1098861/1‑150 (MQ=255)
                                                                                                                                       ccTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTga             >  1:304197/1‑150 (MQ=255)
                                                                                                                                       ccTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATATCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTga             <  2:524983/150‑1 (MQ=255)
                                                                                                                                        ggagttccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTgag            >  1:999748/7‑150 (MQ=255)
                                                                                                                                         gagttccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTgaga           >  1:181511/6‑150 (MQ=255)
                                                                                                                                          cACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgat          <  1:1623931/150‑1 (MQ=255)
                                                                                                                                            ccccGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCACAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgatga        >  2:159296/1‑150 (MQ=255)
                                                                                                                                              ccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGAtt      >  2:99167/1‑150 (MQ=255)
                                                                                                                                               cGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTg     <  1:1016253/150‑1 (MQ=255)
                                                                                                                                               cGCAAGACCTTGCGGAACTCCACGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTg     <  1:547461/150‑1 (MQ=255)
                                                                                                                                                  aaGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGc  <  1:887595/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAGC  >  NC_000913/1520435‑1520730

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