Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,594 G→C 5.4% W98C (TGG→TGC yafO → mRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913252,5940GC 5.4% 449.5 / 11.5 138W98C (TGG→TGCyafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (79/50);  new base C (4/3);  total (84/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 138 total aligned reads displayed.

CCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCT  >  NC_000913/252446‑252741
                                                                                                                                                    |                                                                                                                                                   
ccTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGt                                                                                                                                                    <  1:1670581/150‑1 (MQ=255)
ccTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTAAGGGGGCTTTTGATGAGCAAGCATGGt                                                                                                                                                    >  2:270189/1‑150 (MQ=255)
   ggCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGc                                                                                                                                                 >  1:1502165/1‑150 (MQ=255)
      cATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATCGTTGCTCa                                                                                                                                              >  1:1224182/1‑150 (MQ=255)
       aTTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCAATTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAt                                                                                                                                             <  1:450462/150‑1 (MQ=255)
        ttAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAtt                                                                                                                                            <  2:672686/150‑1 (MQ=255)
         tAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACCGTTGCGCCACTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTg                                                                                                                                           <  2:1056922/150‑1 (MQ=255)
           aTCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGACTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                         <  1:1344744/150‑1 (MQ=255)
           aTCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTTGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                         >  1:1549255/1‑150 (MQ=255)
           aTCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                         >  2:449603/1‑150 (MQ=255)
           aTCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                         >  1:1463740/1‑150 (MQ=255)
            tCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCATTTCAGCAGAACGCATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCCTTGcct                                                                                                                                        >  2:269379/1‑149 (MQ=255)
             cAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAt                                                                                                                                       >  2:724749/1‑150 (MQ=255)
             aaaaTTTGAGCGAGTAGCTCATATTCATCTGGCAAATGAGAATATTCCAATTCCGCCACAGTTGCGCAAATTCAGCAGAACGAATGACGAATCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaacg                                                                                                                                       <  1:145093/149‑7 (MQ=255)
                        gAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTg                                                                                                                            >  2:1626636/1‑150 (MQ=255)
                            tGCTGATATTCATCTGGCAAATGAGAATAATCCATTTCCGCAACAGTTGCGCCAATTCAGCAGAACGTATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAcc                                                                                                                        <  2:1145652/150‑1 (MQ=255)
                              cTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGTCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTc                                                                                                                      >  2:1451988/1‑150 (MQ=255)
                              cTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTc                                                                                                                      <  1:571237/150‑1 (MQ=255)
                              cTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAATAGCGCCAAATCAGCAGAACGCAACACGCCGCGCATTAGCTATAAAGTCCGGGGGCGTTAGATCAGCCAGCAAGGTTACTCATTGCCATTCAGCACCCTGACCCTc                                                                                                                      >  2:1210593/1‑150 (MQ=255)
                               tCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGCCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGGTCATTGCCATTCTGAACCGGGAACCTCa                                                                                                                     >  2:311180/1‑150 (MQ=255)
                               tCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                                                                     >  2:365950/1‑150 (MQ=255)
                                cATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGACCCTCAt                                                                                                                    >  1:605165/1‑150 (MQ=255)
                                 atatTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATa                                                                                                                   <  2:504039/150‑1 (MQ=255)
                                      cATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTg                                                                                                              >  2:865739/1‑150 (MQ=255)
                                      cATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTg                                                                                                              >  1:387409/1‑150 (MQ=255)
                                      cATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTCCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGTGGGCGTTTGATGAGCAAGCAAGGTTGCTCATAGCGATTCTGAAACCTGAACCTCATAAAcag                                                                                                              >  1:276924/1‑148 (MQ=255)
                                      cATCTGGAAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGAGCATTTGGTATATTGTCAGCGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGACATTATGATAACTCAACCTCATAAACtt                                                                                                              >  2:518277/1‑149 (MQ=255)
                                       aTCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTgg                                                                                                             >  1:1556070/1‑150 (MQ=255)
                                        tCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                                                                            <  2:1212333/150‑1 (MQ=255)
                                         cTGGCCAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCACCAGAACGACCGACGCAGCCCACTTGGAATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGGTGCTCATTGCCATTCTGAAACCTGAACCTCATAACCTGGCt                                                                                                           <  2:262656/150‑1 (MQ=255)
                                          tGGCAAAAGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATACTGAAACCTGAACCTCATAAACTGGCTc                                                                                                          >  1:1674260/1‑150 (MQ=255)
                                           ggCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGATCCTCATAAACTGGCTCg                                                                                                         >  2:37883/1‑150 (MQ=255)
                                           ggCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTAAGCAGAACGAATGACGAAGCTCATTTGGTAGATTGTCAGGGGGCGTTTGATGCGCAAGTATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                                                                         >  1:842866/1‑150 (MQ=255)
                                             cAAATGAGAATAATCCATTTGCGCGACAGTTGCGGCAATTCAGCAGAACGAATGACGATGCGCTTTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                                                                       <  2:594493/150‑1 (MQ=255)
                                             cAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                                                                       >  2:839720/1‑150 (MQ=255)
                                              aaaTGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgaga                                                                                                      >  1:458966/1‑150 (MQ=255)
                                              aaaTGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGATCATTGCCATTCTGAAACCTGAACATCATAAAATGGATCgata                                                                                                      >  2:1287913/1‑148 (MQ=255)
                                               aaTGAGAATAATACATTTCCGCCACAGGTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGTCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGAt                                                                                                     <  1:270189/150‑1 (MQ=255)
                                                 tGAGAATAATCCATTTCCGACACAGTTGCGCCATTTCAGCAGAACGAATGACGAATCGCATTTGGTATATTGTAAGGGGGCGTTTGATGACCAACCATGGTTGCTCATTGACAATCTCAAACCTGAACATAATAAACTGGCTCGAGATaa                                                                                                   >  2:237954/1‑150 (MQ=255)
                                                        aatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaa                                                                                            >  2:480758/1‑150 (MQ=255)
                                                          tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                                                                         >  1:575897/1‑135 (MQ=255)
                                                          tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                                                                         <  2:575897/135‑1 (MQ=255)
                                                          tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAAt                                                                                          >  2:815087/1‑150 (MQ=255)
                                                            cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGc                                                                                        >  1:771691/1‑150 (MQ=255)
                                                            cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGc                                                                                        >  1:685391/1‑150 (MQ=255)
                                                               ttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGAATATTGTCATGGGGCGGTTGATGAGCAAGCATGTTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATa                                                                                     <  1:237954/150‑1 (MQ=255)
                                                                tCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaa                                                                                    <  2:1173153/150‑1 (MQ=255)
                                                                    ccACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACAGGCTCGAGATTACAACCAAATGCATAAAAtt                                                                                >  1:473833/1‑150 (MQ=255)
                                                                     cacaGTTGCGCCAATTCAGCAGAACGAATGACGATGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTg                                                                               <  1:1143189/150‑1 (MQ=255)
                                                                     cacaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTg                                                                               >  2:1352408/1‑150 (MQ=255)
                                                                      acaTTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCGCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGATATAACAACCAAATGCATAAAATTgg                                                                              <  2:1299390/150‑1 (MQ=255)
                                                                      acaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                                                                              <  1:1203784/150‑1 (MQ=255)
                                                                       aaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                                                                             <  2:1415472/149‑1 (MQ=255)
                                                                         gTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaa                                                                           >  2:55417/1‑150 (MQ=255)
                                                                           tGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaaa                                                                         >  1:1016257/1‑150 (MQ=255)
                                                                               ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGCAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                                                                     >  1:82082/1‑150 (MQ=255)
                                                                               ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                                                                     <  1:500282/150‑1 (MQ=255)
                                                                                cAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAATGGCAg                                                                   <  2:731208/150‑1 (MQ=255)
                                                                                  aTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGa                                                                  <  1:430316/150‑1 (MQ=255)
                                                                                     cagcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGc                                                               >  2:1583487/1‑150 (MQ=255)
                                                                                     cagcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCACATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGc                                                               >  2:123199/1‑150 (MQ=255)
                                                                                            aCGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGACCCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc                                                        >  1:133820/1‑150 (MQ=255)
                                                                                            aCGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc                                                        >  2:430867/1‑150 (MQ=255)
                                                                                            aCGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc                                                        <  1:480758/150‑1 (MQ=255)
                                                                                             cGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGTCAGAAGCGTTTCGCa                                                       >  2:894376/1‑150 (MQ=255)
                                                                                              gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCAt                                                      >  2:1238286/1‑150 (MQ=255)
                                                                                                aTGACGAAGCTCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAATCCTGACCCTCATAAACTGTCTCTATATAACAACCAACTGCATAAAATTGGGAATCTGGAAGAACCGTTTCGCATGc                                                    >  2:316562/1‑150 (MQ=255)
                                                                                                     gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAGCCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGTAAATGGCAGAAGCGTTTCGCATGCGtttt                                               >  2:677851/1‑150 (MQ=255)
                                                                                                      aaGCGCATTTGGTATATTGTCAGGTGGCGTTTGATGAGCAAGCAAGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAAAAACCAAATGCATAAAATTTGTAAAATGGCAGAATCGTTTCGCATGCGttttt                                              >  2:1630622/1‑150 (MQ=255)
                                                                                                      aaGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGatttt                                              >  2:165403/1‑150 (MQ=255)
                                                                                                          gcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCGCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGGTTCGCATGCGGGTTTGAAt                                          <  2:114348/150‑1 (MQ=255)
                                                                                                               tGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTtata                                     >  2:738866/1‑150 (MQ=255)
                                                                                                                ggTATATTGTCAGGGGGCGGTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTatat                                    <  2:923106/150‑1 (MQ=255)
                                                                                                                  tataTTGTCAGGGGGCGTTTGATGAGCATGCTTGGTTGCTCATTGCCATTCTGAAACGTGAACCTCATAAACTGGATCGAGATAACAACGAAATGCATAAAATTGGGAAATTGGCAGAAGCGTTTCGCATGCGTTTTTGATTTTAtattt                                  >  1:162774/1‑149 (MQ=255)
                                                                                                                  tataTTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATACACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTAtatta                                  <  2:1502165/150‑1 (MQ=255)
                                                                                                                  tataTTGTCAGGGGGCGTTTGATGAGCAACCATGGTTGCTCATTGCAATTCTGAAACCTGAACCTCATTAAGTGGCTCGAGATAACATCCAAATGCTTAATATTGGGAAAATGGCAGAAGCGTTTCGCATGCTTTTTTGTATTTttatta                                  >  1:168294/1‑150 (MQ=255)
                                                                                                                    tatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTAAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATg                                <  2:979360/150‑1 (MQ=255)
                                                                                                                       tGTCAGGGGGCGTTTGATGAGCAAGCATGGTTACTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATCGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAAt                             >  1:224863/1‑150 (MQ=255)
                                                                                                                        gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATa                            >  2:1107960/1‑150 (MQ=255)
                                                                                                                           aGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACa                         >  2:1081508/1‑150 (MQ=255)
                                                                                                                               ggCGTTTGATGAGCAAGCATGGGTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACa                     <  2:1306270/150‑1 (MQ=255)
                                                                                                                                gCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACaa                    >  1:712419/1‑150 (MQ=255)
                                                                                                                                  gTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAGCTGAACCTCATAAACTGGCTCGAGATAACAACCTAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAAt                  <  2:983892/150‑1 (MQ=255)
                                                                                                                                    ttGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATaa                <  1:865739/150‑1 (MQ=255)
                                                                                                                                    gcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATaa                >  2:1673384/4‑150 (MQ=255)
                                                                                                                                     tgatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg               >  2:1533981/1‑150 (MQ=255)
                                                                                                                                     tgatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg               >  2:10322/1‑150 (MQ=255)
                                                                                                                                     caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg               >  1:655671/3‑150 (MQ=255)
                                                                                                                                     caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAAAAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg               >  2:242088/3‑150 (MQ=255)
                                                                                                                                      gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGa              <  2:387409/150‑1 (MQ=255)
                                                                                                                                      aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGa              >  1:1196485/2‑150 (MQ=255)
                                                                                                                                       atgaGCAAGCATGGTAGCTCATTGCCATTCTGAAACCTGAACCTCAAAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATTGAAGAAGCGTTTCGCCTGCGTTATTGAATTTATATTATGAATAACATACAAATAAGaa             >  1:283572/1‑150 (MQ=255)
                                                                                                                                        tgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGACCCTCATAAACTGGCTCGAGATAACACCCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTGACGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGaaa            >  2:616655/1‑150 (MQ=255)
                                                                                                                                         gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAAc           <  2:193034/150‑1 (MQ=255)
                                                                                                                                         gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTGTTGAATTTATATTATGAATTACATACAAATAAGAAAc           >  2:1002542/1‑150 (MQ=255)
                                                                                                                                         gaGCAAGAATGGTTGCTCATTGCCATTCTGAAACATGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATTCGTTGTTGAATTTATATTATGAATAAAATACAAATAATaaaa           >  2:696457/1‑149 (MQ=255)
                                                                                                                                           gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTa         >  1:807432/1‑150 (MQ=255)
                                                                                                                                            caagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAAGTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTAt        >  1:1079788/1‑150 (MQ=255)
                                                                                                                                              agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCa      <  2:545097/150‑1 (MQ=255)
                                                                                                                                                  tGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCt  >  2:133849/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCT  >  NC_000913/252446‑252741

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: