Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,059,703 A→T 6.7% 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.7% 480.1 / 17.1 135intergenic (‑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 (60/66);  new base T (5/4);  total (65/70)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 141 total aligned reads displayed.

ATAACGCCACGCCATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTACACAA  >  NC_000913/3059557‑3059848
                                                                                                                                                  |                                                                                                                                                 
aTAACGCCACGCCATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACaaaa                                                                                                                                                  >  2:828021/1‑148 (MQ=255)
   acgccacgccATTCAACCCATCCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAACTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtg                                                                                                                                               <  1:828021/148‑1 (MQ=255)
   acgccacgccATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtg                                                                                                                                               <  3:454408/148‑1 (MQ=255)
    cgccacgccaTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgt                                                                                                                                              <  4:153906/148‑1 (MQ=255)
               tCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCa                                                                                                                                   <  3:721272/148‑1 (MQ=255)
                   cacaACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAAt                                                                                                                                                 <  3:1112756/130‑1 (MQ=255)
                   cacaACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAAt                                                                                                                                                 >  4:1112756/1‑130 (MQ=255)
                                cgAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTc                                                                                                                                           >  2:498758/1‑123 (MQ=255)
                                cgAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTc                                                                                                                                           <  1:498758/123‑1 (MQ=255)
                                  aaCTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAt                                                                                                                        <  4:420616/140‑1 (MQ=255)
                                  aaCTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAt                                                                                                                        >  4:785929/1‑140 (MQ=255)
                                  aaCTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAt                                                                                                                        <  3:785929/140‑1 (MQ=255)
                                    cTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGc                                                                                                                                    >  3:449382/1‑126 (MQ=255)
                                                aaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGc                                                                                                                          <  2:589329/124‑1 (MQ=255)
                                                aaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGc                                                                                                                          >  1:589329/1‑124 (MQ=255)
                                                  tGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                                                                  >  2:162951/1‑130 (MQ=255)
                                                  tGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                                                                  <  1:162951/130‑1 (MQ=255)
                                                      ttAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                             >  1:5695/1‑99 (MQ=255)
                                                      ttAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                             <  2:5695/99‑1 (MQ=255)
                                                       tAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgt                                                                                                                                              >  4:496274/1‑97 (MQ=255)
                                                       tAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                                                                             <  1:135652/130‑1 (MQ=255)
                                                       tAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                                                                             >  2:135652/1‑130 (MQ=255)
                                                             aaaTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                                                                             <  2:697322/140‑1 (MQ=255)
                                                             aaaTCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaa                                                                                             >  1:697322/1‑140 (MQ=255)
                                                                   cTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACa                                                                               >  2:903725/1‑148 (MQ=255)
                                                                    tGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCgg                                                                                     <  3:470392/141‑1 (MQ=255)
                                                                    tGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCgg                                                                                     >  4:470392/1‑141 (MQ=255)
                                                                     gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTACATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTa                                                                                                                     >  1:341728/1‑108 (MQ=255)
                                                                     gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTACATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTa                                                                                                                     <  2:341728/108‑1 (MQ=255)
                                                                      aTCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtg                                                                                                                                               >  3:537529/1‑81 (MQ=255)
                                                                      aTCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtg                                                                                                                                               <  4:537529/81‑1 (MQ=255)
                                                                      aTCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                                                                         >  2:416654/1‑87 (MQ=255)
                                                                      aTCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGaaaa                                                                                                                             >  3:805660/1‑99 (MQ=255)
                                                                       tCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAAc                                                                           <  1:903725/148‑1 (MQ=255)
                                                                        cTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                                                                         >  3:140400/1‑85 (MQ=255)
                                                                        cTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                                                                         <  4:140400/85‑1 (MQ=255)
                                                                          tttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                             >  3:120542/1‑79 (MQ=255)
                                                                          tttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                 >  4:356123/1‑123 (MQ=255)
                                                                          tttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                 <  3:356123/123‑1 (MQ=255)
                                                                              aGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                                                                         <  4:696594/79‑1 (MQ=255)
                                                                              aGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATGGTGTCAt                                                                                                                                         >  3:696594/1‑79 (MQ=255)
                                                                               gCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                                                                  >  1:141634/1‑101 (MQ=255)
                                                                               gCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                                                                  <  2:141634/101‑1 (MQ=255)
                                                                                  gcgcTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                 >  4:677686/1‑115 (MQ=255)
                                                                                     cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                                                                         >  3:690749/1‑71 (MQ=255)
                                                                                     cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                                                                         <  4:690749/72‑2 (MQ=255)
                                                                                             cACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTa                                                                       >  4:136124/1‑130 (MQ=255)
                                                                                             cACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg                                                     >  1:618141/1‑148 (MQ=255)
                                                                                                 gATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAg                                                                                                    <  2:934335/97‑1 (MQ=255)
                                                                                                  aTACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAcc                                                                                                                   <  1:408045/81‑1 (MQ=255)
                                                                                                  aTACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                 >  2:248231/1‑99 (MQ=255)
                                                                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACt                                                               <  3:642519/130‑1 (MQ=255)
                                                                                                       aaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTc                                                                                                                                           <  4:71841/52‑1 (MQ=255)
                                                                                                         tGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                             <  4:362371/48‑1 (MQ=255)
                                                                                                         tGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                             >  1:531696/1‑48 (MQ=255)
                                                                                                         tGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                             <  2:531696/48‑1 (MQ=255)
                                                                                                         tGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                             >  3:362371/1‑48 (MQ=255)
                                                                                                          gTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACAt                                                                         <  1:404811/115‑1 (MQ=255)
                                                                                                          gTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACAt                                                                         >  2:404811/1‑115 (MQ=255)
                                                                                                          gTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                      <  3:995649/118‑1 (MQ=255)
                                                                                                          gTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGTCACTGGATGCGGTGATACGTGAACGa                                         <  3:759158/147‑1 (MQ=255)
                                                                                                          gTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGa                                         >  4:759158/1‑147 (MQ=255)
                                                                                                              cccGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGa                                                    >  1:1070426/1‑132 (MQ=255)
                                                                                                               ccgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                                                                           >  3:151198/1‑60 (MQ=255)
                                                                                                                cgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAg                                                                    <  4:157469/114‑1 (MQ=255)
                                                                                                                cgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAg                                                                    >  3:157469/1‑114 (MQ=255)
                                                                                                                   cgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGa                                                                             >  4:1003220/1‑102 (MQ=255)
                                                                                                                   cgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGa                                                                             <  3:1003220/102‑1 (MQ=255)
                                                                                                                    gACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGaca                                                                                                        >  1:684227/1‑74 (MQ=255)
                                                                                                                    gACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGaca                                                                                                        <  2:684227/74‑1 (MQ=255)
                                                                                                                      cGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCg                                                        >  1:302860/1‑120 (MQ=255)
                                                                                                                      cGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCg                                                        <  2:302860/120‑1 (MQ=255)
                                                                                                                       ggCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg                                                       <  3:398283/120‑1 (MQ=255)
                                                                                                                        gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAg                                        <  3:743213/134‑1 (MQ=255)
                                                                                                                        gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAg                                        >  4:743213/1‑134 (MQ=255)
                                                                                                                           tCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACg                                                <  1:1124279/123‑1 (MQ=255)
                                                                                                                           tCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACg                                                >  2:1124279/1‑123 (MQ=255)
                                                                                                                           tCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGca                       <  4:1038232/148‑1 (MQ=255)
                                                                                                                               tataTGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGcc                                                                                        >  4:916318/1‑79 (MQ=255)
                                                                                                                                 tatGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                 >  1:255424/1‑68 (MQ=255)
                                                                                                                                 tatGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                 <  2:255424/68‑1 (MQ=255)
                                                                                                                                  atGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGAtt                                    <  4:372877/128‑1 (MQ=255)
                                                                                                                                    gCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                      >  3:1029659/1‑92 (MQ=255)
                                                                                                                                    gCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                      <  4:1029659/92‑1 (MQ=255)
                                                                                                                                       aCCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGTACATTAc                                                                      <  2:221267/89‑1 (MQ=255)
                                                                                                                                       aCCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                      >  1:221267/1‑89 (MQ=255)
                                                                                                                                         atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGa                                                                             <  3:693133/79‑1 (MQ=255)
                                                                                                                                         atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCAt         >  4:955939/2‑148 (MQ=255)
                                                                                                                                           gACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATg                                                          >  2:693306/1‑97 (MQ=255)
                                                                                                                                           gACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATg                                                          <  1:693306/97‑1 (MQ=255)
                                                                                                                                           gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg                                                     >  2:103907/1‑102 (MQ=255)
                                                                                                                                           gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg                                                     <  1:103907/102‑1 (MQ=255)
                                                                                                                                           gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGaa                                            >  1:862143/1‑111 (MQ=255)
                                                                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaac                                                                                                      >  4:407524/1‑49 (MQ=255)
                                                                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaac                                                                                                      <  3:407524/49‑1 (MQ=255)
                                                                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                      >  1:946632/1‑81 (MQ=255)
                                                                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCAc                                                                >  2:740980/1‑87 (MQ=255)
                                                                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCAc                                                                <  1:740980/87‑1 (MQ=255)
                                                                                                                                                aaaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAAc                                                                                           >  1:240202/1‑59 (MQ=255)
                                                                                                                                                 aaaTTGTGTCATAGTGCAGGAAAAAGGATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACaaa                   >  3:1010896/1‑130 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATAACGCCACGCCATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTACACAA  >  NC_000913/3059557‑3059848

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