Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,042,319 T→A 9.2% intergenic (+26/+170) 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,3190TA 9.2% 197.9 / 10.7 85intergenic (+26/+170)ygfZ/yqfAiron‑sulfur cluster repair protein, plumbagin resistance/hemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base A (2/6);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

GATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTAT  >  NC_000913/3042171‑3042466
                                                                                                                                                    |                                                                                                                                                   
gATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATc                                                                                                                                                    >  1:1159238/1‑150 (MQ=255)
    gTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAGCCGGATAGCATCTTCCGCGTACGAGACGATGCGAAAACATTGCAAATCGAGACGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGTCATTTATCGcaa                                                                                                                                                >  2:936152/1‑148 (MQ=255)
    gTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCt                                                                                                                                                >  2:1137095/1‑150 (MQ=255)
      cAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGa                                                                                                                                              <  1:994202/150‑1 (MQ=255)
      cAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGa                                                                                                                                              >  2:738287/1‑150 (MQ=255)
      cAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGa                                                                                                                                              <  1:637494/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:537304/150‑2 (MQ=255)
            gTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCGATAAATGCCTGATGCGat                                                                                                                                        <  1:1050483/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                                          <  1:77586/150‑4 (MQ=39)
                                                                                                          aTTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGCCTTATCCgatc                                          <  2:32945/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)
                                                                                                                                                  tATCGCCGGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTAt  <  1:1097750/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATGGTCAGGTCGTGGTACAGGTCGTCATGAATAACGATATGGAACCGGATAGCATCTTCCGCGTACGCGACGATGCGAATACATTGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTAT  >  NC_000913/3042171‑3042466

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: