Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,234,024 G→A 5.3% L41L (CTG→TTG)  higA ← antitoxinof the HigB‑HigA toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,234,0240GA 5.3% 246.4 / 2.6 75L41L (CTG→TTG) higAantitoxinof the HigB‑HigA toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (34/37);  new base A (2/2);  total (36/39)
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.98e-01

ATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTCCCCTTAGTACGATGAACAGCG  >  NC_000913/3233878‑3234171
                                                                                                                                                  |                                                                                                                                                   
aTCCATAAGGGTACGAATCCCGGCTATACCGCCAGGCATGGCTTGAGCCATGGCTGTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGagcag                                                                                                                                                     >  4:124222/1‑147 (MQ=255)
   cATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAga                                                                                                                                                 <  7:22010/148‑1 (MQ=255)
    aTAAGGGTACGAATCACGGCTATACCGCCAGTCATGTCTTGAGACCTGGCATTAAATTCCTCAAATTCGGTAGATGATTCTTCCCTCGAGGGTATTTTGGCACACAACAGATCAAGCAAGGGGTTTTCAGGATCGTGGAGCAGCAgat                                                                                                                                                >  7:274358/1‑148 (MQ=255)
         ggTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTa                                                                                                                                          >  3:26262/1‑149 (MQ=255)
          gTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTAc                                                                                                                                         <  4:26262/149‑1 (MQ=255)
             cGAATCACGGCTCTACCGCCAGGCATGGCTTGAGCCATGGCAGTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGGTGAGCAGCAGATGATCTACCa                                                                                                                                       <  2:261831/148‑1 (MQ=255)
              gAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGTACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGt                                                                                                                                     >  2:397203/1‑149 (MQ=255)
                   aCGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAg                                                                                                                                >  7:115915/1‑149 (MQ=255)
                        tataCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCt                                                                                                                           >  2:60126/1‑149 (MQ=255)
                          taCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTgg                                                                                                                         <  5:359084/149‑1 (MQ=255)
                          taCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTgg                                                                                                                         <  4:244178/149‑1 (MQ=255)
                           aCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTggg                                                                                                                        >  7:291046/1‑149 (MQ=255)
                            ccGCCAGTCTTGTCTTGAGCCATGGCATTAAATTCATCAAATGCGGGCGCTGATTCTTCCCACTCGGTTATTTTGGCACACACCATATCCAGCAATGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGgt                                                                                                                       <  8:274358/149‑1 (MQ=255)
                            ccGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCTCCAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTGCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGgt                                                                                                                       <  5:239960/149‑1 (MQ=255)
                            ccGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTgg                                                                                                                         <  4:249316/147‑1 (MQ=255)
                             cGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGgtg                                                                                                                      >  3:182606/1‑149 (MQ=255)
                              gCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGgtg                                                                                                                      >  6:288731/1‑148 (MQ=255)
                              gCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGgtg                                                                                                                      >  3:84947/1‑148 (MQ=255)
                               ccAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTa                                                                                                                    <  5:156799/149‑1 (MQ=255)
                                   gCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGt                                                                                                                >  2:293066/1‑149 (MQ=255)
                                      tGGCTTGAGCCATGGCATTAAATTCCGCAAGTTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTcc                                                                                                             >  1:286612/1‑149 (MQ=255)
                                        gCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGTCACACACCAGATCCAGCAAGGGGTTTTCAGGAGCGTTTAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTTTATTGTTCCTc                                                                                                           <  2:221782/149‑1 (MQ=255)
                                            gAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGCGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTACTCGTTc                                                                                                       >  2:101256/1‑149 (MQ=255)
                                            gAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGtt                                                                                                        <  8:115915/148‑1 (MQ=255)
                                             aGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCt                                                                                                      <  4:84947/149‑1 (MQ=255)
                                                cATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGaa                                                                                                   <  5:187640/149‑1 (MQ=255)
                                                cATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGaa                                                                                                   <  3:338639/149‑1 (MQ=255)
                                                  tGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATa                                                                                                 >  7:386992/1‑149 (MQ=255)
                                                      aTTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGGAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCg                                                                                             >  4:155488/1‑149 (MQ=255)
                                                      aTTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCg                                                                                             >  3:214025/1‑149 (MQ=255)
                                                       ttAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACCCACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGc                                                                                            >  3:266781/1‑149 (MQ=255)
                                                       ttAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGc                                                                                            <  1:83324/149‑1 (MQ=255)
                                                         aaaTTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATGCAGCAAGGGGTTTTCAGGTTCGTTGAGGAGCAGATGATCTACCTGTTCCAGGGCCTGGGTGTATTGTTGGTCGTTCTGAATGCCCGcct                                                                                          >  4:197920/1‑148 (MQ=255)
                                                           aatCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGAGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTg                                                                                                     >  7:145123/3‑136 (MQ=255)
                                                           aTTCTGCAAATTCGGGCGCTGATTCTTCCCACGGGGTTATTTTGTCACAGTCCAGATCCAGCAAGGGGTTTTCAGGAGCGTTGAGCTGCAGATGTTCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTg                                                                                                     <  8:145123/136‑1 (MQ=255)
                                                               cGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAg                                                                                    >  3:259563/1‑149 (MQ=255)
                                                               cGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAg                                                                                    <  4:259563/149‑1 (MQ=255)
                                                                gCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAgg                                                                                   >  6:163052/1‑149 (MQ=255)
                                                                 cAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGt                                                                                  <  8:257653/149‑1 (MQ=255)
                                                                       cGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCcac                                                                             <  5:288731/148‑1 (MQ=255)
                                                                        gggCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCAcag                                                                           <  7:261069/149‑1 (MQ=255)
                                                                          gcgcTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCAcag                                                                           >  4:11675/1‑147 (MQ=255)
                                                                            gcTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGt                                                                       <  3:124222/149‑1 (MQ=255)
                                                                            gcTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGt                                                                       >  4:195042/1‑149 (MQ=255)
                                                                               gATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAg                                                                    >  1:40064/1‑149 (MQ=255)
                                                                                aTTCTTCCCACGCGGTTATTTGGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGc                                                                   <  7:368632/149‑1 (MQ=255)
                                                                                   cttcCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCttt                                                                >  8:391084/1‑149 (MQ=255)
                                                                                    ttcCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCtttt                                                               >  3:163550/1‑149 (MQ=255)
                                                                                      cccACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCtttt                                                               <  8:192941/147‑1 (MQ=255)
                                                                                       ccACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGAGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTtctc                                                            >  2:239886/1‑149 (MQ=255)
                                                                                         aCGCGGTTATTTTGGCACACACCAGATCCAACAAGGGGTTTTCAGGATCGGTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGACAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCt                                                          >  5:311298/1‑149 (MQ=255)
                                                                                           gcgGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGc                                                        <  3:195042/149‑1 (MQ=255)
                                                                                           gcgGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGc                                                        <  3:11675/149‑1 (MQ=255)
                                                                                             ggTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCtt                                                      <  1:293066/149‑1 (MQ=255)
                                                                                               ttATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTg                                                     >  1:391599/1‑148 (MQ=255)
                                                                                                  ttttGGCACCCACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCaa                                                  <  4:266781/148‑1 (MQ=255)
                                                                                                       gCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTc                                            <  1:397203/149‑1 (MQ=255)
                                                                                                         acacacCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCgg                                          >  3:177859/1‑149 (MQ=255)
                                                                                                         acacacCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAACAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTctg                                          >  4:809/1‑147 (MQ=255)
                                                                                                         acacacCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAACAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCgg                                          >  2:296448/1‑149 (MQ=255)
                                                                                                                  aTCCTGCAAGGGGTTTTCAGGATGGTTGAGCAGCTGATGATCTAGCAGTTCCAGCGCCTGGGTGTTTTGTTGCTCGTGCTTAATACCCGCCAGAATAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGgcaatcgcaa                                 <  4:177859/149‑1 (MQ=255)
                                                                                                                  aTCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGgcaatcgcaa                                 <  2:40064/149‑1 (MQ=255)
                                                                                                                    ccAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGgcaatcgcaat                                >  1:140946/1‑148 (MQ=255)
                                                                                                                     cAGCAAGGGGTTTTCAGGCTCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGgcaatcgcaat                                <  2:286612/147‑1 (MQ=255)
                                                                                                                      aGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCAt                             >  1:397295/1‑149 (MQ=255)
                                                                                                                       gCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCAtt                            >  1:88905/1‑149 (MQ=255)
                                                                                                                        cAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCAttt                           <  7:266613/149‑1 (MQ=255)
                                                                                                                             ggTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTcc                      <  1:21457/149‑1 (MQ=255)
                                                                                                                              gTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTccc                     <  4:182606/149‑1 (MQ=255)
                                                                                                                              gTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTccc                     <  2:88905/149‑1 (MQ=255)
                                                                                                                              gTTTTCAGGATCGTTGAGCAACAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTccc                     <  3:809/149‑1 (MQ=255)
                                                                                                                              gTTTTCAGGATCGTTGAGCAACAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTccc                     <  1:296448/149‑1 (MQ=255)
                                                                                                                              gTTTTAAGGATCGGTGCGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGGTCCTCGTTCGGCATACCCGCCAGATAAGGTGCCACAGCAGTTAGCTGTTCTCCGGCTTGCAAGATGTCGGCAATCGCAATCATGTTTTccc                     >  4:203497/1‑149 (MQ=255)
                                                                                                                               ttttCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTcccc                    >  7:319627/1‑149 (MQ=255)
                                                                                                                               ttttCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTcccc                    <  8:291046/149‑1 (MQ=255)
                                                                                                                                        tCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTCCCCTTAGTACGa           <  8:386992/149‑1 (MQ=255)
                                                                                                                                                cagcagATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTCCCCTTAGTACGATGAACAGc   <  4:163550/149‑1 (MQ=255)
                                                                                                                                                cagcagATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTCCCCTTAGTACGATGAACAGc   >  8:71002/1‑149 (MQ=255)
                                                                                                                                                 agctgATGATCTACCAGGTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGAGCCACAGCAGTTAGCTTTTCTCCTGCTTGCTAGATGACGGCAATCGCAATGATTTTTTCCCCTTAGTACGAGGAACAgag  >  6:333234/1‑147 (MQ=255)
                                                                                                                                                 agcagATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTCCCCTTAGTACGAt          <  5:289314/141‑1 (MQ=255)
                                                                                                                                                 agcagATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTACTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTCCCCTTAGTACGAt          >  6:289314/1‑141 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTCTGAATACCCGCCAGAAAAGGTGCCACAGCAGTTAGCTTTTCTCCTGCTTGCAAGATGTCGGCAATCGCAATCATTTTTTCCCCTTAGTACGATGAACAGCG  >  NC_000913/3233878‑3234171

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