Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,400 A→G N41N (AAT→AAC)  ettA ← energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,4000AG94.1% 88.6 / ‑3.1 37N41N (AAT→AAC) ettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (15/17);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

CCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCC  >  NC_000913/4630159‑4630640
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
ccaccTTGTCAAACCCGGCACCCGGACCGGCGTACGGCGCATCCCCTCCACCCGGGCGTTCCAGGGGTTTACCCACTCCAGAACCCGCTTCTTCAAGGGCCTCACGCAGGGTGTGTTCCGGGTTCAGCTGCGTTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGCCGCGCTTCCCCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGAccc                                                                                                                                                                                                                                          <  1:531106/247‑1 (MQ=255)
       gTCAAAATCGGCCTCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCATTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAAACAATCTTGATGTCTGGCTGCGGCCGCGCCTCACCTTCGATGTCTTTCTCAATCCCC‑CCCATAATGCGACGCAGGGTAGAATTGCCCGCGCCGTTCAGAcaccgggccc                                                                                                                                                                                                                                   >  2:202344/1‑242 (MQ=255)
                       ggtgcGGCGTACAGCGCATACCCTTCACCCAGGCTTTTCAGGGCTTTAACCACTTCAAAACCCGTTTCTTCAATGGACTCCCGCACGGTGTGTTCCGGGTTCAGCGGCGTTTCCTGCGGCAGTTAACCAATCTTTATGTCTGGCGGCGGCCCCGTTTCCCCTTCGATTTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAgg                                                                                                                                                                                                                  <  1:202344/247‑1 (MQ=255)
                          ccGGCGTACAGCGCAACCCCTTCATCCAGGCGTTCCAGGGCGTTACCCCCTCCAAAACCCTCTTCTTCAATGGACCCCCGCAGGGTGTGTTCCGGGTTCAGCTGCGTTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAAGGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGgaa                                                                                                                                                                                                               <  1:618337/250‑1 (MQ=255)
                               gTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGAAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAAc                                                                                                                                                                                                          >  1:55247/1‑251 (MQ=255)
                                tACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCATTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGCAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGACGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTCGACATTCCCCCACCCATCAAACCCAGGCCCCC‑AAATTTTACCCCAggggaaaaacc                                                                                                                                                                                                         >  2:363447/1‑243 (MQ=255)
                                  cAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAACCAATCTTGATGTCTGGCTGCGGAAGCGCTTCACCTTCGATGTCTTTATTAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCCCCGTTCAGACCCAGGACACC‑CATTTTTTCCCCCGGGGAGAAACTcc                                                                                                                                                                                                       >  2:329570/1‑250 (MQ=255)
                                        aTACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGGCTTTTTCAATGCCC‑GCCATAATGCGCACCAGGGTAGGCTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCAAGGGAAGAAACTCcgag                                                                                                                                                                                                    >  1:275853/1‑248 (MQ=255)
                                        aTACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAAACAATCTTGATGTCTGGCTGCGGACGCGCCTCACCTTCGGTGGCTTTATCAATGCCC‑GCCATAATGCGCCGCAGGGTAGACTTACCCCCGCCGTTCAGGCCCAGGACACC‑AATTTTTGTCCCAGGGAAGaaaatccgagggg                                                                                                                                                                                                 >  2:342870/1‑241 (MQ=255)
                                             cTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt                                                                                                                                                                                            >  1:141965/1‑251 (MQ=255)
                                               ccATCCAGGCGTTTAAGGGCGTTAACCCCTACAGAACCCGCTTCTTCATGGGCCTCACGCGCGGTGTGTTCGGGTTTCGGCTGCGGTTCCTGCGGCAGATAACCATTCTTGAGGTCTGGCTGCGGACGCGCTTCACCTTGGATGTTTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGtttt                                                                                                                                                                                           <  1:328420/249‑1 (MQ=255)
                                               ccATCCAGGCGTTACAGGGCGTAAACCGCTTAAGAAACCTCTTCTTCAGGGGCCTCACGCACGGTGTGTCGCGGGTTCAGCTGGGGTTCCTGGGGCAGGTAACCAATCTGTATGTCGGGCTGCGGCCGCGCTTCACCTTCGTTGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt                                                                                                                                                                                            <  1:367624/248‑1 (MQ=255)
                                                           ttcggggCGTAAACCACTTAAGAACCCGTTTCTTCAATGGACTCCCGCCGGGTTTGTTCCGGTTTCAGCTGCGTTCCCGGCGGCAGGTAACAAATTTTGGTGTCGGGCTGCGGACGCGCTCCCCCTCCGATGTCTTTATCAAGGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGAc                                                                                                                                                                               <  1:631452/246‑1 (MQ=255)
                                                                             taaaaaaCCGCTTCTTCAATGGACTCCCCCCGGGTGTGTCCCGGGTTCAGCGGCGGTTCCTGCGGCAGTTACCCAATTTTGATGTCTGGCGGCGGCCGCCCTTCCCCTTCGATGTCTTAATCAAGGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACtt                                                                                                                                                             <  1:613994/246‑1 (MQ=255)
                                                                                  aaCCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCg                                                                                                                                                        >  1:598542/1‑250 (MQ=255)
                                                                                                     ccccgcgcggggtgtTCGGGGTTCCGCTGGGGTTCCGGCGGCAGGTACCCATTCTTGATGTCTGGCTGCGGACGCGCTTCCCCTGCGATGTCTTTA‑CAAAGCCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCACGGTATAAAc                                                                                                                                    <  1:605224/240‑1 (MQ=255)
                                                                                                                ggtttCCGGTTTCATCTGCGTTCCCTGCGGCAGGTACCCAATCTTGTTGTCTGGCGGCGGCCCCTCTTCCCCTTCAATGTCTTTACCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCa                                                                                                                         <  1:656151/248‑1 (MQ=255)
                                                                                                                        gtttCAGCGGCGGTCCCGGGGGCAGTTACCCAATCTTGATGTCGGGCTGCGGCCGCTCTTCCCCTTCGATGTTTTTATAAAGCCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGAc                                                                                                                 <  1:511637/249‑1 (MQ=255)
                                                                                                                          ttCAGCTGCGGTTTCTGCGGCAGGTAACCAATCTTGCTGTCTGGCTGAGGTCGCGCCTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCCCGCCCTTCAGACCCAAGACACC‑AATTTTTGCCCCAGGGCAGAAACTCAAAGAGATGTTTTTCAAAATATGACCTTTCGGCGGAACAACTTTGCCGACACGGTGAGTGATATTAACGCATTACGCCACGTTGGACtt                                                                                                               >  2:360874/1‑251 (MQ=255)
                                                                                                                                             gCAGGTAACCAATCTTGATGTATGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGGGCCACGTTGGACTTCGCCTCTTTGTTTATCGTg                                                                                            >  2:85201/1‑251 (MQ=255)
                                                                                                                                                 aTAACCAATCTTGTTGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATaa                                                                                        <  1:113674/251‑1 (MQ=255)
                                                                                                                                                          cTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTc                                                                               >  1:416388/1‑251 (MQ=255)
                                                                                                                                                          cTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACCTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTc                                                                               >  1:192713/1‑251 (MQ=255)
                                                                                                                                                            tGATGTCTGGCGGCGGACGCGCTTCCCCTTCGAGGTTTTTACAAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTc                                                                               <  1:85201/249‑1 (MQ=255)
                                                                                                                                                                   gggtcgggagcccgcTTCCCCTTGGATTTCTTTACCAAGCCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGaa                                                                      <  1:363447/239‑1 (MQ=255)
                                                                                                                                                                                    cccctgcGATGTCTTTCTCAAGGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGTAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTcccc                                                     <  1:179906/245‑1 (MQ=255)
                                                                                                                                                                                           gATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGAAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGGCACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGCTGATGAGTTTTCAAAGGCGAAGTGTAGCATTTTTCCCCGCCTAAt                                              >  2:496349/1‑251 (MQ=255)
                                                                                                                                                                                             tttCTTTATCAACCCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCTTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATc                                             <  1:30381/248‑1 (MQ=255)
                                                                                                                                                                                             gttgtaaactaagcccc‑GCCATAATGTGCAGCAGGTGGGTCGACCCCGCGCCGTTCAGACCCAGGACCCC‑AATTTGTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGAGGCTTGGTATAAACAAATTGGGCCCCGTTGGACTTCGCCTCTATTTTTATCGTGATGATGAGTTGTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATc                                             <  2:463205/236‑1 (MQ=255)
                                                                                                                                                                                                 tttATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAgt                                        >  1:640634/1‑250 (MQ=255)
                                                                                                                                                                                                            cc‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCAc                             >  1:60789/1‑251 (MQ=255)
                                                                                                                                                                                                               gCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGCCCCATGACACC‑AATTTTTGCCCCAGTGACGAAACTCAGAGAGATGTTTTTCAACATATGACGTTTCGGCGGAACAACTTTGCCGACCCGATGCATGGTCTAAACGATTTGAGCCTCGTTGGCCTTCGCCTCAATTTTTATCGTGATGCTGAGTTTTCAAAGGCGAAGTGTCTCCTTTTTTCCCGCCTAATCCCAGTTGGTCGAT‑ACCCt                          >  2:294110/1‑251 (MQ=255)
                                                                                                                                                                                                                       gcaccgtggggtagaattACCCGAGCCATTCAGACCCATGA‑ACCTAATTTTTGCCCCAGGGAAGAATCTCACAGATATGTTTTTCAAAATATGTCGGTTCGGCGGAACAATTTTTCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGCCTTCGCCTCTATGTTTATATTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCTTAATCCCAGTTGGTCGATCACCTTCGCagta                   <  2:604898/235‑1 (MQ=255)
                                                                                                                                                                                                                         gcagcaGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCagtaaa                 <  1:329570/251‑1 (MQ=255)
                                                                                                                                                                                                                         gcagcaGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCagtaa                  >  1:589543/1‑250 (MQ=255)
                                                                                                                                                                                                                                     aCTTACCCGCGCCGTTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGtt                                                                                                                     <  1:413013/139‑1 (MQ=255)
                                                                                                                                                                                                                                         aCCCGCGCCGTTCAGACCCAGGAAACC‑AATTTTTGCCCCAGGGAAGAAACTTCGAGAGATGTTTTTAAAAATATGACTGTTCGGCGGAACAACCTGGCCGACCCGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGAGGAGTTTGCAAAGGCGAATTGTAGCATTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTcc  <  2:192713/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
CCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACC‑AATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCC  >  NC_000913/4630159‑4630640

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