Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,178 A→C 5.3% E89A (GAG→GCG) ‡ tqsA → pheromone AI‑2 transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,674,1780AC 5.3% 262.0 / 6.4 76E89A (GAG→GCG) ‡tqsApheromone AI‑2 transporter
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base C (2/2);  total (40/36)
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.10e-01

GCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCT  >  NC_000913/1674030‑1674322
                                                                                                                                                    |                                                                                                                                                 
gCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGa                                                                                                                                                   <  2:240412/149‑1 (MQ=255)
    tatttattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTg                                                                                                                                               <  2:134354/149‑1 (MQ=255)
      tttattGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTggc                                                                                                                                             <  2:360170/149‑3 (MQ=255)
           tGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGa                                                                                                                                       <  4:597274/149‑1 (MQ=255)
            gCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGTCAACTCGTTGAGCGCGGaa                                                                                                                                      <  2:628459/149‑2 (MQ=255)
                ttATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGGGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTa                                                                                                                                  >  1:547489/1‑149 (MQ=255)
                ttATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTa                                                                                                                                  >  1:518596/1‑149 (MQ=255)
                 tATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTAc                                                                                                                                 >  2:95088/1‑149 (MQ=255)
                     cTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCa                                                                                                                             >  1:547276/1‑149 (MQ=255)
                      ttAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCaa                                                                                                                            <  1:98971/149‑1 (MQ=255)
                         aCCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAatat                                                                                                                         >  1:515394/1‑149 (MQ=255)
                              cTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGACGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCaa                                                                                                                    <  4:275557/149‑1 (MQ=255)
                               tGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAAc                                                                                                                   <  2:629876/149‑1 (MQ=255)
                                   gCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCta                                                                                                               >  2:335646/1‑149 (MQ=255)
                                     aaCACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCtatt                                                                                                             <  2:698616/149‑1 (MQ=255)
                                     aaCACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCtatt                                                                                                             >  1:374896/1‑149 (MQ=255)
                                     aaCACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCtatt                                                                                                             <  4:206792/149‑1 (MQ=255)
                                            ggTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgc                                                                                                      <  4:76397/149‑1 (MQ=255)
                                             gTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgcc                                                                                                     >  1:451581/1‑149 (MQ=255)
                                             gTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgcc                                                                                                     >  3:86597/1‑149 (MQ=255)
                                              tCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgccg                                                                                                    >  3:87408/1‑149 (MQ=255)
                                              tCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgccg                                                                                                    >  1:744589/1‑149 (MQ=255)
                                                         gtgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCt                                                                                         <  3:298398/149‑1 (MQ=255)
                                                           gCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTg                                                                                       <  3:619374/149‑1 (MQ=255)
                                                              gcgTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAAc                                                                                    <  1:397372/149‑1 (MQ=255)
                                                                  gtACTGGTGGTGTCGATTTTGATGCCCATCATCGTGCTGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgtt                                                                                <  1:309602/149‑1 (MQ=255)
                                                                   tACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttg                                                                               >  1:448053/1‑149 (MQ=255)
                                                                   tACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttg                                                                               <  2:249723/149‑1 (MQ=255)
                                                                              tCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGt                                                                    >  4:591872/1‑149 (MQ=255)
                                                                               cGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTa                                                                   <  4:96629/149‑1 (MQ=255)
                                                                                gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                  <  2:547276/149‑1 (MQ=255)
                                                                                gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                  <  3:237513/149‑1 (MQ=255)
                                                                                gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                  >  2:517748/1‑149 (MQ=255)
                                                                                 aTTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAgg                                                                 <  2:465362/149‑1 (MQ=255)
                                                                                 aTTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAgg                                                                 <  2:575007/149‑1 (MQ=255)
                                                                                  ttttGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAggg                                                                >  1:159483/1‑149 (MQ=255)
                                                                                     tgatgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGAtt                                                             <  1:492561/149‑1 (MQ=255)
                                                                                       atgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGa                                                           <  3:415432/149‑1 (MQ=255)
                                                                                        tgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTAGCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCTGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGAc                                                          >  2:378888/1‑149 (MQ=255)
                                                                                         gaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTAAGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACg                                                         >  2:467862/1‑149 (MQ=255)
                                                                                          aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAc                                                                                                        <  3:605422/101‑1 (MQ=255)
                                                                                          aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGAc                                                                                                        >  4:605422/1‑101 (MQ=255)
                                                                                          aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGt                                                        <  2:451581/149‑1 (MQ=255)
                                                                                          aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGt                                                        <  2:204861/149‑1 (MQ=255)
                                                                                             atcatcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAActc                                                                                                                 >  2:287023/1‑89 (MQ=255)
                                                                                             atcatcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAActc                                                                                                                 <  1:287023/89‑1 (MQ=255)
                                                                                                 tcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGtt                                                 <  3:658328/149‑1 (MQ=255)
                                                                                                 tcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGtt                                                 <  3:52865/149‑1 (MQ=255)
                                                                                                 tcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGATGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGtt                                                 >  1:214929/1‑149 (MQ=255)
                                                                                                                  ttGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                >  1:723006/1‑149 (MQ=255)
                                                                                                                  ttGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                >  3:362134/1‑149 (MQ=255)
                                                                                                                    gCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATtata                              <  1:467862/149‑1 (MQ=255)
                                                                                                                      tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                             >  1:282150/1‑100 (MQ=255)
                                                                                                                      tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                             <  2:282150/100‑1 (MQ=255)
                                                                                                                       aTTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTg                           >  2:373963/1‑149 (MQ=255)
                                                                                                                        ttAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGa                          <  1:652394/149‑1 (MQ=255)
                                                                                                                          aGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATc                        <  3:591872/149‑1 (MQ=255)
                                                                                                                             ttATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGa                     >  3:669615/1‑149 (MQ=255)
                                                                                                                                tCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcg                  >  3:658791/1‑149 (MQ=255)
                                                                                                                                  tGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcgcg                >  3:416635/1‑149 (MQ=255)
                                                                                                                                  tGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcgcg                >  4:99389/1‑149 (MQ=255)
                                                                                                                                      ttCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGa            >  2:380318/1‑149 (MQ=255)
                                                                                                                                      gtCCGCGTCAACTCGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt           >  4:700619/2‑149 (MQ=255)
                                                                                                                                       tCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt           <  4:508165/149‑1 (MQ=255)
                                                                                                                                       tCCGCGTCAACTCGTTGAGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                  <  2:582516/108‑1 (MQ=255)
                                                                                                                                       tCCGCGTCAACTCGTTGAGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                  >  1:582516/1‑108 (MQ=255)
                                                                                                                                          gtgcTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGGAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGAc        >  3:591967/3‑149 (MQ=255)
                                                                                                                                            gcTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCa                                                    >  3:588360/1‑103 (MQ=255)
                                                                                                                                            gcTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCa                                                    <  4:588360/103‑1 (MQ=255)
                                                                                                                                               cAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGc   >  2:50427/1‑149 (MQ=255)
                                                                                                                                                aaCGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                             <  3:222356/74‑1 (MQ=255)
                                                                                                                                                aaCGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                             >  4:222356/1‑74 (MQ=255)
                                                                                                                                                aaCGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt  >  1:51739/1‑149 (MQ=255)
                                                                                                                                                aaCGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt  >  2:20042/1‑149 (MQ=255)
                                                                                                                                                aaCGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt  >  3:196048/1‑149 (MQ=255)
                                                                                                                                                aaCGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt  >  3:276522/1‑149 (MQ=255)
                                                                                                                                                aaCGAGTTGACGCGGACGGTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt  >  2:287745/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGACGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCT  >  NC_000913/1674030‑1674322

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: