Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,515,968 C→G 9.6% G237A (GGA→GCA)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,515,9680CG 9.6% 265.4 / 16.4 83G237A (GGA→GCA) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (40/35);  new base G (3/5);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGG  >  NC_000913/4515820‑4516114
                                                                                                                                                    |                                                                                                                                                  
tgtgCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTcc                                                                                                                                                   <  1:829644/150‑1 (MQ=255)
 gtgCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCg                                                                                                                                                  >  2:703489/1‑150 (MQ=255)
  tgCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGt                                                                                                                                                 <  2:1606891/150‑1 (MQ=255)
   gccgccCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGATGTTGGGCTGAGATGAACTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTa                                                                                                                                                >  2:1270542/1‑150 (MQ=255)
   gccgccCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTa                                                                                                                                                >  2:266400/1‑150 (MQ=255)
    ccgccCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTaa                                                                                                                                               >  1:270163/1‑150 (MQ=255)
     cgccCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAAc                                                                                                                                              <  1:49265/150‑1 (MQ=255)
     cgccCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAAc                                                                                                                                              >  1:994287/1‑150 (MQ=255)
         caccaTCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAAcgcac                                                                                                                                          >  2:923358/1‑150 (MQ=255)
         caccaTCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAAcgcac                                                                                                                                          <  2:1481459/150‑1 (MQ=255)
           ccaTCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAAcgcaccg                                                                                                                                        <  1:1341513/150‑1 (MQ=255)
            caTCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTGCGTTAcggaccgc                                                                                                                                       <  1:1203461/150‑1 (MQ=255)
                aGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTGCGTTAcggaccgcagag                                                                                                                                   <  1:525911/150‑4 (MQ=255)
                aGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTGCGTTAcggaccgcagag                                                                                                                                   <  1:308940/150‑4 (MQ=255)
                aGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTGCGTTAcggaccgcagag                                                                                                                                   <  2:1031462/150‑4 (MQ=255)
                aGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTGCGTTAcggaccgcagag                                                                                                                                   <  2:491819/150‑4 (MQ=255)
                   gtgtgCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGcaccacca                                                                                                                                <  1:384077/149‑1 (MQ=255)
                       gCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACcgc                                                                                                                            <  1:1194030/150‑1 (MQ=255)
                            tctTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTAcc                                                                                                                       >  2:980018/1‑150 (MQ=255)
                            tctTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTAcc                                                                                                                       <  1:294261/150‑1 (MQ=255)
                             ctTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCAc                                                                                                                     <  2:178012/150‑1 (MQ=255)
                              tttCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCAGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCCCCGCGTACCAc                                                                                                                     >  2:612322/1‑150 (MQ=255)
                                 cGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTc                                                                                                                  <  1:574976/150‑1 (MQ=255)
                                  gggTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCa                                                                                                                 >  1:208312/1‑150 (MQ=255)
                                    gttgttTTGTGAAGAGGTTGGGCTGAGCTGACATTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCAGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAAt                                                                                                               >  2:1026538/1‑150 (MQ=255)
                                       gttTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGCAACGCACCGCACCACCACCGCGTACCACGTCAATGGc                                                                                                            >  2:1030884/1‑150 (MQ=255)
                                        ttttGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCa                                                                                                           >  2:34178/1‑150 (MQ=255)
                                               aaGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGt                                                                                                    >  1:926773/1‑150 (MQ=255)
                                                  agGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGc                                                                                                 <  1:315351/150‑1 (MQ=255)
                                                       gggCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGc                                                                                            <  2:133415/150‑1 (MQ=255)
                                                        ggCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCg                                                                                           >  1:634815/1‑150 (MQ=255)
                                                        ggCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCg                                                                                           <  1:1008158/150‑1 (MQ=255)
                                                        ggCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCAGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCg                                                                                           >  2:900444/1‑150 (MQ=255)
                                                         gctgagctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACcgc                                                                                                                            >  1:1034801/1‑116 (MQ=255)
                                                         gctgagctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACcgc                                                                                                                            <  2:1034801/116‑1 (MQ=255)
                                                              gctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACgg                                                                                     <  1:1622165/150‑1 (MQ=255)
                                                               ctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACggg                                                                                    <  2:1300801/150‑1 (MQ=255)
                                                                tgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGa                                                                                   >  2:702704/1‑150 (MQ=255)
                                                                 gaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAg                                                                                  <  2:257730/150‑1 (MQ=255)
                                                                   ccTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGcc                                                                                >  1:998674/1‑150 (MQ=255)
                                                                      tCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGt                                                                             >  2:1588094/1‑150 (MQ=255)
                                                                        cACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGa                                                                           >  1:104851/1‑150 (MQ=255)
                                                                           gcccgccTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAg                                                                        <  1:636454/150‑1 (MQ=255)
                                                                           gcccgccTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGACATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAg                                                                        <  1:1172207/150‑1 (MQ=255)
                                                                           gcccgccTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCGCTGCGGTCCGTAACGCACCGCACCACCACCGCGGACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAg                                                                        <  2:873814/150‑1 (MQ=255)
                                                                              cgccTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTg                                                                     >  1:1339026/1‑150 (MQ=255)
                                                                               gccTCGATACCAAAGTCCTGCGGAATGGCACGGGGAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATCGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGc                                                                    <  2:808426/150‑1 (MQ=255)
                                                                                ccTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCg                                                                   >  1:190374/1‑150 (MQ=255)
                                                                                 cTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCgg                                                                  >  1:262794/1‑150 (MQ=255)
                                                                                      tACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGAc                                                             <  2:1517010/150‑1 (MQ=255)
                                                                                        ccAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCg                                                           >  1:1162994/1‑150 (MQ=255)
                                                                                        ccAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCg                                                           >  1:1296726/1‑150 (MQ=255)
                                                                                        ccAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCg                                                           <  2:246816/150‑1 (MQ=255)
                                                                                          aaaGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTa                                                         >  2:718807/1‑150 (MQ=255)
                                                                                                cTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGc                                                   >  2:656572/1‑150 (MQ=255)
                                                                                                 tGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCg                                                  >  2:1358877/1‑150 (MQ=255)
                                                                                                  gCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGa                                                 >  2:1035613/1‑150 (MQ=255)
                                                                                                   cGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGaa                                                >  1:791754/1‑150 (MQ=255)
                                                                                                        tGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGa                                           <  2:418616/150‑1 (MQ=255)
                                                                                                         ggCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGAc                                          >  2:859847/1‑150 (MQ=255)
                                                                                                          gCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACg                                         <  2:208312/150‑1 (MQ=255)
                                                                                                           cACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACgg                                        >  2:560481/1‑150 (MQ=255)
                                                                                                            aCGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACggg                                       <  1:923358/150‑1 (MQ=255)
                                                                                                                 tAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGc                                  >  1:967977/1‑150 (MQ=255)
                                                                                                                 tAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGc                                  >  1:961794/1‑150 (MQ=255)
                                                                                                                   aCAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCg                                >  1:792372/1‑150 (MQ=255)
                                                                                                                       aGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCa                            <  2:906062/150‑1 (MQ=255)
                                                                                                                          tCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCa                         <  2:634815/150‑1 (MQ=255)
                                                                                                                          tCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCa                         <  1:656572/150‑1 (MQ=255)
                                                                                                                          tCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCa                         <  2:922291/150‑1 (MQ=255)
                                                                                                                           caccacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAg                        >  1:1551117/1‑150 (MQ=255)
                                                                                                                                   cgccCACGCTCTGCTGTCCGTAACGCACCGCACAACCACCGCGTACCAAGTCAATGGCATCCATGTTGCCGCGCGAAACGGGACCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCg                >  2:401953/1‑150 (MQ=255)
                                                                                                                                     cccACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAg              >  1:1202703/1‑150 (MQ=255)
                                                                                                                                       cACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCg            <  2:270163/150‑1 (MQ=255)
                                                                                                                                         cGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGc          >  2:530495/1‑150 (MQ=255)
                                                                                                                                         cGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGc          <  1:702704/150‑1 (MQ=255)
                                                                                                                                          gCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGTCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCt         <  1:703489/150‑1 (MQ=255)
                                                                                                                                             gtGCGGTGCGTTACGGACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGc      >  2:325418/2‑150 (MQ=255)
                                                                                                                                              tGCGGTGCGTTACGGACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCg     >  1:147329/1‑150 (MQ=255)
                                                                                                                                              tGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCg     <  2:1162994/150‑1 (MQ=255)
                                                                                                                                               gCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGa    <  2:104851/150‑1 (MQ=255)
                                                                                                                                                 ggTGCGTTACGGACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAgg  >  1:764995/1‑150 (MQ=255)
                                                                                                                                                 ggTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAgg  >  2:921169/1‑150 (MQ=255)
                                                                                                                                                 ggTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAgg  <  1:34178/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGTGCCGCCCACCATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGG  >  NC_000913/4515820‑4516114

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