Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,059,704 A→T 6.7% 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 6.7% 482.2 / 17.2 135intergenic (‑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 (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 = 9.49e-01

Only 100 of 145 total aligned reads displayed.

ATAACGCCACGCCATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTACACAAT  >  NC_000913/3059557‑3059849
                                                                                                                                                   |                                                                                                                                                 
   acgccacgccATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtg                                                                                                                                                <  3:454408/148‑1 (MQ=255)
    cgccacgccaTTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgt                                                                                                                                               <  4:153906/148‑1 (MQ=255)
               tCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCa                                                                                                                                    <  3:721272/148‑1 (MQ=255)
                   cacaACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAAt                                                                                                                                                  >  4:1112756/1‑130 (MQ=255)
                               gcgAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAcc                                                                                                                    <  1:796873/148‑1 (MQ=255)
                                cgAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTc                                                                                                                                            >  2:498758/1‑123 (MQ=255)
                                  aaCTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAt                                                                                                                         >  4:785929/1‑140 (MQ=255)
                                  aaCTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAt                                                                                                                         <  4:420616/140‑1 (MQ=255)
                                  aaCTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAt                                                                                                                         <  3:785929/140‑1 (MQ=255)
                                  aaCTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCAt                                                                                                                         >  3:420616/1‑140 (MQ=255)
                                    cTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGc                                                                                                                                     <  4:449382/126‑1 (MQ=255)
                                                aaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGc                                                                                                                           <  2:589329/124‑1 (MQ=255)
                                                aaTGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGc                                                                                                                           >  1:589329/1‑124 (MQ=255)
                                                  tGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                                                                   <  1:162951/130‑1 (MQ=255)
                                                      ttAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                              <  2:5695/99‑1 (MQ=255)
                                                       tAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgt                                                                                                                                               >  4:496274/1‑97 (MQ=255)
                                                       tAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                                                                              >  2:135652/1‑130 (MQ=255)
                                                       tAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGc                                                                                                              <  1:135652/130‑1 (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)
                                                                    tGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAg                                                                               >  2:146571/1‑148 (MQ=255)
                                                                     gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTACATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTa                                                                                                                      >  1:341728/1‑108 (MQ=255)
                                                                     gATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTACATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTa                                                                                                                      <  2:341728/108‑1 (MQ=255)
                                                                      aTCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtg                                                                                                                                                <  4:537529/81‑1 (MQ=255)
                                                                      aTCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtg                                                                                                                                                >  3:537529/1‑81 (MQ=255)
                                                                      aTCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                                                                          >  2:416654/1‑87 (MQ=255)
                                                                      aTCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                                                                          <  1:416654/87‑1 (MQ=255)
                                                                      aTCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaaca                                                                                                      <  2:552786/123‑1 (MQ=255)
                                                                      aTCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaaca                                                                                                      >  1:552786/1‑123 (MQ=255)
                                                                       tCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAAc                                                                            <  1:903725/148‑1 (MQ=255)
                                                                        cTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCAt                                                                                                                                          <  4:140400/85‑1 (MQ=255)
                                                                          tttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                              >  3:120542/1‑79 (MQ=255)
                                                                          tttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                              <  4:120542/79‑1 (MQ=255)
                                                                          tttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                  <  3:356123/123‑1 (MQ=255)
                                                                          tttAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                  >  4:356123/1‑123 (MQ=255)
                                                                               gCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCa                                                                                                                   <  2:141634/101‑1 (MQ=255)
                                                                                 tGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGGGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCa                                                                  <  1:326756/148‑1 (MQ=255)
                                                                                     cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                                                                          <  4:690749/72‑2 (MQ=255)
                                                                                     cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAATTTTGTGTCag                                                                                                                                          >  3:690749/1‑71 (MQ=255)
                                                                                     cTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTgg                                                              <  1:146571/148‑1 (MQ=255)
                                                                                             cACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTa                                                                        <  3:136124/130‑1 (MQ=255)
                                                                                             cACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg                                                      >  1:618141/1‑148 (MQ=255)
                                                                                                 gATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAg                                                                                                     >  1:934335/1‑97 (MQ=255)
                                                                                                 gATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAg                                                                                                     <  2:934335/97‑1 (MQ=255)
                                                                                                  aTACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTAcc                                                                                                                    <  1:408045/81‑1 (MQ=255)
                                                                                                  aTACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                  <  1:248231/99‑1 (MQ=255)
                                                                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACt                                                                >  4:642519/1‑130 (MQ=255)
                                                                                                     cAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACt                                                                <  3:642519/130‑1 (MQ=255)
                                                                                                       aaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTc                                                                                                                                            >  3:71841/1‑52 (MQ=255)
                                                                                                       aaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTc                                                                                                                                            <  4:71841/52‑1 (MQ=255)
                                                                                                       aaTGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                       <  4:17675/121‑1 (MQ=255)
                                                                                                         tGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                              >  1:531696/1‑48 (MQ=255)
                                                                                                         tGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATtgtg                                                                                                                                              >  3:362371/1‑48 (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)
                                                                                                              cccGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGa                                                     <  2:1070426/132‑1 (MQ=255)
                                                                                                               ccgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                                                                            <  4:151198/60‑1 (MQ=255)
                                                                                                               ccgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAg                                                                                                                            >  3:151198/1‑60 (MQ=255)
                                                                                                                cgccgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAg                                                                     >  3:157469/1‑114 (MQ=255)
                                                                                                                   cgACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGa                                                                              <  3:1003220/102‑1 (MQ=255)
                                                                                                                    gACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGaca                                                                                                         <  2:684227/74‑1 (MQ=255)
                                                                                                                    gACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGaca                                                                                                         >  1:684227/1‑74 (MQ=255)
                                                                                                                      cGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCg                                                         <  2:302860/120‑1 (MQ=255)
                                                                                                                      cGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCg                                                         >  1:302860/1‑120 (MQ=255)
                                                                                                                       ggCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg                                                        <  3:398283/120‑1 (MQ=255)
                                                                                                                       ggCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCgg                                                        >  4:398283/1‑120 (MQ=255)
                                                                                                                        gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAg                                         >  4:743213/1‑134 (MQ=255)
                                                                                                                        gCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAg                                         <  3:743213/134‑1 (MQ=255)
                                                                                                                           tCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGca                        >  3:1038232/1‑148 (MQ=255)
                                                                                                                               tataTGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGcc                                                                                         >  4:916318/1‑79 (MQ=255)
                                                                                                                               tataTGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGcc                                                                                         <  3:916318/79‑1 (MQ=255)
                                                                                                                                 tatGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCaa                                                                                                  <  2:255424/68‑1 (MQ=255)
                                                                                                                                  atGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGAtt                                     <  4:372877/128‑1 (MQ=255)
                                                                                                                                    gCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                       <  4:1029659/92‑1 (MQ=255)
                                                                                                                                    gCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                       >  3:1029659/1‑92 (MQ=255)
                                                                                                                                       aCCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGTACATTAc                                                                       <  2:221267/89‑1 (MQ=255)
                                                                                                                                         atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGa                                                                              <  3:693133/79‑1 (MQ=255)
                                                                                                                                         atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGa                                                                              >  4:693133/2‑80 (MQ=255)
                                                                                                                                         atGACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCAt          >  4:955939/2‑148 (MQ=255)
                                                                                                                                           gACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATg                                                           >  2:693306/1‑97 (MQ=255)
                                                                                                                                           gACACAATTTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATg                                                           <  1:693306/97‑1 (MQ=255)
                                                                                                                                           gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTg                                                      <  1:103907/102‑1 (MQ=255)
                                                                                                                                           gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGaa                                             <  2:862143/111‑1 (MQ=255)
                                                                                                                                           gACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGaa                                             >  1:862143/1‑111 (MQ=255)
                                                                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGacaac                                                                                                       <  3:407524/49‑1 (MQ=255)
                                                                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                       <  2:946632/81‑1 (MQ=255)
                                                                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTAc                                                                       >  1:946632/1‑81 (MQ=255)
                                                                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCAc                                                                 >  2:740980/1‑87 (MQ=255)
                                                                                                                                               caAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCAc                                                                 <  1:740980/87‑1 (MQ=255)
                                                                                                                                                aaaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGaatc                                                                                            <  2:240202/59‑3 (MQ=255)
                                                                                                                                                aaaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAAc                                                                                            >  1:240202/1‑59 (MQ=255)
                                                                                                                                                 aaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACaaa                    <  4:1010896/130‑1 (MQ=255)
                                                                                                                                                 aaaTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTACACaac  >  1:134016/1‑147 (MQ=255)
                                                                                                                                                   aTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGa                                          <  2:107427/106‑1 (MQ=255)
                                                                                                                                                   aTTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGa                                          >  1:107427/1‑106 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATAACGCCACGCCATTCAACACAACCAGGGCGCGAACTCGCTGAGCGAAATGCATTAGCGCAAATCCCTGATCTTTTAAGCTGCGCTTTTAACCACGGATACAAATGTTACCCGCCGACGGCTTCGGTATATGCAACCTGACACAAAATTGTGTCATAGTGCAGGAAAAAGCATTTACCAGGAGCAGACAACAGCAATGAAACGCCCGGACTACAGAACATTACAGGCACTGGATGCGGTGATACGTGAACGAGGATTTGAGCGCGCGGCACAAAAGCTGTGCATTACACAAT  >  NC_000913/3059557‑3059849

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