Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 622,725 A→C E142D (GAA→GAC entS → enterobactin exporter, iron‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913622,7250AC98.3% 206.1 / ‑3.3 60E142D (GAA→GACentSenterobactin exporter, iron‑regulated
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (33/25);  total (35/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 = 7.34e-01

ATTGGCTTCATTGGACTGTGCCTTAATGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGAAAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACGGGCTGGCGGCGGCGGGCACGTTTATTACC  >  NC_000913/622579‑622866
                                                                                                                                                  |                                                                                                                                             
aTTGGCTTCATTGGACTGTGCCTTAATGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACt                                                                                                                                           >  2:495979/1‑151 (MQ=255)
aTTGGCTTCATTGGACTGTGCCTTAATGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACt                                                                                                                                           >  1:451200/1‑151 (MQ=255)
    gCTTCATTGGACTGTGCCTTAATGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAAt                                                                                                                                       <  1:1115881/151‑1 (MQ=255)
     cTTCATTGGACTGTGCCTTAATGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAAt                                                                                                                                       >  2:1131014/1‑150 (MQ=255)
      ttCTTTGGACTGGGCCTTAATGCACTGCTGCCGGAGCCGTCAAGGCTGGCAATCTATTTCCTTGGTTTATGGGATGGTTTTTTCGCCTAGCTTGGCGTTACGGCGCTATTGGCCTCGACACCCGCCATGTTCGGGCGTGACGACTTATTgt                                                                                                                                     >  2:716266/1‑150 (MQ=255)
       tCATTGGACTGTGCCTTAATGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATg                                                                                                                                      <  2:1130795/149‑1 (MQ=255)
         aTTGGACTGTGCCTTAATGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAgg                                                                                                                                  <  1:883733/151‑1 (MQ=255)
                 gtgCCTTAATGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGc                                                                                                                           <  1:987940/150‑1 (MQ=255)
                 gtgCCTTAATGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCg                                                                                                                          <  2:66288/151‑1 (MQ=255)
                    ccTTAATGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATc                                                                                                                       >  1:771617/1‑151 (MQ=255)
                        aaTGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCAcc                                                                                                                    >  1:951447/1‑150 (MQ=255)
                           gCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGt                                                                                                                >  2:1006810/1‑151 (MQ=255)
                           gCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGt                                                                                                                >  2:202221/1‑151 (MQ=255)
                                  tgcCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGt                                                                                                         >  1:939306/1‑151 (MQ=255)
                                    ccGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGc                                                                                                       <  2:852283/151‑1 (MQ=255)
                                      ggAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCg                                                                                                      >  2:312548/1‑150 (MQ=255)
                                      ggAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGt                                                                                                     <  2:743579/151‑1 (MQ=255)
                                      ggAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGt                                                                                                     >  1:971106/1‑151 (MQ=255)
                                       gAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTATCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGTGCGATCACCATGTTGGCCGTGCg                                                                                                      <  2:827298/149‑1 (MQ=255)
                                           cGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTggg                                                                                                >  1:908841/1‑151 (MQ=255)
                                           cGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTggg                                                                                                <  1:838792/151‑1 (MQ=255)
                                             tCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTgggg                                                                                               >  1:1113392/1‑150 (MQ=255)
                                             tCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGt                                                                                              >  1:368937/1‑151 (MQ=255)
                                              cATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTc                                                                                             <  2:908780/151‑1 (MQ=255)
                                                       cAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGAttt                                                                                     >  2:785044/1‑150 (MQ=255)
                                                              tttACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATg                                                                             >  1:695617/1‑151 (MQ=255)
                                                              tttACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGAAAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATg                                                                             >  1:733950/1‑151 (MQ=255)
                                                               ttACTTGGTTTATGGGATGGTTTTTTCGCATCCCTTGGCGTTACGGCGCTATTGGCGGCGACACCCGCACTGGGAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCCGCCATg                                                                             >  2:862867/1‑150 (MQ=255)
                                                                 aCTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGAtt                                                                          >  2:101083/1‑151 (MQ=255)
                                                                 aCTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGAtt                                                                          <  1:123685/151‑1 (MQ=255)
                                                                    tGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTgg                                                                        <  2:344926/150‑1 (MQ=255)
                                                                       tttATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCgg                                                                     >  1:216976/1‑150 (MQ=255)
                                                                            gggATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGGTTAtt                                                               >  2:178737/1‑151 (MQ=255)
                                                                            gggATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGATAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTAtt                                                               <  1:196592/151‑1 (MQ=255)
                                                                             ggATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTg                                                              <  2:290184/151‑1 (MQ=255)
                                                                               aTGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGAGACACCAGCACTGGGAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGAACGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGc                                                             >  1:58118/1‑150 (MQ=255)
                                                                                 ggTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTgg                                                          <  2:1053214/151‑1 (MQ=255)
                                                                                       ttCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGAcc                                                     <  1:75898/150‑1 (MQ=255)
                                                                                        tCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGAcc                                                     <  2:451795/149‑1 (MQ=255)
                                                                                         cGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGt                                                  <  2:1177461/151‑1 (MQ=255)
                                                                                         cGCATCGCTTGGCGTTACGGAGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCg                                                    >  2:553083/1‑149 (MQ=255)
                                                                                           cATCGCTTGGCGTTACGGCGCTATTGGCGGCGAACCCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTGGCCCATGCTTGGCGGTTTTCTGCCAGCGGCCGgggg                                                >  2:951120/1‑151 (MQ=255)
                                                                                            aTCGCTTGGCGTTACGGCGCTATTGGCGGCTACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGc                                               <  1:671662/151‑1 (MQ=255)
                                                                                            aTCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGAGAACTTAATGCAGGCCGGGGCGCTCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTTTTGCCGGCGACCCGTgg                                                >  1:1125921/1‑150 (MQ=255)
                                                                                            aTCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGc                                               >  1:1091546/1‑151 (MQ=255)
                                                                                               gCTTGGCGCTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGGGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGAACGGGGGCGt                                             >  2:1186643/1‑150 (MQ=255)
                                                                                                cTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCg                                              >  1:877992/1‑148 (MQ=255)
                                                                                                   ggCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCt                                        <  1:1161206/151‑1 (MQ=255)
                                                                                                             cgcTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACgg                              <  2:942502/151‑1 (MQ=255)
                                                                                                              gcTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACggg                             <  2:708181/151‑1 (MQ=255)
                                                                                                               cTATTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACggg                             >  2:732593/1‑150 (MQ=255)
                                                                                                                 aTTGGCGGCGACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACGGGCt                           <  1:239569/150‑1 (MQ=255)
                                                                                                                      cggcgACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACGGGCTggcggc                     >  1:735297/1‑151 (MQ=255)
                                                                                                                          gACACCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCTCCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTTCGGGCTGgcggcggcg                 <  2:956114/151‑1 (MQ=255)
                                                                                                                            cacCAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACGGGCTGGCGGCGGCggg               >  1:735955/1‑151 (MQ=255)
                                                                                                                               cAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACGGGCTGGCGGCGGCggggcc            >  2:886888/1‑148 (MQ=255)
                                                                                                                               cAGCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACGGGCTGGCGGCGGCGGGCAc            >  1:932083/1‑151 (MQ=255)
                                                                                                                                 gCACTGGTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACGGGCTGGCGGCGGCGGGCACGt          >  2:491473/1‑151 (MQ=255)
                                                                                                                                      ggTAGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGGGGCGGAGCCTGGAACTACGGGCTGGCGGCGGCGGGCACGTttatt     >  2:189333/1‑151 (MQ=255)
                                                                                                                                         aGGGCGTGACAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACGGGCTGGCGGCGGCGGGCACGTTTATTAcc  <  1:319868/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
ATTGGCTTCATTGGACTGTGCCTTAATGCACTGCTGCCGGAGCCGTCATTGCTGGCAATCTATTTACTTGGTTTATGGGATGGTTTTTTCGCATCGCTTGGCGTTACGGCGCTATTGGCGGCGACACCAGCACTGGTAGGGCGTGAAAACTTAATGCAGGCCGGGGCGATCACCATGTTGACCGTGCGTCTGGGGTCGGTGATTTCGCCCATGATTGGCGGTTTATTGCTGGCGACCGGTGGCGTAGCCTGGAACTACGGGCTGGCGGCGGCGGGCACGTTTATTACC  >  NC_000913/622579‑622866

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