Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,042,314 A→G 8.4% intergenic (+21/+175) ygfZ → / ← yqfA iron‑sulfur cluster repair protein, plumbagin resistance/hemolysin III family HylIII inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,042,3140AG 8.4% 234.3 / 11.5 86intergenic (+21/+175)ygfZ/yqfAiron‑sulfur cluster repair protein, plumbagin resistance/hemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base G (2/5);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

CTGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCT  >  NC_000913/3042165‑3042463
                                                                                                                                                     |                                                                                                                                                     
cTGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCa                                                                                                                                                       >  1:1071997/1‑150 (MQ=255)
cTGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCa                                                                                                                                                       <  2:953776/150‑1 (MQ=255)
cTGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCa                                                                                                                                                       <  1:157603/150‑1 (MQ=255)
cTGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCAGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCa                                                                                                                                                       >  2:532679/1‑150 (MQ=255)
 tGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGTGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCAt                                                                                                                                                      <  2:948068/150‑1 (MQ=255)
  ggAAGATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCAtt                                                                                                                                                     >  1:468355/1‑150 (MQ=255)
  ggAAGATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCAtt                                                                                                                                                     >  1:32945/1‑150 (MQ=255)
      gATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATc                                                                                                                                                 >  1:1159238/1‑150 (MQ=255)
          gTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAGCCGGATAGCATCTTCCGCGTACGAGACGATGCGAAAACATTGCAAATCGAGACGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGTCATTTATCGcaa                                                                                                                                             >  2:936152/1‑148 (MQ=255)
          gTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCt                                                                                                                                             >  2:1137095/1‑150 (MQ=255)
            cAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGa                                                                                                                                           >  2:738287/1‑150 (MQ=255)
            cAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGa                                                                                                                                           <  1:637494/150‑1 (MQ=255)
            cAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGa                                                                                                                                           <  1:994202/150‑1 (MQ=255)
              ggTCGTGGTACAGGTTGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATg                                                                                                                                         >  2:4410/1‑150 (MQ=255)
              ggTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCGATAAATGCCTGATg                                                                                                                                         <  1:930215/150‑1 (MQ=255)
              ggTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATg                                                                                                                                         >  1:716345/1‑150 (MQ=255)
              ggTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATg                                                                                                                                         >  1:548745/1‑150 (MQ=255)
              gctCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATg                                                                                                                                         <  2:88093/148‑1 (MQ=255)
               gtgggggTACAGGGCGGCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGc                                                                                                                                        <  2:923474/145‑1 (MQ=255)
               gTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCGATAAATGCCTGATGc                                                                                                                                        <  1:1305378/150‑1 (MQ=255)
               gTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATAAATTGCATATCGAGCCGCTGCCATATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGc                                                                                                                                        <  2:850025/150‑1 (MQ=255)
                tCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTTAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCg                                                                                                                                       >  2:697899/1‑150 (MQ=255)
                  gTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCGATAAATGCCTGATGCGat                                                                                                                                     <  1:1050483/150‑2 (MQ=255)
                  gTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCGATAAATGCCTGATGCGat                                                                                                                                     <  1:537304/150‑2 (MQ=255)
                  gTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCGATAAATGCCTGATGCGat                                                                                                                                     <  1:437464/150‑2 (MQ=255)
                   tGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACg                                                                                                                                    <  2:145467/150‑1 (MQ=255)
                          ggTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGgcgc                                                                                                                             >  2:156506/1‑150 (MQ=255)
                           gtcgtcATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGcgcg                                                                                                                            >  1:1254934/1‑150 (MQ=255)
                            tcgtcATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGTATTTATATAATGTTACGTATTCGATCGATTTGTAAATGTCTTTTTCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCgtgt                                                                                                                           <  2:782367/150‑1 (MQ=255)
                              gtcATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCt                                                                                                                         >  2:1113947/1‑150 (MQ=255)
                                   gAATAACGATATGGAACCGGATAGCATCATCCGCGTACGCGACGATGCGAATACAATGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCCCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                                                    <  2:71305/150‑1 (MQ=255)
                                           atatGGAACCGTATTGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAATAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTggg                                                                                                            <  1:510118/150‑1 (MQ=255)
                                                taaCCGTATAGCTGCTTCCGCGTACGCGACGATGCTAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAAGCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTg                                                                                                       <  2:1211780/149‑1 (MQ=39)
                                                gAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTg                                                                                                       >  1:919186/1‑150 (MQ=255)
                                                   ccGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTTCATTTCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATGTATCGCCTGCTGCGCCGCTGGCGCGTCTTAGCTTGTCTGGGAttggctg                                                                                                    >  2:232682/1‑145 (MQ=255)
                                                    cgtaTAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCTTTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGc                                                                                                   <  1:879826/147‑1 (MQ=255)
                                                    cgcaTAGCAGATTCCGCGTACGCGACGATGCTAATACATTGAATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTTATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGc                                                                                                   <  1:518699/147‑1 (MQ=255)
                                                     ggATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATgac                                                                                                  >  2:651704/1‑148 (MQ=255)
                                                     cgATACCATCTTCCGCGTACCCTACGACGCTAATACATTGCATATCGTCCCGCTGCCGTATTCGATCGAAGAGTAAAACTCTTTATCGCATCAGGCATTTAACGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGcc                                                                                                  <  2:329281/149‑1 (MQ=37)
                                                      tatatCATCTTCCGGGTCCGCGACTATGCTAATACATTGCTTATCGCGCCGGGGCCGAATTCGCTCGATGAGTAAATCTCTTTATCGCATCAGGCTTTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCt                                                                                                 <  2:119896/149‑1 (MQ=255)
                                                      gATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCt                                                                                                 >  1:1133456/1‑150 (MQ=255)
                                                       aTAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCtt                                                                                                <  2:882557/150‑1 (MQ=255)
                                                        tAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCttt                                                                                               >  1:1107507/1‑150 (MQ=255)
                                                         aGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCAAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTa                                                                                              >  2:941978/1‑150 (MQ=255)
                                                          gCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTAt                                                                                             >  2:54781/1‑150 (MQ=255)
                                                          gCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTAt                                                                                             >  1:703127/1‑150 (MQ=255)
                                                           cATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATc                                                                                            <  2:813708/150‑1 (MQ=255)
                                                                tCCGCGTACGCGACGATGCGAATACAATGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCAATTATCGCCTGATGCGACGCTGGCGCGGCTTATCAGGCCTGGGATTTGATGCCTGTATcacatc                                                                                       >  2:630073/1‑148 (MQ=39)
                                                                 ccGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCacaaca                                                                                      <  2:898333/150‑1 (MQ=255)
                                                                        cgcgACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGa                                                                               >  1:363006/1‑150 (MQ=255)
                                                                          cgaagaTGATAATACATTGCATATCGAGACGCTGCCTTATTCGCTCTAAGAGTAAATCTATTATTCGAAACAGGAATTTTTCGCCTGATGCGACGCTGTCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAgg                                                                             <  1:248897/150‑1 (MQ=255)
                                                                              gaTGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCg                                                                         <  1:156506/150‑1 (MQ=255)
                                                                                 gCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAg                                                                      >  1:1116349/1‑150 (MQ=255)
                                                                                  cGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGa                                                                     <  2:74859/150‑1 (MQ=255)
                                                                                     aTACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCg                                                                  >  2:1037584/1‑150 (MQ=255)
                                                                                        cATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGAAGCCCCGGTGGCGCGTCTTATCATGTCTCGGATTAGCTGCCTTTAGCACAACAGCGACGAGGCCCGTAGATCAGAt                                                               >  1:357062/1‑150 (MQ=255)
                                                                                         aTTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAGTCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATa                                                              >  2:1229870/1‑150 (MQ=255)
                                                                                         aTTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATa                                                              >  2:895252/1‑150 (MQ=255)
                                                                                          atGCATATCGTGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTAATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATaa                                                             <  2:715201/149‑1 (MQ=255)
                                                                                            gCATTTCGAGCCGCTGCCGTATTCGCTCGAAGTGTAAATCTCTTTATCGCATCTGGCATTTATCGCCTGATGGGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAgg                                                           <  2:495773/150‑1 (MQ=255)
                                                                                            gCATCTCGAGCCGCTGCCGTATTCTCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAgg                                                           <  2:45644/150‑1 (MQ=255)
                                                                                              atatCGAGCCGCGGCCGTATTCGCTCGAAGCGTAAATCTCTTTATCGCATCAGGCATTAATCGCCTGATGCGACGCTGTCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCg                                                         <  1:235935/150‑1 (MQ=255)
                                                                                                          tGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTAACATGTCTGGGATTAGATGCCTTTATCACAACAGGGACGAGGCCCGAAGATCGGATAAGGAGTTAACGACGAAt                                             >  1:160447/1‑150 (MQ=255)
                                                                                                           gCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATc                                            >  1:779144/1‑150 (MQ=255)
                                                                                                            ccGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATcc                                           >  1:494385/1‑150 (MQ=255)
                                                                                                            ccGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATcc                                           >  1:304539/1‑150 (MQ=255)
                                                                                                             cGTATTCGCTCGAAGAGTAAATCTCTTTATCGCAGCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCg                                          <  1:746463/150‑1 (MQ=255)
                                                                                                                ttgcGCTCGAAGATTAAATCTATATATCCCATACGGCATTTAGCCCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCa                                       <  2:1156047/147‑1 (MQ=255)
                                                                                                                aTTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGCCTTATCCgatc                                       <  2:32945/150‑4 (MQ=39)
                                                                                                                aTTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGCCTTATCCgatc                                       <  1:77586/150‑4 (MQ=39)
                                                                                                                aTTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATGTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGCCTTATCCgatc                                       <  1:935586/150‑4 (MQ=255)
                                                                                                                     cTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTtgt                                  >  2:748163/1‑150 (MQ=255)
                                                                                                                      tCGAAGATTAAATCTCTTTATCGCATCAGGCATTTTTCGCCTTATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTtgtg                                 <  2:1071087/150‑1 (MQ=255)
                                                                                                                      tCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTGGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTtgtg                                 >  2:658496/1‑150 (MQ=255)
                                                                                                                       cGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGc                                <  2:475267/150‑1 (MQ=255)
                                                                                                                       cGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCGTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGc                                >  2:322045/1‑150 (MQ=255)
                                                                                                                         aagataaaaTCTCTTTATCGAATCAGGAATTTATCGCCTTATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTc                              <  2:642581/144‑1 (MQ=255)
                                                                                                                                  tctcTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCATATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTg                     <  2:694281/150‑1 (MQ=255)
                                                                                                                                     ctTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATg                  >  2:592487/1‑150 (MQ=255)
                                                                                                                                            gCATCAGGCGATAAATGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCt           >  1:812033/1‑150 (MQ=38)
                                                                                                                                            gCATCAGGCGATAAATGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCt           >  2:1134156/1‑150 (MQ=38)
                                                                                                                                            gCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCt           <  2:816229/150‑1 (MQ=255)
                                                                                                                                             cATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTAACGCCGAATCCGCCAGTTGTGCTCCGATGCCTGATGCGCCGCTg          >  1:887689/1‑150 (MQ=255)
                                                                                                                                               tCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATCCCTGATTCGACGCTggg        >  1:821580/1‑149 (MQ=255)
                                                                                                                                                ccggCATATATCGCCTGAGGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgcg       <  1:291178/148‑1 (MQ=255)
                                                                                                                                                cAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgcg       <  2:919186/150‑1 (MQ=255)
                                                                                                                                                cAGGCATTACTCCCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgcg       <  2:438041/150‑1 (MQ=255)
                                                                                                                                                  ggCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGcgcg     >  2:1123426/1‑150 (MQ=255)
                                                                                                                                                    cATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTc   >  2:150433/1‑150 (MQ=255)
                                                                                                                                                    cATTTATCGCCTGATGCGACGCTGGCGCGGCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTc   <  1:52995/150‑1 (MQ=255)
                                                                                                                                                     aTTTATCTCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCt  <  2:884333/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGAAGATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCT  >  NC_000913/3042165‑3042463

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