Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 733,108 C→T 5.6% L1176L (CTG→TTG)  rhsC → Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913733,1080CT 5.6% 218.6 / 3.0 71L1176L (CTG→TTG) rhsCRhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (31/36);  new base T (2/2);  total (33/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 = 9.59e-01

CGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGG  >  NC_000913/732962‑733253
                                                                                                                                                  |                                                                                                                                                 
cgcAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACctg                                                                                                                                                 >  1:375173/1‑149 (MQ=255)
 gcAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAACTTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACctgc                                                                                                                                                <  5:71845/149‑1 (MQ=255)
 gcAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACctgc                                                                                                                                                >  6:88991/1‑149 (MQ=255)
  cAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACctgct                                                                                                                                               <  5:337618/149‑1 (MQ=255)
   aGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACctgct                                                                                                                                               >  2:342157/1‑148 (MQ=255)
      tGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCtgaat                                                                                                                                           <  8:431151/149‑1 (MQ=255)
        cAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACTGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCtgaatga                                                                                                                                         <  6:141209/149‑1 (MQ=255)
         aaaGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCCGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCtgaatgaa                                                                                                                                        <  5:12458/149‑1 (MQ=255)
         aaaGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACTTGCtgaatg                                                                                                                                          >  1:194567/1‑147 (MQ=255)
         aaaGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACTTGCtgaatg                                                                                                                                          >  1:188234/1‑147 (MQ=255)
         aaaGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCtgaatgaa                                                                                                                                        <  8:449633/149‑1 (MQ=255)
         aaaGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCtgaatgaa                                                                                                                                        <  7:377949/149‑1 (MQ=255)
         aaaGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCtgaatg                                                                                                                                          >  6:132803/1‑147 (MQ=255)
          aaGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTTGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAag                                                                                                                                       <  7:192565/149‑1 (MQ=255)
          aaGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAag                                                                                                                                       <  5:88991/149‑1 (MQ=255)
          aaGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAag                                                                                                                                       <  2:375173/149‑1 (MQ=255)
            gCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAagag                                                                                                                                     <  7:210159/149‑1 (MQ=255)
                gATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACTTGCTGAATGAAGAGAAcc                                                                                                                                 <  2:194567/149‑1 (MQ=255)
                gATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACTTGCTGAATGAAGAGAAcc                                                                                                                                 <  2:188234/149‑1 (MQ=255)
                 aTGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAAccc                                                                                                                                <  5:132803/149‑1 (MQ=255)
                   ggACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACctgct                                                                                                                                               >  6:50629/1‑132 (MQ=255)
                   ggACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACctgct                                                                                                                                               <  5:50629/132‑1 (MQ=255)
                     aCCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCAt                                                                                                                            >  7:308793/1‑149 (MQ=255)
                       ccGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCAt                                                                                                                            >  6:437419/1‑147 (MQ=255)
                           tataCACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAGTGAAGAGAACCCGCATCAGCTg                                                                                                                      <  2:92060/149‑1 (MQ=255)
                           tataCACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTg                                                                                                                      <  1:342157/149‑1 (MQ=255)
                           tataCACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTg                                                                                                                      <  8:236337/149‑1 (MQ=255)
                              acacGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcag                                                                                                                   >  8:243076/1‑149 (MQ=255)
                                acGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcag                                                                                                                   >  7:275708/1‑147 (MQ=255)
                                     ggCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTa                                                                                                            <  3:293367/149‑1 (MQ=255)
                                        gcgAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATcc                                                                                                         <  2:456950/149‑1 (MQ=255)
                                                 tCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCgg                                                                                                >  6:230548/1‑149 (MQ=255)
                                                  cACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCgg                                                                                                <  5:437419/148‑1 (MQ=255)
                                                     cTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCaggg                                                                                              >  8:61930/1‑147 (MQ=255)
                                                        tACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcag                                                                                                                   >  7:72463/1‑123 (MQ=255)
                                                        tACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcag                                                                                                                   <  8:72463/123‑1 (MQ=255)
                                                                  aCCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTat                                                                                     <  1:63240/143‑1 (MQ=255)
                                                                  aCCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTat                                                                                     >  2:63240/1‑143 (MQ=255)
                                                                         cgGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGATTGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCg                                                                        <  1:349421/149‑1 (MQ=255)
                                                                         cgGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCACATCAGCTgc                                                                                                                     <  8:140577/104‑1 (MQ=255)
                                                                         cgGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCACATCAGCTgc                                                                                                                     >  7:140577/1‑104 (MQ=255)
                                                                           gCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGc                                                                      <  7:243076/149‑1 (MQ=255)
                                                                              tGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTg                                                                   >  1:449239/1‑149 (MQ=255)
                                                                                ccGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTa                                                                 >  2:419869/1‑149 (MQ=255)
                                                                                 cGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTAt                                                                <  8:308793/149‑1 (MQ=255)
                                                                                  gCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTAtt                                                               <  2:449239/149‑1 (MQ=255)
                                                                                      gCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACTCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACaa                                                           >  6:92463/1‑149 (MQ=255)
                                                                                        ccTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAAc                                                          >  1:199899/1‑148 (MQ=255)
                                                                                        ccTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATAACAAcc                                                         >  5:125696/1‑149 (MQ=255)
                                                                                               aGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCg                                                  >  5:403928/1‑149 (MQ=255)
                                                                                               aGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCg                                                  >  7:78637/1‑149 (MQ=255)
                                                                                                   cGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATCACAACCGCCACCGCtat                                              <  6:125696/149‑1 (MQ=255)
                                                                                                     gAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCta                                               <  5:277422/146‑1 (MQ=39)
                                                                                                     gAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCta                                               >  6:277422/1‑146 (MQ=39)
                                                                                                      aaGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATgag                                                                               <  8:214586/113‑1 (MQ=255)
                                                                                                      aaGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATgag                                                                               >  7:214586/1‑113 (MQ=255)
                                                                                                         ggACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGAc                                        <  1:419869/149‑1 (MQ=39)
                                                                                                             aacaGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGc                                    >  5:173942/1‑149 (MQ=39)
                                                                                                             aacaGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGc                                    >  8:468068/1‑149 (MQ=39)
                                                                                                                      gTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCACCACCGCTATTATGACCCGCTGCAGGGGc                           <  7:61930/149‑1 (MQ=21)
                                                                                                                         cGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGat                        <  8:275708/149‑1 (MQ=33)
                                                                                                                         cGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGa                         <  7:468068/148‑1 (MQ=32)
                                                                                                                          gCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                       >  8:372688/1‑149 (MQ=34)
                                                                                                                            aGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatata                     >  7:468129/1‑149 (MQ=34)
                                                                                                                              aaTACGCTGAATGGGGTAACCTGCTGACTGCAGCGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAAcc                                                         <  4:278625/111‑1 (MQ=11)
                                                                                                                              aaTACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAAcc                                                         >  3:278625/1‑111 (MQ=32)
                                                                                                                                tACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCAc                 >  8:112015/1‑149 (MQ=34)
                                                                                                                                tACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCAc                 >  8:472741/1‑149 (MQ=34)
                                                                                                                                   gATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCa                  >  3:188346/1‑145 (MQ=34)
                                                                                                                                   gATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCa                  <  4:188346/145‑1 (MQ=34)
                                                                                                                                   gATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACt                >  7:285473/1‑147 (MQ=34)
                                                                                                                                   gATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCa              <  6:173942/149‑1 (MQ=34)
                                                                                                                                    aTGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAg             >  8:234230/1‑149 (MQ=34)
                                                                                                                                           gggTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGa      >  6:247688/1‑149 (MQ=34)
                                                                                                                                           gggTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGa      >  7:386465/1‑149 (MQ=34)
                                                                                                                                             gTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCGTGTATTAGAACCGCCAGCGCTATTATGCCCGGCTGCAGGGGCGATATATCACTCAGGATCCGAtt    <  5:230548/149‑1 (MQ=11)
                                                                                                                                             gTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAtt    <  8:386465/149‑1 (MQ=33)
                                                                                                                                                aCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCACCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTgg  <  2:199899/148‑1 (MQ=17)
                                                                                                                                                  |                                                                                                                                                 
CGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGG  >  NC_000913/732962‑733253

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