Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,958 A→T 6.4% intergenic (+11/+14) 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,9580AT 6.4% 214.9 / 4.5 62intergenic (+11/+14)tqsA/pntBpheromone AI‑2 transporter/pyridine nucleotide transhydrogenase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base T (2/2);  total (32/30)
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.86e-01

TTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCACGTT  >  NC_000913/1674812‑1675106
                                                                                                                                                  |                                                                                                                                                    
tttttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCaga                                                                                                                                                    <  1:657710/149‑1 (MQ=255)
 ttttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCagag                                                                                                                                                   <  3:467441/149‑1 (MQ=255)
 ttttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCagag                                                                                                                                                   <  2:795634/149‑1 (MQ=255)
   ttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGc                                                                                                                                                 <  4:579600/149‑1 (MQ=255)
    tGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGcc                                                                                                                                                <  4:254828/149‑1 (MQ=255)
             ttgttAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGtt                                                                                                                                       <  1:716382/149‑1 (MQ=255)
                  aGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACaga                                                                                                                                  >  1:222504/1‑149 (MQ=255)
                     aCCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCt                                                                                                                               <  2:222504/149‑1 (MQ=255)
                         gTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCa                                                                                                                           <  3:86077/149‑1 (MQ=255)
                              tATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGAtt                                                                                                                      <  4:511305/149‑1 (MQ=255)
                                tgctgcTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGc                                                                                                                    >  2:730087/1‑149 (MQ=255)
                                 gctgctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAgg                                                                                                                                          <  2:40372/126‑1 (MQ=255)
                                 gctgctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAgg                                                                                                                                          >  1:105653/1‑126 (MQ=255)
                                 gctgctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAgg                                                                                                                                          >  1:40372/1‑126 (MQ=255)
                                 gctgctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAgg                                                                                                                                          <  2:105653/126‑1 (MQ=255)
                                   tgctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATc                                                                                                                 <  1:128344/149‑1 (MQ=255)
                                    gctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGCTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATcc                                                                                                                >  1:691888/1‑149 (MQ=255)
                                        tCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGc                                                                                                            >  4:694563/1‑149 (MQ=255)
                                           gTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGGATCCACGCTgg                                                                                                         <  1:426741/149‑1 (MQ=255)
                                              ccGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCtt                                                                                                      <  4:278039/149‑1 (MQ=255)
                                                    aCAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGcc                                                                                                                                                <  2:206088/101‑1 (MQ=255)
                                                    aCAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGcc                                                                                                                                                >  1:206088/1‑101 (MQ=255)
                                                              tCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAg                                                                                      >  3:248695/1‑149 (MQ=255)
                                                              tCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAg                                                                                      >  2:484468/1‑149 (MQ=255)
                                                                aaaaTTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCa                                                                                    >  4:299053/1‑149 (MQ=255)
                                                                     tGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTgg                                                                               <  3:107341/149‑1 (MQ=255)
                                                                                aaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCt                                                                     >  1:500944/1‑148 (MQ=255)
                                                                                aaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCt                                                                     >  1:790702/1‑148 (MQ=255)
                                                                                aaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCt                                                                     <  2:790702/148‑1 (MQ=255)
                                                                                aaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCt                                                                     <  2:500944/148‑1 (MQ=255)
                                                                                aaacaGCGGGAGGTCAAAGCATCGCCGTTCTGCTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCtt                                                                    <  2:691888/149‑1 (MQ=255)
                                                                                               aaaGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACaga                                                                                                                                  >  1:29030/1‑72 (MQ=255)
                                                                                               aaaGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACaga                                                                                                                                  <  2:29030/72‑1 (MQ=255)
                                                                                                aaGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGttt                                                    <  4:370683/149‑1 (MQ=255)
                                                                                                 aGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCAtgt                                                                                 >  4:517562/1‑119 (MQ=255)
                                                                                                 aGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCAtgt                                                                                 <  3:517562/119‑1 (MQ=255)
                                                                                                  gCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTg                                                  <  1:730087/149‑1 (MQ=255)
                                                                                                  gCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTg                                                  <  3:299053/149‑1 (MQ=255)
                                                                                                   cATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGc                                                 >  1:413178/1‑149 (MQ=255)
                                                                                                           ttCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGca                                                <  4:254362/142‑1 (MQ=255)
                                                                                                           ttCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGca                                                >  3:254362/1‑142 (MQ=255)
                                                                                                               gTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAg                                     >  3:201967/1‑149 (MQ=255)
                                                                                                                       aTCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGtt                             >  1:588601/1‑149 (MQ=255)
                                                                                                                          tcAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCAt                          >  3:457909/1‑149 (MQ=255)
                                                                                                                            aaTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGc                                    >  1:173414/1‑137 (MQ=255)
                                                                                                                            aaTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGCGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGc                                    <  2:173414/137‑1 (MQ=255)
                                                                                                                              tAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGaa                      <  4:457909/149‑1 (MQ=255)
                                                                                                                                aaGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGtttt                                                   >  1:221970/1‑118 (MQ=255)
                                                                                                                                aaGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGtttt                                                   <  2:221970/118‑1 (MQ=255)
                                                                                                                                  gagTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGtt                  >  1:639064/1‑149 (MQ=255)
                                                                                                                                   cctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAgtg                               <  1:13503/133‑1 (MQ=255)
                                                                                                                                   cctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAgtg                               >  2:13503/3‑135 (MQ=255)
                                                                                                                                   agTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGttt                 >  3:53739/1‑149 (MQ=255)
                                                                                                                                      gACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTaaa              <  4:248695/149‑1 (MQ=255)
                                                                                                                                       aCGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAg             >  2:147662/1‑149 (MQ=255)
                                                                                                                                        cGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGa            <  1:147662/149‑1 (MQ=255)
                                                                                                                                         ggCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGtttt                                                   >  1:161469/1‑109 (MQ=255)
                                                                                                                                         ggCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGtttt                                                   <  2:161469/109‑1 (MQ=255)
                                                                                                                                          gCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACa          >  2:381818/1‑149 (MQ=255)
                                                                                                                                           ccTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACaa         >  1:36917/1‑149 (MQ=255)
                                                                                                                                            cTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaa                                                                 >  4:589188/1‑92 (MQ=255)
                                                                                                                                            cTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaa                                                                 <  3:589188/92‑1 (MQ=255)
                                                                                                                                             tCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATc       >  4:210839/1‑149 (MQ=255)
                                                                                                                                             tCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATc       >  3:682472/1‑149 (MQ=255)
                                                                                                                                             tCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATc       >  4:753208/1‑149 (MQ=255)
                                                                                                                                                  agagGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCACGtt  >  4:186947/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCACGTT  >  NC_000913/1674812‑1675106

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