Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,531,391 G→A intergenic (‑158/‑72) zipA ← / → cysZ cell division protein ZipA/sulfate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,531,3910GA98.5% 224.5 / ‑3.7 67intergenic (‑158/‑72)zipA/cysZcell division protein ZipA/sulfate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/37);  minor base . (0/1);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTG  >  NC_000913/2531256‑2531521
                                                                                                                                       |                                                                                                                                  
ccTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCa                                                                                                                                <  1:779477/140‑1 (MQ=255)
  ttGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTTCTTGTTGTTGTTTTTTAGAACCAATCCTAACTATCCAtt                                                                                                                              >  2:380255/1‑140 (MQ=255)
   tGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTg                                                                                                                             >  2:639546/1‑140 (MQ=255)
    gCCACCACGGCAAACATTTACTCACTAAGAGCATTTGCCGATTACCTCAAGTGCAAGCGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTgc                                                                                                                            >  2:175981/1‑140 (MQ=255)
        ccacGGCGAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCaa                                                                                                                        <  2:169954/140‑1 (MQ=255)
        ccacGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCaa                                                                                                                        <  1:41912/140‑1 (MQ=255)
             gCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTAcc                                                                                                                   <  2:591489/140‑1 (MQ=255)
               aaaCATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCg                                                                                                                 <  2:497203/140‑1 (MQ=255)
               aaaCATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCg                                                                                                                 <  1:474112/140‑1 (MQ=255)
                    tttACTCACTAAGAGGATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTaaa                                                                                                            <  1:292899/140‑1 (MQ=255)
                      tACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAAt                                                                                                          <  1:653170/140‑1 (MQ=255)
                       aCTCAATAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATa                                                                                                         >  2:507900/1‑140 (MQ=255)
                             tAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTGCCCGGTAAAATACGCAGa                                                                                                   <  1:507900/140‑1 (MQ=255)
                             tAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGa                                                                                                   >  2:583/1‑140 (MQ=255)
                             tAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGa                                                                                                   >  1:539782/1‑140 (MQ=255)
                               agagTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAAt                                                                                                 <  1:64885/140‑1 (MQ=255)
                                     tttGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATAGAGAAATTACCCGGTAAATTACGCAGTATTTTCCt                                                                                           >  2:365803/1‑140 (MQ=255)
                                       tGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTgg                                                                                         >  1:141441/1‑140 (MQ=255)
                                       tGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTgg                                                                                         >  2:708731/1‑140 (MQ=255)
                                         ccGATTACCTCAAGTGCGAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGa                                                                                       >  2:35204/1‑140 (MQ=255)
                                         ccGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGa                                                                                       <  1:547063/140‑1 (MQ=255)
                                         ccGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATAGGGCAATTAACCGGTAGAATACGAAGAATTTTCCTGGGa                                                                                       >  1:241266/1‑140 (MQ=255)
                                          cGATTACCTCAAGTGTAAGTGCGCTATTATCTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGAt                                                                                      <  1:380255/140‑1 (MQ=255)
                                           gATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGAtt                                                                                     >  2:193449/1‑140 (MQ=255)
                                           gATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGAtt                                                                                     >  2:468184/1‑140 (MQ=255)
                                           gATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGAtt                                                                                     >  1:598825/1‑140 (MQ=255)
                                                 cTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCTGAATTTTCCTGGGATTGGt                                                                                  >  2:293003/1‑137 (MQ=255)
                                                 cTCAAGTGCAAGTGCACTATTAACTTTCACAGCAAAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGt                                                                                  <  1:293003/137‑1 (MQ=255)
                                                  tCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCaaa                                                                              <  1:481999/140‑1 (MQ=255)
                                                       tgcaagtgcaCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAg                                                                         <  1:639546/140‑1 (MQ=255)
                                                       tgcaagtgcaCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAg                                                                         <  1:193449/140‑1 (MQ=255)
                                                       tgcaagtgcaCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAg                                                                         <  2:638430/140‑1 (MQ=255)
                                                       tgcaagtgcaCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAg                                                                         <  1:525223/140‑1 (MQ=255)
                                                         caagtgcaCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCt                                                                       >  2:426320/1‑140 (MQ=255)
                                                           agtgcaCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCt                                                                       >  2:603331/1‑138 (MQ=255)
                                                           agtgcaCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCt                                                                       <  1:603331/138‑1 (MQ=255)
                                                             tgcaCTATTAACTTTAACAGCAAAATGATCGATGAAAGCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATc                                                                   <  1:363004/140‑1 (MQ=255)
                                                              gcaCTATTAACTTTCGCAGCACAAAGATAAATGAAATAGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATcc                                                                  <  1:399653/140‑1 (MQ=255)
                                                              gcaCTATTAACTTTCACAGCACAAATATATATGAAATGGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATcc                                                                  <  1:365803/140‑1 (MQ=255)
                                                               caCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCCGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCt                                                                 >  2:445342/1‑140 (MQ=255)
                                                                aCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAAATGACCTGGGACTGGTCAAAAGGAGATCATCCTg                                                                >  2:399653/1‑140 (MQ=255)
                                                                    ttAACTTTCACAGCACAAAGATGGATGAAATTGTGCTTTGTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTa                                                            <  2:257799/140‑1 (MQ=255)
                                                                    ttAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTa                                                            >  2:343849/1‑140 (MQ=255)
                                                                            cacaGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtca                                                    <  2:539782/140‑1 (MQ=255)
                                                                            cacaGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtca                                                    <  1:175981/140‑1 (MQ=255)
                                                                            cacaGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtca                                                    <  1:583/140‑1 (MQ=255)
                                                                            cacaGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCAGGGCTTGGTCAAAAGGTGCTCCTCCTGATTATGGTTtca                                                    >  2:208597/1‑140 (MQ=255)
                                                                              caGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCATAAGTAGGTCATCCGGATTATGGTTtcatc                                                  >  1:257799/1‑140 (MQ=255)
                                                                              caGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcatc                                                  >  1:721526/1‑140 (MQ=255)
                                                                                 cacaAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCAtt                                               <  1:708731/140‑1 (MQ=255)
                                                                                 cacaAAGATAGATGAAATCGTGCTTTTTGCCGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCAtt                                               <  1:445342/140‑1 (MQ=255)
                                                                                     aaGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATA‑‑‑‑‑‑‑‑‑TCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCa                                  <  2:644581/140‑1 (MQ=255)
                                                                                      agatagatGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACAt                                          <  1:343849/140‑1 (MQ=255)
                                                                                                aaTCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCCTCCTGATTATGGTTCCATCATACACATCTGCCCCACg                                >  1:289525/1‑140 (MQ=255)
                                                                                                  tCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCa                              >  2:31279/1‑140 (MQ=255)
                                                                                                  tCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCa                              <  2:747617/140‑1 (MQ=255)
                                                                                                   cGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCa                                  <  2:122169/135‑1 (MQ=255)
                                                                                                   cGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCa                                  >  1:122169/1‑135 (MQ=255)
                                                                                                   cGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCa                              >  2:804836/1‑139 (MQ=255)
                                                                                                   cGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCa                              <  1:804836/139‑1 (MQ=255)
                                                                                                     tGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCg                           >  2:775515/1‑140 (MQ=255)
                                                                                                     tGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCg                           <  1:35204/140‑1 (MQ=255)
                                                                                                     tGCTTTTTGCTGTTTTTTCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCg                           >  2:88268/1‑140 (MQ=255)
                                                                                                                      tCGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTgcg          <  1:775515/140‑1 (MQ=255)
                                                                                                                       cGAACATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTgcgc         <  2:88774/140‑1 (MQ=255)
                                                                                                                          aCATATCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAg      <  2:141441/140‑1 (MQ=255)
                                                                                                                              atCCTAACTATCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTg  <  2:598825/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
CCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTG  >  NC_000913/2531256‑2531521

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