Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 281,599 C→T G129S (GGC→AGC)  yagA ← CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,5990CT95.7% 441.0 / 5.1 140G129S (GGC→AGC) yagACP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (5/1);  new base T (63/71);  total (68/72)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 140 total aligned reads displayed.

TGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGC  >  NC_000913/281502‑281691
                                                                                                 |                                                                                            
  gCAGCGTCCACTGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGt                                                                                          <  2:3265396/100‑1 (MQ=255)
  gCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGt                                                                                          >  1:746109/1‑100 (MQ=255)
   cAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTg                                                                                         >  1:2195356/1‑100 (MQ=255)
   cAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTg                                                                                         <  2:2638968/100‑1 (MQ=255)
   cAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTg                                                                                         <  1:416053/100‑1 (MQ=255)
    aGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGa                                                                                        <  2:1594891/100‑1 (MQ=255)
      cGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAg                                                                                      <  1:3245667/100‑1 (MQ=255)
      cGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAg                                                                                      >  2:3127650/1‑100 (MQ=255)
      cGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAg                                                                                      <  1:415863/100‑1 (MQ=255)
       gTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAgc                                                                                     >  1:1291246/1‑100 (MQ=255)
       gTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAgc                                                                                     >  1:1723169/1‑100 (MQ=255)
       gTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAgc                                                                                     <  1:708009/100‑1 (MQ=255)
         ccACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAgcgc                                                                                   >  1:2543133/1‑100 (MQ=255)
         ccACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAgcgc                                                                                   >  2:1401230/1‑100 (MQ=255)
          cACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAgatcg                                                                                  >  1:2188639/1‑96 (MQ=255)
          cACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGcc                                                                                  >  2:162765/1‑100 (MQ=255)
          cACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGcc                                                                                  >  2:2238320/1‑100 (MQ=255)
           aCCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGccc                                                                                 >  2:356527/1‑100 (MQ=255)
              gCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGgc                                                                              <  2:1459478/100‑1 (MQ=255)
              gCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGgc                                                                              >  1:2993653/1‑100 (MQ=255)
               ccAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGgca                                                                             <  2:2726417/100‑1 (MQ=255)
                cAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGgcag                                                                            <  1:1014825/100‑1 (MQ=255)
                cAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGgcag                                                                            <  1:1633739/100‑1 (MQ=255)
                cAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGgcag                                                                            >  2:678878/1‑100 (MQ=255)
                 aaaaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGgcagc                                                                           >  1:3175748/1‑100 (MQ=255)
                  aaaGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGgcagca                                                                          >  1:1756097/1‑100 (MQ=255)
                     ggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGc                                                                       >  2:2005620/1‑100 (MQ=255)
                     ggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGc                                                                       <  1:1098312/100‑1 (MQ=255)
                     ggAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGc                                                                       <  1:3057768/100‑1 (MQ=255)
                      gAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGcc                                                                      <  2:2954930/100‑1 (MQ=255)
                          gTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATgg                                                                  <  1:1652296/100‑1 (MQ=255)
                             gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCTATGGCgg                                                               <  1:1772879/100‑1 (MQ=255)
                              cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCggg                                                              <  2:594010/100‑1 (MQ=255)
                              cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCggg                                                              <  1:1917214/100‑1 (MQ=255)
                              cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCggg                                                              >  1:2105432/1‑100 (MQ=255)
                              cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCggg                                                              >  2:2544305/1‑100 (MQ=255)
                              cccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCggg                                                              <  1:1356657/100‑1 (MQ=255)
                               ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                             <  1:1068276/100‑1 (MQ=255)
                               ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                             >  2:1602534/1‑100 (MQ=255)
                               ccTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGc                                                             <  2:2051209/100‑1 (MQ=255)
                                cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGcc                                                            >  1:1730061/1‑100 (MQ=255)
                                 ttAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCa                                                           <  2:1409475/100‑1 (MQ=255)
                                    aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCa                                                        <  2:742757/100‑1 (MQ=255)
                                     aaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAg                                                       >  2:1307843/1‑100 (MQ=255)
                                      aTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAgg                                                      <  2:1291242/100‑1 (MQ=255)
                                      aTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAgg                                                      <  2:2984895/100‑1 (MQ=255)
                                        ccATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGtt                                                    >  2:2002142/1‑100 (MQ=255)
                                         cATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTa                                                   >  1:2550701/1‑100 (MQ=255)
                                          aTCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTAt                                                  <  1:3247668/100‑1 (MQ=255)
                                           tCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATg                                                 <  2:1720049/100‑1 (MQ=255)
                                            cTGCCAGAGGCGGTTCGGCGCGTCGTCTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATgg                                                >  2:797170/1‑100 (MQ=255)
                                             tGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGa                                               <  1:2406012/100‑1 (MQ=255)
                                              gCCAGAGGCGGTTCGGCGCGTCGTGTTCGGACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCCGCAGGCCATGGCGGGCCATCAGGATACGGAc                                              >  2:1158265/1‑100 (MQ=255)
                                              gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGAc                                              >  1:1099000/1‑100 (MQ=255)
                                              gCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGAc                                              >  2:1672890/1‑100 (MQ=255)
                                                 agagGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGt                                           >  1:1852798/1‑100 (MQ=255)
                                                   agGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGc                                         <  1:2416207/100‑1 (MQ=255)
                                                    ggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCt                                        >  1:1304350/1‑100 (MQ=255)
                                                       ggTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGaa                                     >  2:3259374/1‑100 (MQ=255)
                                                        gTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAg                                    >  2:1765811/1‑100 (MQ=255)
                                                            ggCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCgg                                >  2:601186/1‑100 (MQ=255)
                                                             gcgcGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGTAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggg                               <  2:2359734/100‑1 (MQ=255)
                                                             gcgcGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggg                               <  1:60742/100‑1 (MQ=255)
                                                              cgcgTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGc                              <  2:970162/100‑1 (MQ=255)
                                                              cgcgTCGTGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGc                              >  1:1990749/1‑100 (MQ=255)
                                                                 gtcgtGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATg                           >  1:2252238/1‑100 (MQ=255)
                                                                   cgtGTTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgt                         <  2:712065/100‑1 (MQ=255)
                                                                    gtgtTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtg                        >  2:2740931/1‑100 (MQ=255)
                                                                     tgtTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGgtgt                       >  1:1879962/1‑100 (MQ=255)
                                                                      gtTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGtgtg                      <  2:1756092/100‑1 (MQ=255)
                                                                      gtTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGtgtg                      >  2:64881/1‑100 (MQ=255)
                                                                      gtTCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGtgtg                      >  1:182901/1‑100 (MQ=255)
                                                                      gtTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGATATGGACGGTGCTGAAGGCGGGCATGGtgtg                      >  2:148615/1‑100 (MQ=255)
                                                                       ttCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGACATGGTGTGc                     <  2:1090687/100‑1 (MQ=255)
                                                                        tCGAACCGGCGCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGCTATGGACGGTGCTGAAGGCGGGCATGGTGTGcc                    <  2:2201888/100‑1 (MQ=255)
                                                                        tCGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGcc                    >  2:2788534/1‑100 (MQ=255)
                                                                         cGAACCGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGccc                   >  2:1719078/1‑100 (MQ=255)
                                                                           aaCCGGGCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCt                 <  1:10419/100‑1 (MQ=255)
                                                                             ccGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTgg               >  1:667485/1‑100 (MQ=255)
                                                                             ccGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTgg               <  2:2993636/100‑1 (MQ=255)
                                                                             ccGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTgg               <  2:10781/100‑1 (MQ=255)
                                                                             ccGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTgg               >  2:365998/1‑100 (MQ=255)
                                                                              cGGCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGt              <  1:1401234/100‑1 (MQ=255)
                                                                               ggCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTc             >  2:2827034/1‑100 (MQ=255)
                                                                               ggCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTc             >  1:2624502/1‑100 (MQ=255)
                                                                               ggCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTc             >  1:2116342/1‑100 (MQ=255)
                                                                               ggCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTc             >  2:2183904/1‑100 (MQ=255)
                                                                                gCCCGTGGCGGGAGTGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTcc            >  2:2078641/1‑100 (MQ=255)
                                                                                gCCCGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTcc            <  2:1865026/100‑1 (MQ=255)
                                                                                  ccGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTc          >  2:1655725/1‑100 (MQ=255)
                                                                                  ccGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTc          <  1:1307634/100‑1 (MQ=255)
                                                                                  ccGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTc          <  1:2002150/100‑1 (MQ=255)
                                                                                   cGTGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCg         <  1:1655729/100‑1 (MQ=255)
                                                                                    gTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGa        >  2:648615/1‑100 (MQ=255)
                                                                                     tGGCGGGAATGCTCGGTGAGGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTGTGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAg       <  2:1507834/100‑1 (MQ=255)
                                                                                     tGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAg       <  2:448344/100‑1 (MQ=255)
                                                                                     tGGCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAg       <  1:410550/100‑1 (MQ=255)
                                                                                      ggCGGGAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGc      >  1:2500926/1‑100 (MQ=255)
                                                                                         gggAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAg   <  1:678879/100‑1 (MQ=255)
                                                                                          ggAATGCTCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAgc  <  1:3070101/100‑1 (MQ=255)
                                                                                                 |                                                                                            
TGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGC  >  NC_000913/281502‑281691

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