Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,605,791 A→G intergenic (‑128/‑13) fhuF ← / → yjjZ ferric‑siderophore reductase FhuF/DUF1435 domain‑containing protein YjjZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,605,7910AG97.5% 129.1 / ‑4.3 40intergenic (‑128/‑13)fhuF/yjjZferric‑siderophore reductase FhuF/DUF1435 domain‑containing protein YjjZ
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (19/20);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

TCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAG  >  NC_000913/4605570‑4606034
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
tcgCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCt                                                                                                                                                                                                                        >  1:76124/1‑251 (MQ=255)
    ccGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGc                                                                                                                                                                                                                    >  1:299925/1‑251 (MQ=255)
                                                  cGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTg                                                                                                                                                                       <  1:478187/250‑1 (MQ=255)
                                                   gCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATAATTATCAATGAGAGAGATGAAATCATGTTGCAACGTAAGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGGTTGCTGATTGCAGGGCTGAt                                                                                                                                                                     >  1:194560/1‑251 (MQ=255)
                                                          gACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCg                                                                                                                                                              <  1:417604/251‑1 (MQ=255)
                                                             aTCTTCATAGAGCGGTGCGGAACGCTAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGGGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGGTGGATGCGGAt                                                                                                                                                           >  2:283049/1‑251 (MQ=255)
                                                               cTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGAttt                                                                                                                                                         >  1:60767/1‑251 (MQ=255)
                                                                         cggggCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTc                                                                                                                                                   >  2:70257/5‑247 (MQ=255)
                                                                                              aaTCGGGATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGGTGTATGCGGATTTATCGTCAGATCAGTGGCGGGTTGTCATTCtt                                                                                                                          >  2:450046/1‑250 (MQ=255)
                                                                                                   ggATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGGTTgggatcggaagagg                                                                                                                     >  2:163162/1‑238 (MQ=255)
                                                                                                   ggATAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTg                                                                                                                                  <  1:163162/238‑1 (MQ=255)
                                                                                                      tAGTAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAACATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGCTTTATCGTCAGATCAGTGGCGGATTGTCCTTCTGGTggggtt                                                                                                                  >  2:130377/1‑248 (MQ=255)
                                                                                                         tAATCTAAATGATAATGATTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTACTTCTTATGGGGTTAGt                                                                                                               >  2:284049/1‑251 (MQ=255)
                                                                                                                  tgataatgatTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGGTGGGGTTAGTATTGGCGcc                                                                                                      >  1:332461/1‑251 (MQ=255)
                                                                                                                  tgataatgatTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGcc                                                                                                      >  1:55007/1‑251 (MQ=255)
                                                                                                                   gataatgatTGCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCg                                                                                                     <  1:450046/251‑1 (MQ=255)
                                                                                                                             gCTAATCCTAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGt                                                                                           <  1:418404/251‑1 (MQ=255)
                                                                                                                                        cGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTgg                                                                                                                                         >  1:470547/1‑194 (MQ=255)
                                                                                                                                        cGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTgg                                                                                                                                         <  2:470547/194‑1 (MQ=255)
                                                                                                                                            aGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGc                                                                            <  1:295634/251‑1 (MQ=255)
                                                                                                                                                   cccGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGGTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTa                                                                     <  2:76124/251‑1 (MQ=255)
                                                                                                                                                    ccGATAGCAAGGTATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGGTATCTGTATCAAGGAGAGAGATGAAATCAGGTTTCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGAGGATTTATCCTTAGAGCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGAGGATGCTGTATCACAAACATTTGCTGCATTat                                                                    <  2:21620/251‑2 (MQ=255)
                                                                                                                                                     cGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACa                                                                   <  1:283049/251‑1 (MQ=255)
                                                                                                                                                                         ggCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTg                                               <  1:329180/251‑1 (MQ=255)
                                                                                                                                                                           cTTGCAAATGATAAAAATTATCATATGCTATTGGATATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGt                                                                                                                                                                   >  1:126716/1‑133 (MQ=255)
                                                                                                                                                                           cTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGt                                                                                                                                                                   <  2:126716/133‑1 (MQ=255)
                                                                                                                                                                                 aattaTAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTTTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTAt                                       <  1:439357/247‑1 (MQ=255)
                                                                                                                                                                                  aTGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTAtt                                      >  1:80118/1‑251 (MQ=255)
                                                                                                                                                                                    gATAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGCTTGTCATTCTGATGGGATTAGTATTGACGCCCATGATGCTGTCTCACAAACAGTTGCGGCATTCCATTTTGCTACCATCGTGCCTGGCACTTATTGc                                    >  1:18443/1‑251 (MQ=255)
                                                                                                                                                                                      tAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCgatgat                                                                                                <  2:456198/189‑1 (MQ=255)
                                                                                                                                                                                      tAAAAATTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCgatgat                                                                                                >  1:456198/1‑189 (MQ=255)
                                                                                                                                                                                       aaaaaTTATCATATGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTgg                                 <  2:100157/251‑1 (MQ=255)
                                                                                                                                                                                                 atatGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTTCCATCGTGCCTGGCACTTATTGCTGGCATCATGCtg                       <  2:425294/251‑1 (MQ=255)
                                                                                                                                                                                                  tattatatTGGTTATCATTATCAAGGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGTAGTGTTGCGGCCGGGATTTCTGATTGCAGGGCTGATGTATGCGTATTTATCGTCAGTTCAGTTGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACATTTGCGGCATTACATTTTGCTTCCATCGTGCCTGGCACTGATTGCTGGCATCATg                          <  1:301684/243‑1 (MQ=255)
                                                                                                                                                                                                   atGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGGTGGGATTAGTATTGACGCCGATGATGCTGTATCCCAAACCGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATg                          >  2:432816/1‑246 (MQ=255)
                                                                                                                                                                                                   atGATATTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATg                          <  1:432816/246‑1 (MQ=255)
                                                                                                                                                                                                        atTGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATa                 >  1:321253/1‑250 (MQ=255)
                                                                                                                                                                                                          tGGTTATCATTATCAATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAAt               <  2:321253/250‑1 (MQ=255)
                                                                                                                                                                                                            gTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGAGCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGaa            <  2:252208/251‑1 (MQ=255)
                                                                                                                                                                                                                       caATGAGAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAg  <  2:299925/250‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
TCGCCCGCACGGCTTGCGCAAGCGTGGGATCCTGCGGCTGGAGATGCGTTCGCCAGATGACATCTTCATAGAGCGGTGCGGAACGATAGGCCATAATCGGGATAGTAATCTAAATGATAATGATTGCTAATCATAGCGATAGGTTTACCCGATAGCAAGGGATTTATCTGGCTTGCAAATGATAAAAATTATCATATGATATTGGTTATCATTATCAATGAAAGAGATGAAATCATGTTGCAACGTACGCTGGGCAGTGGCTGGGGAGTGTTGCTGCCGGGATTGCTGATTGCAGGGCTGATGTATGCGGATTTATCGTCAGATCAGTGGCGGATTGTCATTCTGATGGGATTAGTATTGACGCCGATGATGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAG  >  NC_000913/4605570‑4606034

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