Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,059,703 A→T 6.2% intergenic (‑2/‑50) yqfE ← / → argP pseudogene, LysR family/transcriptional regulator for arginine transport and DNA replication genes; replication initiation inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,059,7030AT 6.2% 224.4 / 4.9 81intergenic (‑2/‑50)yqfE/argPpseudogene, LysR family/transcriptional regulator for arginine transport and DNA replication genes; replication initiation inhibitor
Reads supporting (aligned to +/- strand):  ref base A (35/41);  new base T (2/3);  total (37/44)
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.22e-01

TAACGCCACGCCATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTACACAATC  >  NC_000913/3059558‑3059850
                                                                                                                                                 |                                                                                                                                                   
tAACGCCCCGCCATTCAAGACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTACCGCAAATCCCTGATCTTTTAAGCTGCTCGTTTAACCACGGCAACAAATGTTACCCGCCGCCGGCTGCGGTATAAGCAACCTGACACaaaaat                                                                                                                                                  >  1:1166236/1‑147 (MQ=255)
tAACGCCACGCCATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAAtt                                                                                                                                                  <  1:166019/149‑1 (MQ=255)
  acgccacgccATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgt                                                                                                                                                <  1:20107/149‑1 (MQ=255)
  acgccacgccATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgt                                                                                                                                                <  1:905334/149‑1 (MQ=255)
       acgccaTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCTAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                                                                           <  1:636468/149‑2 (MQ=255)
       acgccaTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                                                                           <  2:22390/149‑2 (MQ=255)
          ccaTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGt                                                                                                                                        <  2:717122/149‑1 (MQ=255)
          ccaTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTATAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCCGACACAATTTTGTGTCaggtt                                                                                                                                        <  1:181828/149‑5 (MQ=255)
             ttCAACACAACCTGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCa                                                                                                                                     >  1:417987/1‑149 (MQ=255)
             ttCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCa                                                                                                                                     <  2:1149904/149‑1 (MQ=255)
                aaCACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGa                                                                                                                                  <  2:110721/149‑1 (MQ=255)
                  cacaACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaa                                                                                                                                <  1:227923/149‑1 (MQ=255)
                    caACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATTGTGCAGGataaa                                                                                                                              >  2:1123345/1‑149 (MQ=255)
                    caACCAGGGCGAGAACTCGCTGAGAGAAATGAATTATCGAAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaaaa                                                                                                                              <  1:1215201/149‑1 (MQ=255)
                     aaCCAGGGCGCGAACTCGCTGACCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTATCCACGGATACCAATGTTACCCTCCAACGGCTTCGGAATATGCAACCTGACAGAAAATAGTGTCATAGTGAAGGAAAAAg                                                                                                                             >  1:1029821/1‑149 (MQ=255)
                        cctggCGCGAACTCGCTGAGCGTAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAt                                                                                                                          <  2:791909/146‑1 (MQ=255)
                        cAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAt                                                                                                                          <  1:59574/149‑1 (MQ=255)
                        cAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCAGACGGCTTCGGTATATGCATCCTGACACAAAATTGTGTCATAGTGCAGGAATAAGCAt                                                                                                                          >  1:1144067/1‑149 (MQ=255)
                         cggccgcgAACTCGCTGAGCGAAATGCATTAGCGCCACTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAtt                                                                                                                         <  2:1420617/145‑1 (MQ=255)
                         aGGGCGCGAACTCTCTGAGCGAAATGCATTAGCGCAAATCACTGATCTTTTAAGCTGCGCTTATAACCAAGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCAGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAtt                                                                                                                         >  1:911308/1‑149 (MQ=255)
                         aGGGCGCGAACTCGCTGAGCGAAATGCGTTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTTACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCACCCTGACCCAAAATTGTGGCATAGTGCAGGAAAAAGCttt                                                                                                                         >  1:1158917/1‑149 (MQ=255)
                            gcgcGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAc                                                                                                                      <  2:248006/149‑1 (MQ=255)
                            gcgcGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAc                                                                                                                      >  2:366177/1‑149 (MQ=255)
                                 aaCTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGa                                                                                                                 <  1:406433/149‑1 (MQ=255)
                                  aCTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAg                                                                                                                <  1:364117/149‑1 (MQ=255)
                                  aCTCGCGGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAAACACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAg                                                                                                                <  2:1193732/149‑1 (MQ=255)
                                         gAGCGAAATGCATTAGCGGAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTCTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaa                                                                                                         <  1:616940/149‑1 (MQ=255)
                                           gCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCGTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaaca                                                                                                       <  2:137119/149‑1 (MQ=255)
                                            cGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACATAATTGTGTAATAGTGCAGGAAAAAGCATTTACAAGGAGCAGACAACAg                                                                                                      >  2:364292/1‑149 (MQ=255)
                                               aaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                   <  1:425326/149‑1 (MQ=255)
                                               aaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                   <  1:4260/149‑1 (MQ=255)
                                               aaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                   >  1:1215779/1‑149 (MQ=255)
                                                                    gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGaa                                                                              >  1:10619/1‑149 (MQ=255)
                                                                        ttttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGCGTCAAAATTGTGTCATAGTGCAGGAAAAAGCTTTTACCAGGAGCAGCCAACAGCAATGAAACGCCCGGACGACAGCACAtt                                                                          >  1:328719/1‑149 (MQ=255)
                                                                         tttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTa                                                                         <  1:1189560/149‑1 (MQ=255)
                                                                         tttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATGTACCAGGAGCAGCCAACAGCAATGAAACGCCCGGACTACAGAACAttc                                                                         >  2:330179/1‑148 (MQ=255)
                                                                          ttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                        <  1:366177/149‑1 (MQ=255)
                                                                            aaGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAg                                                                      >  1:1216602/1‑149 (MQ=255)
                                                                            aaGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAg                                                                      <  2:925378/149‑1 (MQ=255)
                                                                              gCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGc                                                                    >  1:16960/1‑149 (MQ=255)
                                                                                tGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAATTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCAc                                                                  <  1:675586/149‑1 (MQ=255)
                                                                                     ttttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGAt                                                             >  1:1263953/1‑149 (MQ=255)
                                                                                     ttttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGAt                                                             >  2:76345/1‑149 (MQ=255)
                                                                                     ttttAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGAACTACAGAACATTACAGGCACTGGAt                                                             >  2:583986/1‑149 (MQ=255)
                                                                                       ttAACCACGGATTCAAATGTTACCCACCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAATTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGc                                                           <  2:338251/149‑1 (MQ=255)
                                                                                         aaCCACGGATTCAAATGTTACCCGCCGACGGATTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg                                                         <  1:1043891/149‑1 (MQ=255)
                                                                                         aaCCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg                                                         >  1:681662/1‑149 (MQ=255)
                                                                                         aaCCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg                                                         <  1:375892/149‑1 (MQ=255)
                                                                                         aaCCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg                                                         >  2:589843/1‑149 (MQ=255)
                                                                                           ccACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACAAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCggag                                                       >  2:1303903/1‑147 (MQ=255)
                                                                                           ccACGGATACAAATGTTACCCGCCGACGGCGTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg                                                       <  2:1100865/149‑1 (MQ=255)
                                                                                            cACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGa                                                      <  1:818442/149‑1 (MQ=255)
                                                                                             aaggATACAAATGTTACCAGCAGAAGTCTTGTGTATATGAAACCTGACCCAAACTTTTGTAATATTGTAGTAAAAAGCATTTACAATGTGCAGACACCAGCAATGAAACGCCCGGACTACATAACATTACATGCACTGGATGCGGTGAt                                                     <  2:1029821/147‑1 (MQ=255)
                                                                                                     aaaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCCGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGaaa                                             >  1:870020/1‑148 (MQ=255)
                                                                                                     aaaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACATCAGCAATGAAACGCCCGGACTACAGAACAGTACAGGCACTGGATGAGGTGTGACGAGAAc                                             >  1:196216/1‑149 (MQ=255)
                                                                                                     aaaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAAc                                             <  1:646229/149‑1 (MQ=255)
                                                                                                     aaaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTCTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATCAAACGCCCCGACTACAGAACATTACAGGCACTGGATGCCGTGATACGTGAAc                                             >  2:217596/1‑149 (MQ=255)
                                                                                                       aTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGa                                           <  2:970542/149‑1 (MQ=255)
                                                                                                         gTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGCCAACAGCAATGAAACGCCCGGAATACAGAACATTACAGGCACAGGATGCGGTGATACGTGAACGAgg                                         >  2:1205364/1‑149 (MQ=255)
                                                                                                          ttACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGa                                        <  1:1354569/149‑1 (MQ=255)
                                                                                                          ttACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACAAGGAGCAGACAACAGCAATGATACGCCCGGACTACAGAACATTACAGGGACTGGATGCGGTGCTACGTGAACGAGGa                                        >  2:1057291/1‑149 (MQ=255)
                                                                                                           tACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGAt                                       <  1:589843/149‑1 (MQ=255)
                                                                                                               cgccgACGGCTTCGGTATTTGCTAACTGACACAAAATTGTGTCTTAGTGCAGTAAAAAGCATTTACGAGGAGGTGACAACAGGAATGAAACGCCCGTACTACAGATCATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGa                                   <  1:988649/149‑1 (MQ=255)
                                                                                                               cgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGa                                   <  1:699672/149‑1 (MQ=255)
                                                                                                                gccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAg                                  <  1:698159/149‑1 (MQ=255)
                                                                                                                     cGGCTTCGGTATATGCAACCTGACACAAAATGGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAgcgcgc                             <  2:890091/149‑1 (MQ=255)
                                                                                                                      ggCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGcgcgcg                            >  2:928904/1‑149 (MQ=255)
                                                                                                                       gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCgg                           <  2:681662/149‑1 (MQ=255)
                                                                                                                        cTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGTATTTGAGCGCGCGGc                          >  2:150823/1‑149 (MQ=255)
                                                                                                                        cTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGc                          >  1:1481342/1‑149 (MQ=255)
                                                                                                                         ttCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACAAGGAGCAGACAACAGCAATGAAACTCGCGGACTACATAACATTACAGGCACTGGATGCGGTGATACGTGAAAGAGGATTTTAGCGCGCGGca                         >  2:854316/1‑149 (MQ=255)
                                                                                                                             gTATATGCAACCTGACACAAAATTGGGGCATAGTGAAGGAAACAGCATGTACCCGGAGCAGACAACACCAAGGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGGGAACGAGGAATTGAGCGCGCGGCACaaa                     <  2:1172319/149‑1 (MQ=255)
                                                                                                                              tataTGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACACCAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACaaaa                    <  2:431777/149‑1 (MQ=255)
                                                                                                                                 atGCAACCTGACACAAAATTGTCTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCt                 <  2:203147/149‑1 (MQ=255)
                                                                                                                                  tGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGCAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGAAACGTGAACGAGGATTTGATCGCGCGGCACAAAAGCtg                >  1:379606/1‑149 (MQ=255)
                                                                                                                                      aCCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCa            >  1:291935/1‑149 (MQ=255)
                                                                                                                                       tatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGCACATTACAGGCACTGTATGCGGTGATACGAGAACGAGGATTTGAGCGCGCGGCACTAACGCTGTGCAt           >  1:1051741/3‑149 (MQ=255)
                                                                                                                                       tatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCAt           >  1:1003385/3‑149 (MQ=255)
                                                                                                                                        cTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCAtt          >  1:509256/1‑149 (MQ=255)
                                                                                                                                           acacAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCCCTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTaca       >  2:1277825/1‑149 (MQ=255)
                                                                                                                                             acaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTAcaca     >  2:200778/1‑149 (MQ=255)
                                                                                                                                             acaAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGCCACAAAAGCTGTGCATTAcaca     >  1:648922/1‑149 (MQ=255)
                                                                                                                                               aaaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTACACAAt   <  1:76345/149‑1 (MQ=255)
                                                                                                                                                aaaTTGTTTCATAGTGCAGGAAAAACCATTTTCCATGTGAAGACATCAGAAATGAAACTCCCGGACTACATAACATTACAGTCAATGGATGCGGTGATACGTGAACGAGGATTTGAGAGCGCGGCACAAAAGCTGTGCATTACAGAATc  <  2:1331677/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TAACGCCACGCCATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTACACAATC  >  NC_000913/3059558‑3059850

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