Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 611,329 G→A 6.4% T389M (ACG→ATG)  fepA ← ferrienterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913611,3290GA 6.4% 201.7 / 3.6 62T389M (ACG→ATG) fepAferrienterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base A (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-01

GTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAATTTCCCTTCGGTACCGCCCGCCAGACCTTC  >  NC_000913/611181‑611474
                                                                                                                                                    |                                                                                                                                                 
gttTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCg                                                                                                                                                   >  3:481293/1‑149 (MQ=255)
 ttttCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGt                                                                                                                                                  >  1:229932/1‑149 (MQ=255)
 ttttCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACACACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCGTCCGTTGCTGATTCCACTCCGt                                                                                                                                                  >  1:501445/1‑149 (MQ=255)
       gCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCcag                                                                                                                                            >  1:385320/1‑149 (MQ=255)
        cAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCcagc                                                                                                                                           <  2:229932/149‑1 (MQ=255)
         aaaCAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCcagcg                                                                                                                                          >  1:248686/1‑149 (MQ=255)
                 agaAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCagcgtcagcgtc                                                                                                                                  <  2:20596/149‑1 (MQ=255)
                 agaAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCagcgtcagcgtc                                                                                                                                  <  2:574886/149‑1 (MQ=255)
                     aaTTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGt                                                                                                                              >  1:90841/1‑149 (MQ=255)
                        ttCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTaa                                                                                                                           <  4:145629/149‑1 (MQ=255)
                           tGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGa                                                                                                                        <  1:448389/149‑1 (MQ=255)
                           tGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGa                                                                                                                        <  2:245238/149‑1 (MQ=255)
                           tGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGa                                                                                                                        >  4:311700/1‑149 (MQ=255)
                              ttttGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGa                                                                                                                     >  2:183235/1‑149 (MQ=255)
                               tttGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGaa                                                                                                                    <  3:270391/149‑1 (MQ=255)
                                   aaTACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAaatca                                                                                                                >  2:159653/1‑149 (MQ=255)
                                      aCGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAaatcaatc                                                                                                             <  1:504761/149‑1 (MQ=255)
                                       cGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCg                                                                                                            <  1:159653/149‑1 (MQ=255)
                                         ggCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGc                                                                                                          >  3:40382/1‑149 (MQ=255)
                                                    gtggtACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTc                                                                                               >  2:95408/1‑149 (MQ=255)
                                                    gtggtACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTc                                                                                               <  3:537642/149‑1 (MQ=255)
                                                      ggtACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTAATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGc                                                                                             >  4:74005/1‑149 (MQ=255)
                                                       gtACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCt                                                                                            >  2:343542/1‑149 (MQ=255)
                                                         aCTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCtgt                                                                                          <  4:186802/149‑1 (MQ=255)
                                                          cTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCtgtg                                                                                         >  1:66299/1‑149 (MQ=255)
                                                          cTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCATACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCtgtg                                                                                         >  1:242246/1‑149 (MQ=255)
                                                          cTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCATACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCtgtg                                                                                         >  1:497526/1‑149 (MQ=255)
                                                           tCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTgc                                                                                        >  1:215024/1‑149 (MQ=255)
                                                           tCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTgc                                                                                        >  4:31369/1‑149 (MQ=255)
                                                           tCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTgc                                                                                        >  4:109722/1‑149 (MQ=255)
                                                             aCGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTgcag                                                                                      >  2:414425/1‑149 (MQ=255)
                                                                   tCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCAc                                                                                >  4:131427/1‑149 (MQ=255)
                                                                     aaTAGCGCCACCGGTATTCGTTCCGGTCAGTGCCAGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGt                                                                              <  3:311700/149‑1 (MQ=255)
                                                                      aTAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGtc                                                                             <  4:125257/149‑1 (MQ=255)
                                                                         gcgcCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGtcatc                                                                          <  2:16490/149‑1 (MQ=255)
                                                                                 ggTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAGGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGa                                                                  <  2:5991/149‑1 (MQ=255)
                                                                                 ggTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGa                                                                  >  3:542792/1‑149 (MQ=255)
                                                                                       cGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATcga                                                            <  2:385320/149‑1 (MQ=255)
                                                                                        gTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATcgac                                                           >  1:428110/1‑149 (MQ=255)
                                                                                        gTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGCAAATCAATCGGCAGGTTAACTACGATGTGCAGCATCACGTCATCAAGATCGATATcgac                                                           <  4:102054/149‑1 (MQ=255)
                                                                                               tCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGa                                                                                                                        <  2:364839/81‑1 (MQ=255)
                                                                                               tCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGa                                                                                                                        >  1:364839/1‑81 (MQ=255)
                                                                                                cAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACATCATCAAGATCGATATCGACGAAATCCt                                                   <  1:555585/149‑1 (MQ=255)
                                                                                                 aGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCtg                                                  <  1:126674/149‑1 (MQ=255)
                                                                                                               ggAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGt                                    <  3:109722/149‑1 (MQ=255)
                                                                                                                gAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCATACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGtt                                   <  2:242246/149‑1 (MQ=255)
                                                                                                                gAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCATACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGtt                                   <  2:497526/149‑1 (MQ=255)
                                                                                                                   cTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTaaa                                <  1:180727/149‑1 (MQ=255)
                                                                                                                   cTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTaaa                                >  4:395396/1‑149 (MQ=255)
                                                                                                                    ttAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAAt                               >  4:216765/1‑149 (MQ=255)
                                                                                                                    ttAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAAt                               >  2:364762/1‑149 (MQ=255)
                                                                                                                    ttAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGGTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAAt                               >  2:530775/1‑149 (MQ=255)
                                                                                                                      aaGTCCTTCATCCGTTGCTGATTCCAGTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAAttt                             >  4:484933/1‑149 (MQ=255)
                                                                                                                        gTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATGGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAATTTcc                           >  3:352434/1‑149 (MQ=255)
                                                                                                                        gTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAATTTcc                           <  3:460126/149‑1 (MQ=255)
                                                                                                                        gTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAATTTcc                           <  3:412075/149‑1 (MQ=255)
                                                                                                                        gTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAATTTcc                           <  3:216765/149‑1 (MQ=255)
                                                                                                                            ttCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGATGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCACCAAGATCGATATCGACGAAATCCTGTGTCGCTCTTTCGTTAAATTTCCCTTc                       <  4:481191/149‑1 (MQ=255)
                                                                                                                             tCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAATTTCCCTTCg                      <  1:183235/149‑1 (MQ=255)
                                                                                                                                  cGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTGTTTCGTTAAATTTACCTTCGGTAcc                 >  3:371230/1‑149 (MQ=255)
                                                                                                                                        tGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAATTTCCCTTCGGTAccgcccgc           >  3:420157/1‑149 (MQ=255)
                                                                                                                                               aCTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAATTTCCCTTCGGTACCGCCCGCCAGACCt    >  2:45626/1‑149 (MQ=255)
                                                                                                                                                 tCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAATTTCCCTTCGGTACCGCCCGCCAGACCTTc  >  2:26817/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCATCAAGATCGATATCGACGAAATCCTGTGTCGCTTTTTCGTTAAATTTCCCTTCGGTACCGCCCGCCAGACCTTC  >  NC_000913/611181‑611474

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