Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,955 A→T 5.7% intergenic (+8/+17) tqsA → / ← pntB pheromone AI‑2 transporter/pyridine nucleotide transhydrogenase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,674,9550AT 5.7% 218.3 / 2.8 70intergenic (+8/+17)tqsA/pntBpheromone AI‑2 transporter/pyridine nucleotide transhydrogenase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (36/29);  new base T (2/2);  total (39/31)
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.92e-01

TTGATTTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAAT  >  NC_000913/1674807‑1675100
                                                                                                                                                    |                                                                                                                                                 
ttGATTTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTca                                                                                                                                                   <  7:25247/149‑1 (MQ=255)
   aTTTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTATGTTAAGCGATCTCAATAAAGAGTGACGGCCTcag                                                                                                                                                  >  3:250454/1‑147 (MQ=255)
       tttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCagag                                                                                                                                             <  5:121396/148‑1 (MQ=255)
           gggATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGt                                                                                                                                        >  2:291486/1‑149 (MQ=255)
                  ttgttAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGtt                                                                                                                                 <  4:250454/149‑1 (MQ=255)
                    gttAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTa                                                                                                                                >  1:300498/1‑148 (MQ=255)
                         gACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGc                                                                                                                          >  3:322277/1‑149 (MQ=255)
                            cGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTcagc                                                                                                                                                 >  6:94821/1‑123 (MQ=255)
                            cGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTcagc                                                                                                                                                 <  5:94821/123‑1 (MQ=255)
                                  gTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGAt                                                                                                                 >  7:229061/1‑149 (MQ=255)
                                   tATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACag                                                                                                                             >  1:216628/1‑136 (MQ=255)
                                   tATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACag                                                                                                                             <  2:216628/136‑1 (MQ=255)
                                          ctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCa                                                                                                         >  6:6145/1‑149 (MQ=255)
                                               cGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCt                                                                                                     <  1:144943/148‑1 (MQ=255)
                                                    cGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCttt                                                                                               >  2:359738/1‑149 (MQ=255)
                                                      ttGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTgg                                                                                             <  6:110386/149‑1 (MQ=255)
                                                         aCAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAg                                                                                                                    <  8:166236/123‑1 (MQ=255)
                                                         aCAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAg                                                                                                                    >  7:166236/1‑123 (MQ=255)
                                                         aCAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGt                                                                                          >  4:136935/1‑149 (MQ=255)
                                                         aCAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGt                                                                                          >  3:236544/1‑149 (MQ=255)
                                                           aaTTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCGAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCa                                                                                        >  6:148346/1‑149 (MQ=255)
                                                           aaTTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCa                                                                                        >  6:205763/1‑149 (MQ=255)
                                                            attattGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCCGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCAc                                                                                       >  5:312308/1‑149 (MQ=255)
                                                                ttGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCaaa                                                                                   >  1:222081/1‑149 (MQ=255)
                                                                   tCAAAATTGCGCTTGAACAAACAGCGGGATGTCAAAGCATCGCCGTACTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGCATTGCATCCACGCTTGCTTTGGAGTCACCAAACAg                                                                                >  6:102766/1‑149 (MQ=255)
                                                                   tCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACa                                                                                 <  5:205763/148‑1 (MQ=255)
                                                                   tCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTCCAACCACGCTGGCTTTGGCGTCACCAAACAg                                                                                >  7:345196/1‑149 (MQ=255)
                                                                        aTTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCAtgt                                                                           <  1:291486/149‑1 (MQ=255)
                                                                               ttGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGtt                                                                    >  8:332252/1‑149 (MQ=255)
                                                                                 gAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGtttt                                                                  >  5:346876/1‑149 (MQ=255)
                                                                                  aacaaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTc                                                                 >  6:103179/1‑149 (MQ=255)
                                                                                   acaaacaGCGGGAGGTCAAAGCATCGCCGGTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGt                                                                                                                                        <  3:83841/77‑1 (MQ=255)
                                                                                     aaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCtt                                                              >  2:221636/1‑149 (MQ=255)
                                                                                        caGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaa                                                           <  4:236544/149‑1 (MQ=255)
                                                                                        caGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaa                                                           <  1:359738/149‑1 (MQ=255)
                                                                                           cGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaagaa                                                        <  4:322277/149‑1 (MQ=255)
                                                                                             ggAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTAGGTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACa                                                      <  6:239379/149‑1 (MQ=255)
                                                                                                   cAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGc                                                                                                  <  2:366292/99‑1 (MQ=255)
                                                                                                    aaaGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGtt                                               >  4:155128/1‑149 (MQ=255)
                                                                                                     aaGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGttt                                              <  2:222081/149‑1 (MQ=255)
                                                                                                     aaGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGttt                                              <  2:300498/149‑1 (MQ=255)
                                                                                                      aGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCAc                                                                                       >  6:140343/1‑107 (MQ=255)
                                                                                                      aGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGtttt                                             <  5:103179/149‑1 (MQ=255)
                                                                                                      aGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGtttt                                             <  6:30182/149‑1 (MQ=255)
                                                                                                       gCATCGCCGTTCTGTTAAGCGATCTCAATATAGAGTGACGGCCTGAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTg                                            <  6:51596/149‑1 (MQ=255)
                                                                                                         aTCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGca                                          <  8:345196/149‑1 (MQ=255)
                                                                                                          tCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGcac                                         <  5:148346/149‑1 (MQ=255)
                                                                                                          tCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGca                                          >  1:22355/1‑148 (MQ=255)
                                                                                                          tCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGca                                          >  3:346824/1‑148 (MQ=255)
                                                                                                                ttCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCtt                                                              >  6:263986/1‑122 (MQ=255)
                                                                                                                ttCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCtt                                                              <  5:263986/122‑1 (MQ=255)
                                                                                                                     ttAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCTCCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATa                                <  5:102766/147‑1 (MQ=255)
                                                                                                                        aGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCACACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATCGCCAg                           >  5:51062/1‑149 (MQ=255)
                                                                                                                        aGCGATCGCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCa                            <  4:243256/148‑1 (MQ=255)
                                                                                                                            aTCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGtt                       >  5:195225/1‑149 (MQ=255)
                                                                                                                            aTCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCAAGTGtt                       >  4:230120/1‑149 (MQ=255)
                                                                                                                            aTCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGGCACCAAACAGCATGTGGGTGTTTTCCtt                                                              >  1:35172/1‑110 (MQ=255)
                                                                                                                            aTCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGCCACCAAACAGCATGTGGGTGTTTTCCtt                                                              <  2:35172/110‑1 (MQ=255)
                                                                                                                             tctcAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGtt                       >  6:302084/1‑148 (MQ=255)
                                                                                                                              ctcAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCa                     >  2:406600/1‑149 (MQ=255)
                                                                                                                                 aaTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCa                     <  1:20855/146‑1 (MQ=255)
                                                                                                                                 aaTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCa                     >  2:20855/1‑146 (MQ=255)
                                                                                                                                 aaTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATc                   >  1:235747/1‑148 (MQ=255)
                                                                                                                                 aaTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATc                   <  2:22355/148‑1 (MQ=255)
                                                                                                                                       ccctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGtt            >  1:108180/4‑149 (MQ=255)
                                                                                                                                         ctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTCTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTa          <  6:388107/148‑1 (MQ=255)
                                                                                                                                          tGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTa          >  6:175770/1‑148 (MQ=255)
                                                                                                                                          tGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTaa         <  3:155128/149‑1 (MQ=255)
                                                                                                                                          tGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTaa         <  7:332252/149‑1 (MQ=255)
                                                                                                                                            aCGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCTTGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTaaa        <  3:204771/148‑1 (MQ=255)
                                                                                                                                              ggCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGa      >  1:170608/1‑148 (MQ=255)
                                                                                                                                              ggCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGa      >  8:60706/1‑148 (MQ=255)
                                                                                                                                                 cTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCCCCAAACAGCATGTGGGTGTTTTCCTTGAAGAACACCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAAt  >  1:328923/1‑149 (MQ=255)
                                                                                                                                                 cTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAAt  >  4:189236/1‑149 (MQ=255)
                                                                                                                                                    agcagAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAg                               <  3:100543/117‑1 (MQ=255)
                                                                                                                                                    agcagAGGCCGTCAGGGTTACAGAGCTTTCAGAATTGCATCCACGCTGGCTTTGTCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAg                               >  4:100543/1‑117 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTGATTTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAAT  >  NC_000913/1674807‑1675100

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