Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,521,409 C→T 5.0% A553T (GCA→ACA)  yncD ← putative iron outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,521,4090CT 5.0% 256.8 / 2.1 79A553T (GCA→ACA) yncDputative iron outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (38/37);  new base T (2/2);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

TCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAAGCAGCCCATCACCAATACGCGTTTTAC  >  NC_000913/1521262‑1521552
                                                                                                                                                   |                                                                                                                                               
tCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGc                                                                                                                                                 <  3:324880/148‑1 (MQ=255)
tCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCAATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGc                                                                                                                                                 <  1:216149/148‑1 (MQ=255)
    cGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCTTTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGtt                                                                                                                                            <  2:167177/149‑1 (MQ=255)
    cGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGtt                                                                                                                                            >  4:373021/1‑149 (MQ=255)
     gCAACGCCCATTTTACGGGCGATCCCTGGCATCCGAGTACCGTTACACTCCTGTTCATTGCAAACATTGCTGCGATAGGTCACATCAAGCCCGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTc                                                                                                                                           >  1:252675/1‑149 (MQ=255)
      cAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCa                                                                                                                                          <  5:386796/149‑1 (MQ=255)
       aaaGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAg                                                                                                                                         >  8:231054/1‑149 (MQ=255)
       aaaGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCAAGGTCCAGGACGCGGTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAg                                                                                                                                         >  5:361702/1‑149 (MQ=255)
              atatTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCtt                                                                                                                                   >  7:139412/1‑148 (MQ=255)
                atTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTgac                                                                                                                                >  7:282342/1‑149 (MQ=255)
                     gggCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACgg                                                                                                                           <  4:127556/149‑1 (MQ=255)
                      ggCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACggg                                                                                                                          <  7:231054/149‑1 (MQ=255)
                           tCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCttt                                                                                                                     >  6:2025/1‑149 (MQ=255)
                             ccTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTc                                                                                                                    <  6:372543/148‑1 (MQ=255)
                              cTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCg                                                                                                                  <  5:298006/149‑1 (MQ=255)
                              cTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCg                                                                                                                  <  7:375087/149‑1 (MQ=255)
                               tGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCgg                                                                                                                 >  8:155769/1‑149 (MQ=255)
                                   aTCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACgg                                                                                                                           <  8:322571/135‑1 (MQ=255)
                                   aTCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACgg                                                                                                                           >  7:322571/1‑135 (MQ=255)
                                      cGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCAttttt                                                                                                          >  4:97307/1‑149 (MQ=255)
                                      cGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCAttttt                                                                                                          <  1:144082/149‑1 (MQ=255)
                                      cGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCAttttt                                                                                                          <  3:97307/149‑1 (MQ=255)
                                      cGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCAttttt                                                                                                          <  3:373021/149‑1 (MQ=255)
                                      cGATTACAGTTACACTCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACTTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCAttttt                                                                                                          <  5:95827/149‑1 (MQ=255)
                                              gTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCttt                                                                                                                     >  5:149929/1‑130 (MQ=255)
                                              gTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCttt                                                                                                                     <  6:149929/130‑1 (MQ=255)
                                                ttctttCCTGTTCTTTTGAAACATTGCTGGGATAGGTAGCATCTAGCCAGGTCCAGGACGAGTTTACGCGATATTGTCGTGCGAAACGTTGATCCCATTCCAGTTCAGAGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCTTa                                                                                                <  6:14310/144‑1 (MQ=255)
                                                  cAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACTGGTCTTTCCGGCATTTTTGTAAGTCGTAc                                                                                               <  3:208663/148‑1 (MQ=255)
                                                   aaTCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGc                                                                                             >  8:97297/1‑149 (MQ=255)
                                                     tCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCTGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTa                                                                                                >  5:14310/1‑144 (MQ=255)
                                                       cTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCAc                                                                                         <  7:336160/149‑1 (MQ=255)
                                                        tGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCAAGGACGCGTTTAAGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGACTTGACGACTGGTCTTTACGGCATTTTTGTAAGTCGTACGaccacc                                                                                        >  8:316430/1‑149 (MQ=255)
                                                                gCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTa                                                                                    >  1:86562/1‑145 (MQ=255)
                                                                gCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTa                                                                                    <  2:86562/145‑1 (MQ=255)
                                                                gCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGc                                                                                >  3:147135/1‑149 (MQ=255)
                                                                  aaaCATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTa                                                                              >  4:44764/1‑149 (MQ=255)
                                                                       ttGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTGTCGac                                                                         <  5:2025/149‑1 (MQ=255)
                                                                       ttGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGac                                                                         <  3:44764/149‑1 (MQ=255)
                                                                           tgcGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGacaac                                                                      >  7:241181/1‑148 (MQ=255)
                                                                              gATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAAtt                                                                  >  7:161333/1‑149 (MQ=255)
                                                                                 aGGTGGCATGAAGCCAGGTCCAGGACGCGTTTTCGCGAAAATCTCCTTCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTtc                                                                <  8:161333/148‑1 (MQ=255)
                                                                                 aGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGTCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTg                                                                                 <  5:190703/131‑1 (MQ=255)
                                                                                 aGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGTCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTg                                                                                 <  5:196633/131‑1 (MQ=255)
                                                                                 aGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGTCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTg                                                                                 >  6:196633/1‑131 (MQ=255)
                                                                                 aGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGTCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTg                                                                                 >  6:190703/1‑131 (MQ=255)
                                                                                   gTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTgacgac                                                                                                                             >  6:61016/1‑85 (MQ=255)
                                                                                   gTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTgacgac                                                                                                                             <  5:61016/85‑1 (MQ=255)
                                                                                    tCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTtcatca                                                            >  2:410690/1‑149 (MQ=255)
                                                                                       cATCTAGCCAGGTCCAGGACGCGTTTACGCGAAACTCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCGATCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACACTTccctcatca                                                         <  2:252675/149‑1 (MQ=255)
                                                                                         tCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAgt                                                       >  8:42377/1‑149 (MQ=255)
                                                                                         tCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAgt                                                       >  1:222206/1‑149 (MQ=255)
                                                                                             gCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCg                                                   >  6:109606/1‑149 (MQ=255)
                                                                                               cAGGTCCAGGACGCGTTTACGCGAAAATCTCCAGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGt                                                 >  5:2861/1‑149 (MQ=255)
                                                                                                  gTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAACCGTTGATCCCATGCCAGGTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCCCCGCTACTGCTATCGACAACAATTTCATCATCATGGTCGGTCt                                               >  3:21374/1‑148 (MQ=255)
                                                                                                     cAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGTGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGaaa                                           <  8:282342/149‑1 (MQ=255)
                                                                                                      aGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGaaa                                           >  1:77411/1‑148 (MQ=255)
                                                                                                             gTTTACGCGAAAATCTCCTGCGAAAGGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGTGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCga                                   <  7:316430/149‑1 (MQ=255)
                                                                                                             gTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGc                                     >  4:382898/1‑147 (MQ=255)
                                                                                                             gTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCga                                   <  4:147135/149‑1 (MQ=255)
                                                                                                                tACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCgaga                                 <  2:77411/148‑1 (MQ=255)
                                                                                                                tACGCGAAAATCTCCAGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCgaga                                 <  6:2861/148‑1 (MQ=255)
                                                                                                                   gcgAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCg                                                   >  4:330272/1‑127 (MQ=255)
                                                                                                                   gcgAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCg                                                   <  3:330272/127‑1 (MQ=255)
                                                                                                                   gcgAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACt                               >  6:406137/1‑147 (MQ=255)
                                                                                                                    cgAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAag                            <  3:382898/149‑1 (MQ=255)
                                                                                                                     gAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTaa                             <  7:101888/147‑1 (MQ=255)
                                                                                                                     gAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTaa                             >  8:101888/1‑147 (MQ=255)
                                                                                                                     gAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAag                            >  8:334552/1‑148 (MQ=255)
                                                                                                                        aaTCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAagcagc                        <  7:97297/149‑1 (MQ=255)
                                                                                                                            tcCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAAGCAGCCCAt                    >  5:121189/1‑149 (MQ=255)
                                                                                                                              cTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAAGCAGCCCATc                   <  1:192819/148‑1 (MQ=255)
                                                                                                                                  gAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAAGCAGCCCATCACCaa              >  2:194075/1‑149 (MQ=255)
                                                                                                                                     aCGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAAGCAGCCCATCACCAATAc           <  1:258427/149‑1 (MQ=255)
                                                                                                                                          gATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAAGCAGCCCATCACCAATACGCGtt      >  6:155555/1‑149 (MQ=255)
                                                                                                                                           aTCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCAttttt                                                                                                          <  3:288669/48‑1 (MQ=255)
                                                                                                                                           aTCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCAttttt                                                                                                          >  4:288669/1‑48 (MQ=255)
                                                                                                                                           aTCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACt                               <  3:68171/123‑1 (MQ=255)
                                                                                                                                           aTCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACt                               >  4:68171/1‑123 (MQ=255)
                                                                                                                                           aTCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAAGCAGCCCATCACCAATACGCGttt     >  7:281357/1‑149 (MQ=255)
                                                                                                                                              ccATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAAGCAGCCCATCACCAATACGCGTTTTa   <  8:281357/148‑1 (MQ=255)
                                                                                                                                              ccATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAAGCAGCCCATCACCAATACGCGTTTTAc  <  7:42377/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TCGACGCAAAGCCCATATTACGGGCGATCCCTGGCATCCGATTACCGTTACAATCCTGTTCATTGCAAACATTGCTGCGATAGGTCGCATCAAGCCAGGTCCAGGACGCGTTTACGCGAAAATCTCCTGCGAAACGTTGATCCCATGCCAGTTCAGCGCCTTGACGACGGGTCTTTCCGGCATTTTTGTAAGTCGTACGCCCACCGCTACTGCTATCGACAACAATTTCATCATCAGTGTCGGTCTGAAACAATGCGAGACTAAGCAGCCCATCACCAATACGCGTTTTAC  >  NC_000913/1521262‑1521552

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