Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,420,963 A→G 5.1% T94A (ACC→GCC)  ydaT → Rac prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,420,9630AG 5.1% 243.3 / 2.8 78T94A (ACC→GCC) ydaTRac prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base G (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

AAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCCGACGATTCTGTCATTGTTCATTAAGC  >  NC_000913/1420819‑1421108
                                                                                                                                                |                                                                                                                                                 
aaGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCt                                                                                                                                               >  2:256740/1‑149 (MQ=255)
 aGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTg                                                                                                                                              >  5:365836/1‑149 (MQ=255)
          aCAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGt                                                                                                                                     <  8:52132/149‑1 (MQ=255)
              tCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCAt                                                                                                                                 <  7:235820/149‑1 (MQ=255)
                 gcagAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTAc                                                                                                                              <  2:121420/149‑1 (MQ=255)
                  cagAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACt                                                                                                                             >  3:33226/1‑149 (MQ=255)
                    gAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAAGGCAAAGAGAAAAAATTCAGAATCTGATACCAGCAATTCTGGAAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACACTGGAACGCCGTGCATTACTgc                                                                                                                           >  1:226512/1‑148 (MQ=255)
                    gAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTgg                                                                                                                           <  8:79050/149‑1 (MQ=255)
                     aaTATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCt                                                                                                                                               >  5:312918/1‑128 (MQ=255)
                     aaTATCTATCACCGTAGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCt                                                                                                                                               <  6:312918/128‑1 (MQ=255)
                          ctatcACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGgcggcg                                                                                                                      >  3:97662/1‑148 (MQ=255)
                          ctatcACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGgcggcg                                                                                                                      >  1:206588/1‑148 (MQ=255)
                          ctatcACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGgcggcg                                                                                                                      >  2:196152/1‑148 (MQ=255)
                             tcACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGGCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAgg                                                                                                                  <  1:215316/149‑1 (MQ=255)
                               aCCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAgg                                                                                                                  >  1:120316/1‑147 (MQ=255)
                               aCCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGaa                                                                                                                <  1:126443/149‑1 (MQ=255)
                                ccGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCTTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAg                                                                                                               <  1:283239/149‑1 (MQ=255)
                                 cGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGc                                                                                                              <  6:129412/149‑1 (MQ=255)
                                  gttggttgAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCg                                                                                                             >  2:187101/1‑149 (MQ=255)
                                      tttgAAAGGTGAAACGACAACGCAAAGAGAAAAAATTCAGAAGCGGTTCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCATGAAGCGTTa                                                                                                          <  1:110789/147‑1 (MQ=255)
                                      gttgAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGtt                                                                                                           >  5:189473/1‑147 (MQ=255)
                                      gttgAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTaa                                                                                                         <  3:289903/149‑1 (MQ=255)
                                       ttgAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCCCCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAg                                                                                                        >  6:338849/1‑149 (MQ=255)
                                       ttgAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAg                                                                                                        >  3:269073/1‑149 (MQ=255)
                                         gAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGt                                                                                                                                     <  7:367952/118‑1 (MQ=255)
                                         gAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGt                                                                                                                                     >  8:367952/1‑118 (MQ=255)
                                         gAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTa                                                                                                      >  1:291891/1‑149 (MQ=255)
                                           aaGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACg                                                                                                    >  4:47058/1‑149 (MQ=255)
                                                 aaaCGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAAtt                                                                                              >  3:39592/1‑149 (MQ=255)
                                                  aaCGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTGCTGGCGGCGCAGGAAGCGTTAAGTACGGCAAtt                                                                                              >  3:328555/1‑148 (MQ=255)
                                                        aaCGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATgcgc                                                                                       <  8:276441/149‑1 (MQ=255)
                                                         aCGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCa                                                                                      <  3:80010/149‑1 (MQ=255)
                                                          cGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCat                                                                                     <  7:133470/149‑1 (MQ=255)
                                                          cGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCat                                                                                     <  1:196152/149‑1 (MQ=255)
                                                             aaaGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCatgat                                                                                  >  4:129801/1‑149 (MQ=255)
                                                              aaGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCatgat                                                                                  >  8:17872/1‑148 (MQ=255)
                                                               agagAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATGCCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCatgatga                                                                                >  7:179372/1‑149 (MQ=255)
                                                               agagAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATGCCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCatgatga                                                                                >  7:184909/1‑149 (MQ=255)
                                                               agagAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCatgatga                                                                                <  2:120316/149‑1 (MQ=255)
                                                                gagaAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCTCGAGCTGCGGCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCatgatgat                                                                               <  8:79965/149‑1 (MQ=255)
                                                                   aaaaaaTTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCa                                                                            <  8:177733/149‑1 (MQ=255)
                                                                        aTTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTaa                                                                                                         >  6:224553/1‑115 (MQ=255)
                                                                        aTTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTaa                                                                                                         <  5:224553/115‑1 (MQ=255)
                                                                            aGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGcc                                                                   >  2:88493/1‑149 (MQ=255)
                                                                              aaGCTGATACCAGCAATTCTGGCAATCCTTCCGCTCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGt                                                                 >  2:264108/1‑149 (MQ=255)
                                                                                   gATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTAcc                                                            >  1:339409/1‑149 (MQ=255)
                                                                                   gATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTAcc                                                            >  7:1241/1‑149 (MQ=255)
                                                                                    aTCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATGCCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTAc                                                             <  8:184909/147‑1 (MQ=255)
                                                                                    aTCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATGCCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTAc                                                             <  8:179372/147‑1 (MQ=255)
                                                                                          gCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAgc                                                     <  4:97662/149‑1 (MQ=255)
                                                                                             aTTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCa                                                  <  1:196021/149‑1 (MQ=255)
                                                                                             aTTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCa                                                  >  6:140018/1‑149 (MQ=255)
                                                                                              ttCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGc                                                                             >  3:50490/1‑121 (MQ=255)
                                                                                              ttCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGc                                                                             <  4:50490/121‑1 (MQ=255)
                                                                                                 tGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATgc                                                                                         <  5:113745/106‑1 (MQ=255)
                                                                                                 tGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATgc                                                                                         >  6:113745/1‑106 (MQ=255)
                                                                                                  ggCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTc                                              <  4:269073/148‑1 (MQ=255)
                                                                                                    cAATCCTTCCGCGCGCTCTGCGTCACCGACTCTGCATCTTCGATACCCTGGACCGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGCTGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTc                                              <  6:223782/146‑1 (MQ=255)
                                                                                                     aaTCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAgcg                                          <  2:291891/149‑1 (MQ=255)
                                                                                                     aaTCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAgcg                                          >  8:12710/1‑149 (MQ=255)
                                                                                                     aaTCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAgcg                                          >  4:323752/1‑149 (MQ=255)
                                                                                                            ccGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAgcggcggcg                                    <  3:211179/148‑1 (MQ=255)
                                                                                                            ccGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAgcggcggcg                                    <  4:144233/148‑1 (MQ=255)
                                                                                                                  gAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGttcttc                             <  1:27983/149‑1 (MQ=255)
                                                                                                                       gCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCcgac                         >  7:263741/1‑148 (MQ=255)
                                                                                                                          tccccGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCCGACGAtt                     <  2:142662/146‑1 (MQ=255)
                                                                                                                          tCACCGACTCTGCATCTTCGATACCCTGTAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGTCAATTTATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCCGACTAtt                     <  2:226512/149‑1 (MQ=255)
                                                                                                                          tCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCCGACGAtt                     >  6:370507/1‑149 (MQ=255)
                                                                                                                             ccGACTCTGTATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCCGACGATTCTg                  >  1:15649/1‑149 (MQ=255)
                                                                                                                             ccGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCCGACGATTCTg                  <  4:39592/149‑1 (MQ=255)
                                                                                                                                    tGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCCGACGATTCTGTCATTGt           >  5:388035/1‑149 (MQ=255)
                                                                                                                                    tGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCCGACGATTCTGTCATAGt           >  5:306458/1‑149 (MQ=255)
                                                                                                                                      cATCTTAGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCCGACGATTCTGTCATTGTTc         <  1:264108/149‑1 (MQ=255)
                                                                                                                                         cTTCGATACCCTGGAACGCCGTGGATTACTGGCGGAGCAGGAAGCGTTAAGTACTGCAATTGATGCGCATGATGATGAAGTCCTAGCCGTTTTCCGGAAAGGGCATTTCAGCGGCGTCTTTTCTTCCGACGATTCTGTCTTTGTTCa        <  8:12113/147‑1 (MQ=255)
                                                                                                                                            cGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCGGCGGCGGCGGTTCTTCCGACGATTCTGTCATTGTTCATTaa    <  3:47058/148‑1 (MQ=255)
                                                                                                                                            cGATACACTGGACCGCCGTGCATTACTGGCGGCGCACGAAGCGGTAAGTAAGGCAATTGATGCGCATGATGATGCAGTCCACGCCGGTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCATCCGACGCTTCAGTCATTGGTCAGTagg   >  3:262856/1‑147 (MQ=255)
                                                                                                                                             gATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCAGCGGCGGTTCTTCCGACGATTCTGTCATTGTTCATTAAg   >  4:180280/1‑148 (MQ=255)
                                                                                                                                              aTACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCCGACGATTCTGTCATTGTTCATTAAGc  >  5:340890/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGCGAGCTGCGTCACCGACTCTGCATCTTCGATACCCTGGAACGCCGTGCATTACTGGCGGCGCAGGAAGCGTTAAGTACGGCAATTGATGCGCATGATGATGCAGTCCAAGCCGTTTACCGGAAAGCGCATTTCAGCGGCGGCGGTTCTTCCGACGATTCTGTCATTGTTCATTAAGC  >  NC_000913/1420819‑1421108

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: ‑