Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,639,440 C→T T421T (ACC→ACT)  creD → putative inner membrane protein CreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,639,4400CT97.6% 111.7 / NA 42T421T (ACC→ACT) creDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/24);  minor base G (0/1);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCC  >  NC_000913/4639210‑4639659
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                           
tagtcaatcGGTTTTTCCGTGGTATGGTTGATCGCCATGCTGATTGGGGCGTTAACTATCGGTTTTTATTTGCAGTCGGTATTGAAACGCTG‑TGCCAACAGCATGTTGGTTTCCCTCGCTGTGTTGTTGCTGGATGGTGTTATGTGGTGACTGCTCATCTCGGCCGATAGCGTGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGTCATGATGTTTGTGACTCGGAATATCGACTGGTATg                                                                                                                                                                                                          <  2:203731/242‑1 (MQ=255)
 aaCATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGGTGGTGTGATGTGGGGGCTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGGCTGGTATGCg                                                                                                                                                                                                        >  1:127403/1‑251 (MQ=255)
          gTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTg                                                                                                                                                                                               <  1:474330/251‑1 (MQ=255)
               aCCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGGTGGGGTGATGTGGGGGCTGCTCAACTCTGCCGCTAGCGCGCTGTTGTTTGGGACCGGTGTGCTGGTGGTGGGGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCggg                                                                                                                                                                                          >  1:60846/1‑249 (MQ=255)
               aCCGTGGCATGGATAATCGCCAGTCGGATTGGTGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCg                                                                                                                                                                                            <  2:60846/249‑1 (MQ=255)
                  gTGGCATGGATAATCGGCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCATGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAAt                                                                                                                                                                                       <  2:274055/251‑1 (MQ=255)
                       attgATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGAGTGGTGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGGTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGGATGCGGTTTCACTGCCGAAAATGaaa                                                                                                                                                                                   >  1:641644/4‑250 (MQ=255)
                       aTGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTTCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGGACCAGTGGGCTGGAGGTGGCGCTGGCCGGCCTGATGTTTGTGACTCGTAATATCGGCTGGGATGCGTTTTTACTGCCGAAAATGaaa                                                                                                                                                                                   >  1:87648/1‑250 (MQ=255)
                         ggATAATCTCAAGTCTGATTGGGGCGTTAATGAACGGTATGTATTTGCAGGCGGTATTGAAAGGCTGGCGT‑AATAGCATGTTGTTTTCCCTCGCGCTGTTGTTGCTGGAGGGTGTGATGTGGGGACTGCTGAAATCTTCCGATAGCGCTCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCGGGCCGGCATGATGTTTGTGACTCGTAATAGCGACTGGTATGCGTTTTCACTGCTTAAAATGAAAGc                                                                                                                                                                                 <  2:509598/250‑1 (MQ=255)
                             aaTCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTGCACTGCCGAAAATGAAAGCCAGt                                                                                                                                                                             <  1:246533/250‑1 (MQ=255)
                                cGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTaaa                                                                                                                                                                          <  1:587575/250‑1 (MQ=255)
                                              gggCGTTAATGAACGGTAGTTATTTGCAGGCGGTAGTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGGGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGAc                                                                                                                                                           <  1:347379/251‑1 (MQ=255)
                                                    tAATTAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGAtca                                                                                                                                                      <  1:495889/250‑3 (MQ=255)
                                                       ggAACGGTATTTATTTGCTGGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGTTGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGAtcaatta                                                                                                                                                  <  2:344289/250‑7 (MQ=255)
                                                               atttatttGCAGGCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCt                                                                                                                                           <  1:236991/250‑1 (MQ=255)
                                                                      ttcaggcgttATTTAAAGGCTGGCGC‑AACCGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGTTGGTGTGATGTGGGGACTTCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTTGTGCTGGCCGGCAAGATGTTTTTGACTCGTAATATCGACTGGTATGCGTTTTCTCTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGAGCAATTACGTATCTGGAAATa                                                                                                                                    <  2:528808/241‑1 (MQ=255)
                                                                         aggCGGTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGt                                                                                                                                <  1:291061/249‑1 (MQ=255)
                                                                          ggctgTATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCGGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGAGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGtt                                                                                                                               <  2:87648/247‑1 (MQ=255)
                                                                               tATTGAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAAAGGCCGATCAATTACGTATCTGGAAATAAGGTTGaaa                                                                                                                           >  2:507914/1‑250 (MQ=255)
                                                                                   gAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCTAAAATGAAAGCCAGTAAAGAAGTTAAAATGGACGATCAATTACGTATCTGTAAATAAGGTTGAAAAATa                                                                                                                       >  2:356523/1‑250 (MQ=255)
                                                                                   gAAAGGCTGGCGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCGGTATGGTGTGATGTGGGGACTTCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGGTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACTTATCTGGACATAAGGTTGAAAAATaa                                                                                                                      <  2:86822/251‑1 (MQ=255)
                                                                                             cgC‑AACAGCATGTTGTTTTCCCTCTCGCTGTTGTTGCGGGATGGGGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCt                                                                                                            <  2:629753/250‑1 (MQ=255)
                                                                                                                 cccTCGCGCTGTTGTTGCTGGATGGGGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTCTCTGGAAATAAGGTTGAAAAATAACACCGGCGCTAATAAGCGCCGtttttttt                                                                                         >  2:120998/1‑251 (MQ=255)
                                                                                                                         cTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACggtgg                                                                                 >  1:107832/1‑251 (MQ=255)
                                                                                                                         cTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACggtgg                                                                                 >  1:144326/1‑251 (MQ=255)
                                                                                                                          tgttgtTGCTGGATTGAGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCCCTAAAAAGCGCCGTTTTTTTTGACggtgg                                                                                 >  2:484887/1‑250 (MQ=255)
                                                                                                                              tttgcaggatggaggggtgtgGGGACTGCTCAACTCTGCCGCTAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGGGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGAATGCGGTTTCAATGACGAAAATGAAAGCCAGTAAAGAAGTTACAAGGGCCGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGACGATAAAAAGCGCCGTTTTTTTTGAAGGTGGTaaa                                                                             >  2:633477/18‑250 (MQ=255)
                                                                                                                                 gCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCg                                                                         >  1:167574/1‑251 (MQ=255)
                                                                                                                                   tGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGAtcaa                                                                                                                                                     >  1:245667/1‑170 (MQ=255)
                                                                                                                                   tGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGAtcaa                                                                                                                                                     <  2:245667/173‑4 (MQ=255)
                                                                                                                                      atggtgtgtcgtgtgtacgtcTCAACCCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTTCTGGTGGTGGCGCTGGACGGCATGATGTTTGTGACTCGTAATATCGACTGTTATGTGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAAGTACGGATCTGGAAATAAGGTTGAAAAATAAAAACTTCGCTAAAAAGCGCCGTTTTTTTTGACGGTTGTAAAGCTGATTa                                                                     <  2:598895/230‑1 (MQ=255)
                                                                                                                                       tGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTTGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTa                                                                                                                                                                            <  2:479656/146‑1 (MQ=255)
                                                                                                                                       tGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTa                                                                                                                                                                            >  1:479656/1‑146 (MQ=255)
                                                                                                                                          tgtgATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCtt                                                                >  1:411397/1‑251 (MQ=255)
                                                                                                                                                           ccaattctGCCGAGAGCGCGGTGTTGTGGGGAACCAGTGTGCTGGTGGTGGCGCTTGCCGGCATGATGTTTGTGAATCGTAATATCGACTGGTAAGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGCCGATCAATTTCGTAGCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGGTTTTGTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCTCCACAGAAg                                               <  2:428324/246‑1 (MQ=255)
                                                                                                                                                                   gCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAAcc                                       <  1:41578/251‑1 (MQ=255)
                                                                                                                                                                          gcgcgcTTTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCg                                <  1:474040/251‑1 (MQ=255)
                                                                                                                                                                                           actagtgttgtgggggtggtGCTGGCCGGCAGGATGTTTGTGCAGCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCtatgat               <  2:171693/237‑1 (MQ=255)
                                                                                                                                                                                              attgtgCTGGGGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCTAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACTGCGCTAAAAAgc                                                                                                     <  2:494020/160‑1 (MQ=255)
                                                                                                                                                                                              aGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAACTAAGGTTGAAAAATAAAAACGGCGCTAAAAAgc                                                                                                     >  1:494020/1‑162 (MQ=255)
                                                                                                                                                                                                 gtgCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTtccgg         <  1:42172/251‑1 (MQ=255)
                                                                                                                                                                                                        tggtggCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTtccggcgtatcc  >  1:15708/1‑251 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                           
GAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGC‑AACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCC  >  NC_000913/4639210‑4639659

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