Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,898,305 A→G intergenic (+9/‑122) sdaA → / → pdeD L‑serine deaminase I/putative c‑di‑GMP phosphodiesterase PdeD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,898,3050AG99.7% 1521.3 / ‑5.0 359intergenic (+9/‑122)sdaA/pdeDL‑serine deaminase I/putative c‑di‑GMP phosphodiesterase PdeD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (185/173);  minor base T (1/0);  total (186/173)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 359 total aligned reads displayed.

AAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCG  >  NC_000913/1898208‑1898402
                                                                                                 |                                                                                                 
aaaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTc                                                                                                <  2:3475552/101‑1 (MQ=255)
aaaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTc                                                                                                <  1:3475899/101‑1 (MQ=255)
 aaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCg                                                                                               <  1:3475647/101‑1 (MQ=255)
   cGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGcc                                                                                             >  1:3476144/1‑101 (MQ=255)
    gATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGccc                                                                                            >  2:3476147/1‑101 (MQ=255)
    gATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGccc                                                                                            <  2:3475587/101‑1 (MQ=255)
     aTGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCa                                                                                           >  2:3476152/1‑101 (MQ=255)
     aTGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCa                                                                                           >  1:3476151/1‑101 (MQ=255)
       gTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATc                                                                                         <  2:3475706/101‑1 (MQ=255)
        tACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCt                                                                                        >  1:3476156/1‑101 (MQ=255)
        tACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCt                                                                                        <  1:3475237/101‑1 (MQ=255)
          cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGc                                                                                      >  2:3476160/1‑101 (MQ=255)
          cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGc                                                                                      <  1:3475520/101‑1 (MQ=255)
          cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGc                                                                                      >  1:3476163/1‑101 (MQ=255)
          cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGc                                                                                      <  2:3475810/101‑1 (MQ=255)
           gAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCa                                                                                     <  2:3475903/101‑1 (MQ=255)
            aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCaa                                                                                    >  2:3476168/1‑101 (MQ=255)
            aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCaa                                                                                    >  1:3476169/1‑101 (MQ=255)
              aCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACg                                                                                  >  2:3476177/1‑101 (MQ=255)
              aCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACg                                                                                  >  1:3476175/1‑101 (MQ=255)
                 ggTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCa                                                                                     >  1:3476184/1‑95 (MQ=255)
                 ggTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGAt                                                                               <  1:3475636/101‑1 (MQ=255)
                 ggTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGAt                                                                               >  2:3476183/1‑101 (MQ=255)
                   tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATgg                                                                             <  1:3475486/101‑1 (MQ=255)
                      ggACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCg                                                                          <  2:3475417/101‑1 (MQ=255)
                       gACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGa                                                                         >  1:3476190/1‑101 (MQ=255)
                        aCATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGaa                                                                        >  2:3476191/1‑101 (MQ=255)
                            gAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTta                                                                    <  2:3475734/101‑1 (MQ=255)
                            gAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTta                                                                    <  1:3475740/101‑1 (MQ=255)
                            gAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTta                                                                    >  1:3476198/1‑101 (MQ=255)
                              aCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTtata                                                                  <  1:3475603/101‑1 (MQ=255)
                                gCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATAcc                                                                >  2:3476201/1‑101 (MQ=255)
                                gCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATAcc                                                                >  1:3476202/1‑101 (MQ=255)
                                 ccAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATAccc                                                               <  2:3475980/101‑1 (MQ=255)
                                  cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCg                                                              >  2:3476205/1‑101 (MQ=255)
                                   aaaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGc                                                             >  2:3476214/1‑101 (MQ=255)
                                   aaaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGc                                                             >  1:3476211/1‑101 (MQ=255)
                                    aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCt                                                            <  2:3475663/101‑1 (MQ=255)
                                    aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCt                                                            >  2:3476217/1‑101 (MQ=255)
                                      tACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATg                                                                              <  2:3475470/81‑1 (MQ=255)
                                      tACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCttt                                                          >  1:3476220/1‑101 (MQ=255)
                                       aCCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTtc                                                         >  1:3476225/1‑101 (MQ=255)
                                         cgcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCg                                                                                               <  2:3475568/61‑1 (MQ=255)
                                         cgcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTtctc                                                       <  2:3475757/101‑1 (MQ=255)
                                         cgcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTtctc                                                       <  2:3475609/101‑1 (MQ=255)
                                            gAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTc                                                    >  2:3476228/1‑101 (MQ=255)
                                            gAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTc                                                    <  2:3475528/101‑1 (MQ=255)
                                            gAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTc                                                    >  2:3476233/1‑101 (MQ=255)
                                            gAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTc                                                    >  1:3476235/1‑101 (MQ=255)
                                             aaaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTct                                                   <  2:3475431/101‑1 (MQ=255)
                                             aaaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTct                                                   >  1:3476238/1‑101 (MQ=255)
                                             aaaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTct                                                   <  1:3475542/101‑1 (MQ=255)
                                             aaaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTct                                                   <  1:3475887/101‑1 (MQ=255)
                                              aaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTct                                                   >  1:3476242/1‑100 (MQ=255)
                                               aCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctgc                                                 <  1:3475681/101‑1 (MQ=255)
                                                ccTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctgct                                                <  2:3476022/101‑1 (MQ=255)
                                                 cTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctgctg                                               >  1:3476247/1‑101 (MQ=255)
                                                  tCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGACTTTATACCCGCTTTCTCTTCTGCTGt                                              <  1:3475347/101‑1 (MQ=255)
                                                  tCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGt                                              >  1:3476250/1‑101 (MQ=255)
                                                     cgcgGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAat                                           <  2:3475320/101‑1 (MQ=255)
                                                      gcgGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAata                                          >  2:3476255/1‑101 (MQ=255)
                                                      gcgGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAata                                          >  1:3476256/1‑101 (MQ=255)
                                                        ggtggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATAtt                                        <  1:3475723/101‑1 (MQ=255)
                                                        ggtggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATAtt                                        >  2:3476263/1‑101 (MQ=255)
                                                        ggtggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATAtt                                        >  1:3476262/1‑101 (MQ=255)
                                                         gtggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTc                                       >  2:3476266/1‑101 (MQ=255)
                                                           ggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTccc                                     <  2:3475594/101‑1 (MQ=255)
                                                            gtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctgctg                                               <  2:3475635/90‑1 (MQ=255)
                                                             tCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTACTATTCCCCa                                   <  1:3476035/101‑1 (MQ=255)
                                                             tCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCa                                   <  2:3475726/101‑1 (MQ=255)
                                                             tCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCa                                   <  1:3475586/101‑1 (MQ=255)
                                                               tGGCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACt                                 >  1:3476271/1‑101 (MQ=255)
                                                                ggCAATCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTa                                <  1:3475819/101‑1 (MQ=255)
                                                                   aaTCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTacac                             >  2:3476277/1‑101 (MQ=255)
                                                                   aaTCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTacac                             >  1:3476276/1‑101 (MQ=255)
                                                                     tCAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACtt                           <  2:3475773/101‑1 (MQ=255)
                                                                      cAAAGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCAc                                  >  2:3476284/1‑93 (MQ=255)
                                                                       aaaGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTcc                         <  1:3475615/101‑1 (MQ=255)
                                                                       aaaGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTcc                         >  1:3476287/1‑101 (MQ=255)
                                                                        aaGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCa                        >  2:3476288/1‑101 (MQ=255)
                                                                         aGTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCAc                       <  1:3475751/101‑1 (MQ=255)
                                                                          gTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACt                      >  1:3476294/1‑101 (MQ=255)
                                                                          gTCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACt                      <  2:3475836/101‑1 (MQ=255)
                                                                           tCCAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTg                     <  2:3475815/101‑1 (MQ=255)
                                                                             cAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGtt                   >  1:3476300/1‑101 (MQ=255)
                                                                             cAGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGtt                   >  2:3476298/1‑101 (MQ=255)
                                                                              aGTGTGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTg                  <  1:3475776/101‑1 (MQ=255)
                                                                               gtgtGACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGc                 >  2:3476303/1‑101 (MQ=255)
                                                                                  tgACTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTc              >  1:3476307/1‑101 (MQ=255)
                                                                                    aCTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAg            <  2:3475867/101‑1 (MQ=255)
                                                                                    aCTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAg            >  1:3476310/1‑101 (MQ=255)
                                                                                    aCTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAg            <  2:3475511/101‑1 (MQ=255)
                                                                                    aCTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAg            <  1:3475547/101‑1 (MQ=255)
                                                                                     cTAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAgg           <  1:3475961/101‑1 (MQ=255)
                                                                                      tAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGc          >  1:3476315/1‑101 (MQ=255)
                                                                                      tAATACTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGc          <  2:3475852/101‑1 (MQ=255)
                                                                                         tACTTTTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGtt       <  1:3475848/101‑1 (MQ=255)
                                                                                          aCTTCTTGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGttt      <  1:3475765/101‑1 (MQ=255)
                                                                                           cttcttGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTg     >  2:3476317/1‑101 (MQ=255)
                                                                                              cttGCTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCg  >  1:3476323/1‑101 (MQ=255)
                                                                                                 |                                                                                                 
AAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCG  >  NC_000913/1898208‑1898402

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: