Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 169,749 T→A 6.1% intergenic (+22/‑29) fhuA → / → fhuC ferrichrome outer membrane transporter/iron(3+)‑hydroxamate import ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913169,7490TA 6.1% 422.0 / 8.8 142intergenic (+22/‑29)fhuA/fhuCferrichrome outer membrane transporter/iron(3+)‑hydroxamate import ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (68/65);  new base A (2/7);  total (70/72)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

GCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGACCG  >  NC_000913/169600‑169898
                                                                                                                                                     |                                                                                                                                                     
  tCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTc                                                                                                                                                     <  1:1683787/150‑1 (MQ=255)
   ccAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTcc                                                                                                                                                    >  2:1769055/1‑150 (MQ=255)
   ccAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTcc                                                                                                                                                    >  2:774577/1‑150 (MQ=255)
     aaCGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGATTACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt                                                                                                                                                  >  1:2058723/1‑150 (MQ=255)
          ggCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTGcccc                                                                                                                                             <  1:1892135/150‑2 (MQ=255)
          ggCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTAGGGGCACGGATTTCCGTGCCCa                                                                                                                                             >  1:1416800/1‑150 (MQ=255)
           gcgcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAAATCCGTGcccca                                                                                                                                            <  2:792718/150‑3 (MQ=255)
            cgcgGCATGTTAACAACCTGTTCGATCGTGAATACGGCGCCAGCTGCTTTAGCACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTACTTTCTCTTTTGGGGCACGGAAATCCGTGccccaa                                                                                                                                           <  1:2092316/150‑4 (MQ=255)
            cgcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAtt                                                                                                                                           >  1:220999/1‑150 (MQ=255)
              cTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTc                                                                                                                                         <  2:168168/150‑1 (MQ=255)
              cTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGCACGGAAATCCGTGccccaaaag                                                                                                                                        <  1:1737400/150‑7 (MQ=255)
               tGCATGTTAACAACCTGTTCGATCCTGAATCCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGGCAGGTCGTTGCAACCGCAACCATCCGAATCTAATTTCTCTTTTGGGGCACGGAAATCCGTGccccaaaag                                                                                                                                        <  2:146204/150‑7 (MQ=255)
                   tcttAACAACCTGTTCCATCGAGAATACGTCGCCAGCTGCTTTAACACTTATTGCTGCTGCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCATCCTTCCGTTTCTAATTTCTCTTTTGTGGCACGGATTTCCGTGCCCATTTCACAAg                                                                                                                                    <  2:1070671/148‑1 (MQ=255)
                   tGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAg                                                                                                                                    >  2:561301/1‑150 (MQ=255)
                   tGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAg                                                                                                                                    >  2:25368/1‑150 (MQ=255)
                    gTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGt                                                                                                                                   >  1:1025018/1‑150 (MQ=255)
                       aacaacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg                                                                                                                                >  2:1098965/1‑150 (MQ=255)
                        acaacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGc                                                                                                                               >  2:1374739/1‑150 (MQ=255)
                        acaacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGc                                                                                                                               <  1:1229668/150‑1 (MQ=255)
                         caacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCt                                                                                                                              <  2:1249304/150‑1 (MQ=255)
                          aacCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTACCACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTg                                                                                                                             <  2:1808449/150‑1 (MQ=255)
                             cTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTa                                                                                                                          >  1:542919/1‑150 (MQ=255)
                             cTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTa                                                                                                                          <  1:604830/150‑1 (MQ=255)
                               gTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATg                                                                                                                        >  1:1625696/1‑150 (MQ=255)
                                 tCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTATTTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCa                                                                                                                      >  1:1945118/1‑150 (MQ=255)
                                  cGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAAGCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAg                                                                                                                     <  1:906635/150‑1 (MQ=255)
                                   gATCGTGACTACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTTCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAgg                                                                                                                    <  1:903809/150‑1 (MQ=255)
                                    aTCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGa                                                                                                                   >  1:1548952/1‑150 (MQ=255)
                                     tCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGaa                                                                                                                  <  1:2029736/150‑1 (MQ=255)
                                       gttaaTACGTCGCAAGCTGCTTTAACACTTTTGGCTGCTTCTGGGGCGAAGAACGTCAGGTCGTTTCAAGCGCAACCTTCCCTTTCTAATTTCTCTTTTGGGGCACGGTTTTCCGTGCCCATTTCAGAAGTTGGCTGTTATGCATTAATa                                                                                                                <  2:2073107/147‑1 (MQ=255)
                                        ggAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATac                                                                                                               <  2:77289/149‑1 (MQ=255)
                                            tACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGa                                                                                                           <  2:862835/150‑1 (MQ=255)
                                              cGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCTACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACATGTTGGCGGTTTTGCAGGAATACACGAAt                                                                                                         >  2:414919/1‑150 (MQ=255)
                                              cGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAAt                                                                                                         >  1:957284/1‑150 (MQ=255)
                                               gTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATc                                                                                                        <  2:820484/150‑1 (MQ=255)
                                                  gCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCAtt                                                                                                     >  1:2057386/1‑150 (MQ=255)
                                                  gCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCAtt                                                                                                     >  1:40848/1‑150 (MQ=255)
                                                  gCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGGTGGCTGTTATGCAGGAATACACGAATCAtt                                                                                                     >  1:48207/1‑150 (MQ=255)
                                                   ccAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTc                                                                                                    <  2:625663/150‑1 (MQ=255)
                                                     aGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCg                                                                                                  <  2:1785653/150‑1 (MQ=255)
                                                       ctgctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGAt                                                                                                <  2:608994/150‑1 (MQ=255)
                                                        tgctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATa                                                                                               >  1:2145172/1‑150 (MQ=255)
                                                         gctTTAACACTTATGGCTGCTTCTGTGGCGCAGAACGTCAGGTCGTTGCAACAGCAACATTAAGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCACATTTCACAAGTTGTCTGTTATGCAGGAATACACGAATCATTCCGATAc                                                                                              <  1:1435613/150‑1 (MQ=255)
                                                         gctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGTGGCACGGATTTCCGTGCCCATTTCACAAGTAGGCTGTTATGCAGGAATACACGAATCATTCCGATAc                                                                                              >  1:1059844/1‑150 (MQ=255)
                                                          ctTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCACAGTTCACCAGTTGGCTGTTATGCAGGAATACACGAATAAATcccctaac                                                                                             >  2:1116360/1‑144 (MQ=255)
                                                            ttAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCAc                                                                                           >  2:1620452/1‑150 (MQ=255)
                                                              aaCACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACtt                                                                                         >  2:938475/1‑150 (MQ=255)
                                                              aaCACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACtt                                                                                         <  1:25368/150‑1 (MQ=255)
                                                              aaCACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCAATCCGAAACCACtt                                                                                         >  1:692799/1‑150 (MQ=255)
                                                                cacTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACtttt                                                                                       >  1:1845367/1‑150 (MQ=255)
                                                                 acTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTg                                                                                      <  2:557542/150‑1 (MQ=255)
                                                                  cTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGTGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGc                                                                                     <  1:1582835/150‑1 (MQ=255)
                                                                  cTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTTGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGc                                                                                     <  2:1567640/150‑1 (MQ=255)
                                                                           gctTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTaa                                                                            <  2:489565/150‑1 (MQ=255)
                                                                            ctTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAat                                                                           >  2:373958/1‑150 (MQ=255)
                                                                            ctTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCAGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGATTCATTCCGATACCACTTTTGCACTGCGTAat                                                                           <  1:249230/150‑1 (MQ=255)
                                                                             ttCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAata                                                                          <  1:759085/150‑1 (MQ=255)
                                                                             ttCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAata                                                                          <  1:2027583/150‑1 (MQ=255)
                                                                              tCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAatat                                                                         >  1:1841746/1‑150 (MQ=255)
                                                                               cTGGGGCGCAGAACGTCAGCTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAATTCACAAGTTGGCTGTTAAGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGAAataac                                                                        >  2:934106/1‑148 (MQ=255)
                                                                                 ggggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATATCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATAtctc                                                                      >  2:339777/1‑150 (MQ=255)
                                                                                  gggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTcc                                                                     <  1:1332244/150‑1 (MQ=255)
                                                                                   ggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCt                                                                    <  2:2201836/150‑1 (MQ=255)
                                                                                      gcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTc                                                                 >  2:1451600/1‑150 (MQ=255)
                                                                                      gcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATATCCTTTc                                                                 >  2:1599536/1‑150 (MQ=255)
                                                                                       cAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCg                                                                <  1:711288/150‑1 (MQ=255)
                                                                                       cAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCg                                                                <  2:1841746/150‑1 (MQ=255)
                                                                                       cAGAACGTCAGGTCGTTGAAACCGCAACCTTCAGTTTCTAATTTCTCTTTTGGGGCTCGGATTTCCGTGCCCATTTCACAAGTTTGCTGATATGCAGGAATACACGAATCAATCCGATACCACTTTTGCACTGCGTTATATCTCCTTTCg                                                                >  2:507587/1‑150 (MQ=255)
                                                                                         taaCGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCgtg                                                              <  2:740152/149‑1 (MQ=255)
                                                                                           aCGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGtgtg                                                            >  1:309131/1‑150 (MQ=255)
                                                                                           aCGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGtgtg                                                            >  2:1289657/1‑150 (MQ=255)
                                                                                           aCGTCAGGTCGTTGCAACCGCAACCGTCGGTTTCTAATTTCTCTTTTGGGGCCCGGATTTCCGTGACCATTTCACAAGTTCGCTGTTATGCAGGAATACACGAATCATTCCGATACCTCTTTTGCACTGCGTAATATCTCCTTTCGtgtg                                                            <  1:344742/150‑1 (MQ=255)
                                                                                               cAGGTCGTTGAAACCGGAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCg                                                        <  2:1422177/150‑1 (MQ=255)
                                                                                                          aaccgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGCTACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtt                                             >  1:1206780/1‑150 (MQ=255)
                                                                                                           accgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTTACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCttt                                            <  2:2205957/150‑1 (MQ=255)
                                                                                                            ccgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCtttt                                           <  1:1769055/150‑1 (MQ=255)
                                                                                                             cgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTg                                          <  1:360056/150‑1 (MQ=255)
                                                                                                             cgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTg                                          <  2:1014304/150‑1 (MQ=255)
                                                                                                             cgcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATAAACGAATCATTCCGATACAACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGTGCGCACGATTTTg                                          >  2:1178440/1‑150 (MQ=255)
                                                                                                              gcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGc                                         <  2:1180921/150‑1 (MQ=255)
                                                                                                              gcaaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGCGCCCGGGCACACgaaaggag                                         <  1:1165609/150‑8 (MQ=255)
                                                                                                                aaccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCAt                                       >  1:1107877/1‑150 (MQ=255)
                                                                                                                 accTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATc                                      >  2:1863260/1‑150 (MQ=255)
                                                                                                                  ccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTTCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATcc                                     >  1:1142415/1‑150 (MQ=255)
                                                                                                                  ccTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATcc                                     <  1:2105213/150‑1 (MQ=255)
                                                                                                                               aTTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGTATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAAcc                        <  2:2026434/150‑1 (MQ=255)
                                                                                                                               aTTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACAGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAAcc                        >  2:2075520/1‑150 (MQ=255)
                                                                                                                                tttCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCCCAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCt                       >  1:852420/1‑150 (MQ=255)
                                                                                                                                  tctcTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCttt                     >  1:1541254/1‑150 (MQ=255)
                                                                                                                                   ctctTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGTAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTc                    >  2:37314/1‑150 (MQ=255)
                                                                                                                                    tctTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTcc                   <  2:347878/150‑1 (MQ=255)
                                                                                                                                     tgaaatgggCACGGAAATCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCt                  >  1:780438/7‑150 (MQ=255)
                                                                                                                                     tgaaatgggCACGGAAATCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCt                  >  1:22024/7‑150 (MQ=255)
                                                                                                                                          ggggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGGGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCgg             <  1:1374739/150‑1 (MQ=255)
                                                                                                                                          ggggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTGATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCgg             <  1:1294960/150‑1 (MQ=255)
                                                                                                                                           gggCACGGATTTCCGTGCCCATTTCACTAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCTCTGTCGTTAACCTTTCCTGCCggg            <  2:863800/150‑1 (MQ=255)
                                                                                                                                           gggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCggg            <  2:1229616/150‑1 (MQ=255)
                                                                                                                                            ggCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGa           <  2:1945118/150‑1 (MQ=255)
                                                                                                                                             gCACGGATTTCCGTGCCCCTTTCACAAGGTGGCTGTTATGCAGGAATACACGCATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCCCGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGaa          >  2:104490/1‑150 (MQ=255)
                                                                                                                                                     ttCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTCCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCCGGCAGGAAAggttaac  <  1:993450/150‑7 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCACGCTTTTGCATCCGCTGTCGTTAACCTTTCCTGCCGGGAAAGTGACCG  >  NC_000913/169600‑169898

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