Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,059,704 A→T 8.3% intergenic (‑3/‑49) 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,7040AT 8.3% 180.9 / 2.1 70intergenic (‑3/‑49)yqfE/argPpseudogene, LysR family/transcriptional regulator for arginine transport and DNA replication genes; replication initiation inhibitor
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base T (5/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCACGCCATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTACACAAT  >  NC_000913/3059561‑3059849
                                                                                                                                               |                                                                                                                                                 
cgccacgccaTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                              <  2:1020616/149‑1 (MQ=255)
cgccacgccaTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                              <  1:1138771/149‑1 (MQ=255)
  ccacgccaTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTc                                                                                                                                            <  1:634843/149‑1 (MQ=255)
     cgccaTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGAACTATTAAGCTGAGAATATACCCACAGATAAAAATGTTACCCGCCGACCGATTCGGTATATTCAACCTGACAAAAAAGTGTTTCATa                                                                                                                                         >  1:1149477/1‑149 (MQ=255)
      gccaTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCaggt                                                                                                                                        <  1:893273/149‑4 (MQ=255)
        caTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTg                                                                                                                                      <  1:395871/149‑1 (MQ=255)
        caTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTg                                                                                                                                      <  1:594893/149‑1 (MQ=255)
               cacaACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaa                                                                                                                               >  1:748931/1‑149 (MQ=255)
               cacaACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaa                                                                                                                               >  2:685/1‑149 (MQ=255)
               cacaACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTGCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaa                                                                                                                               <  1:50373/149‑1 (MQ=255)
                acaACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGCCGGCTTCGGTATATGCACCCTGACACAAAAATGTGTCATAGTGCAGGaaaa                                                                                                                              >  2:460660/1‑149 (MQ=255)
                  aaacAGGGCGCTAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACGACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                                                                            <  2:328746/146‑1 (MQ=255)
                  aaCCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGGGTCATAGTGCAGGAAAAAg                                                                                                                            >  2:1142314/1‑149 (MQ=255)
                      aGGGCGCGAACTCGCTGAGCGAAATGCATTATCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGCATACAAATGTTACCCGCCGACGGCTTCGGTATATACAACCTGAAACAAAATTGTGTCATAGTGCAGGAAAAAGCAtt                                                                                                                        >  2:365906/1‑149 (MQ=255)
                         gtgcgAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGGTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTGCGGTATTTGCAACCTGACACAAAATTGTGGCATAGTGCAGGAAAACGCAGTTAc                                                                                                                     >  2:219954/3‑149 (MQ=255)
                         gcgcGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAc                                                                                                                     >  1:511078/1‑149 (MQ=255)
                           gcgAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                                                                   <  1:730020/149‑1 (MQ=255)
                            cgAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAg                                                                                                                  <  1:368268/149‑1 (MQ=255)
                              aaCTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGa                                                                                                                <  1:272876/149‑1 (MQ=255)
                                cTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                                                                              <  1:1176276/149‑1 (MQ=255)
                                 tCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCa                                                                                                             <  1:22105/149‑1 (MQ=255)
                                     tGAGCGAAATGCATTAGCGCAATTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAGAGTGCAGGAAAAAGCATTTACCAGGAGCAGaca                                                                                                         >  1:920102/1‑149 (MQ=255)
                                     tGAGCGAAATGCATTAGCGCAACTCACTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAATTGTTACCCGCCGCACGATTCGGTATATGAAGCCTGCCAGAAATATGTCTCAGAGAGCAGGAAACAGGATTTACCAGGCGCCGaca                                                                                                         >  2:816634/1‑149 (MQ=255)
                                            aaTGCATTAGACCAAATCCCTTATCTATTAAGCTGCGCTTTTAACCACGGAAACAAATGTTACACGCCGACGGCTTCTGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                  <  2:769387/149‑1 (MQ=255)
                                              tGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGAAGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGCGCAGACAACAGCAATg                                                                                                <  1:892369/149‑1 (MQ=255)
                                                cATTAGCGCAAATCCCTGATCTTTTAAGCTGCGGTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGAAACCTGACACAAATTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                                                                              <  2:402071/149‑1 (MQ=255)
                                                cATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                                                                              <  1:588702/149‑1 (MQ=255)
                                                cATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                                                                              >  1:1131314/1‑149 (MQ=255)
                                                 aTTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCCTTTACCAGGAGCAGACAACAGCAATGaaa                                                                                             >  2:147901/1‑149 (MQ=255)
                                                 aTTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaaa                                                                                             <  1:1171609/149‑1 (MQ=255)
                                                  ttAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAAc                                                                                            <  2:273207/149‑1 (MQ=255)
                                                       gcAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCg                                                                                       <  2:740809/149‑1 (MQ=255)
                                                       gcAAATCCCTGATCTTTTAAGCTGCGCTTTAAACCACGGATACAAAAGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCg                                                                                       <  2:259991/149‑1 (MQ=255)
                                                              ccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAgg                                                                                                                 <  2:750017/116‑1 (MQ=255)
                                                              ccTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAgg                                                                                                                 >  1:750017/1‑116 (MQ=255)
                                                               cTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGATACGCCCGGACTACAg                                                                               >  1:986579/1‑149 (MQ=255)
                                                                tGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAg                                                                                                     <  2:729556/126‑1 (MQ=255)
                                                                tGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAg                                                                                                     >  1:729556/1‑126 (MQ=255)
                                                                 gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATCTTACCAGCCGACGGCTACGGTATATGCAACCTGACACATAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGaa                                                                             >  2:885162/1‑149 (MQ=255)
                                                                     ttttAAGCTGCGCTTTTTACCACGGATACAACTGTTACCCGCCGACGTCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGACAAGACAACAGCAATGATACGCCCGGACTACAGAACAtt                                                                         >  1:1051473/1‑149 (MQ=255)
                                                                     ttttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACAAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACagt                                                                         >  1:281097/1‑147 (MQ=255)
                                                                      tttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTa                                                                        >  2:775796/1‑149 (MQ=255)
                                                                       ttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                       <  1:1014122/149‑1 (MQ=255)
                                                                       ttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                       >  1:340946/1‑149 (MQ=255)
                                                                              gcgcTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATAGGCAACCTGACACAAAATTGGGTCATAGTGCAGGAAAAAGCAGTTACCAGGAGCAGACAACAGCAAAGAAACGCCCGGAGTACAGTACATGACAGGcccg                                                                >  2:514467/1‑146 (MQ=255)
                                                                                 gttttTACCACGGATACAAATGTTACCGGCCGACGTGTTCGGTATATGCAACCTGTCACAAAATTGTGTCTTAGTGCAGTAAAAAGCATTTAGCAGGAGCAGACAACAGCAATGAAACGCCCTGACTACAGAACATTACAGGCTCTGGa                                                             <  1:17418/148‑1 (MQ=255)
                                                                                  ttttACCCAGGGATACAAAGGTTACCCGCCGCCGGCTTCGGTATATGCTACCTGACACAAATTTGTGTCTTAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGCCTACAGAACATTACAGGCTCTggtt                                                            >  2:614879/1‑147 (MQ=255)
                                                                                      aaCCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg                                                        <  1:685/149‑1 (MQ=255)
                                                                                       aCCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGt                                                       >  2:304733/1‑149 (MQ=255)
                                                                                          aCGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGAt                                                    <  1:469336/149‑1 (MQ=255)
                                                                                              atagtaaTGTTACCCGACGACGTCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGAATTTACCTTGAGCAGACTACAGCAATGATAGGCCCGGACTACAGAACATTTCAGTCACTGGATGCGGTGATACGt                                                <  1:365906/144‑1 (MQ=255)
                                                                                                  aaaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAAc                                            >  1:584324/1‑149 (MQ=255)
                                                                                                  aaaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAAc                                            >  1:180913/1‑149 (MQ=255)
                                                                                                  aaaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGGGTCTTAGTGCTGGAAAAGGCATTTACCAGGAGCAGACAACAGCAATGAACCGCCCGGCCTACAGAACTTTACAGGCCCTGGATGCGGTGATCCGTGAAc                                            >  2:104156/1‑149 (MQ=255)
                                                                                                        gttcacccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGCCTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGAt                                      >  1:965644/7‑149 (MQ=255)
                                                                                                         aCCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGAtt                                     >  1:960726/1‑149 (MQ=255)
                                                                                                                    gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGAGGCGGTGATACGTGAACGAGGATTTGAGCGCGCgg                          <  2:607946/149‑1 (MQ=255)
                                                                                                                      ttCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTTCAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGca                        >  2:1054044/1‑149 (MQ=255)
                                                                                                                      ttCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGAAGACAACAGCAATGAAAAGCACGGACTACGGAACATTACAGGAAATGGATGCCGTGATACGTGAACGAGGATTTGAGCGCGCGGca                        >  1:863305/1‑149 (MQ=255)
                                                                                                                         ggTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACaa                     >  1:344505/1‑149 (MQ=255)
                                                                                                                              atGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAACCGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCt                >  2:1133487/1‑149 (MQ=255)
                                                                                                                                gCATCCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCtgt              <  2:1125870/149‑1 (MQ=255)
                                                                                                                                gCAACCTGACACAAAATTGTGTCTTAGTGCAGGAAAAATCAATTACCAGGAGCAGACATCAGCTATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGATGATTTGAGCGCGCGGCACAAAAGCtgt              <  1:885162/149‑1 (MQ=255)
                                                                                                                                 cAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCtgtg             <  2:180913/149‑1 (MQ=255)
                                                                                                                                 cAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCtgtg             <  2:748931/149‑1 (MQ=255)
                                                                                                                                  actatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCTGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGc            >  1:1076695/5‑149 (MQ=255)
                                                                                                                                    tatGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCAt          >  1:452816/3‑149 (MQ=255)
                                                                                                                                     cTGACTCAAAATTTTGTCTTAGTTCAGGAAAAAGCTTTTACACGGCGCAGACATCAGCAATGTAACTCCCGGACTTCAGAACATAACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACTAAAGCTGTGCAtt         <  1:816634/149‑1 (MQ=255)
                                                                                                                                     atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCAtt         >  2:673191/2‑149 (MQ=255)
                                                                                                                                        acacAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTaca      <  2:511078/149‑1 (MQ=255)
                                                                                                                                            aaaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATATGAGCGCGCGCCACAAAAGCTGTGCATTACACaaa  >  1:404289/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CGCCACGCCATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTACACAAT  >  NC_000913/3059561‑3059849

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: