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

GCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCT  >  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)
           tGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGa                                                                                                                                       <  4:597274/149‑1 (MQ=255)
            gCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCG‑TCAACTCGTTGAGCGCGGaa                                                                                                                                      <  2:628459/149‑2 (MQ=255)
                ttATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGGGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTa                                                                                                                                  >  1:547489/1‑149 (MQ=255)
                ttATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTa                                                                                                                                  >  1:518596/1‑149 (MQ=255)
                 tATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTAc                                                                                                                                 >  2:95088/1‑149 (MQ=255)
                     cTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCa                                                                                                                             >  1:547276/1‑149 (MQ=255)
                      ttAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCaa                                                                                                                            <  1:98971/149‑1 (MQ=255)
                         aCCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAatat                                                                                                                         >  1:515394/1‑149 (MQ=255)
                              cTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGACGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCaa                                                                                                                    <  4:275557/149‑1 (MQ=255)
                               tGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAAc                                                                                                                   <  2:629876/149‑1 (MQ=255)
                                   gCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCta                                                                                                               >  2:335646/1‑149 (MQ=255)
                                     aaCACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCtatt                                                                                                             <  2:698616/149‑1 (MQ=255)
                                     aaCACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCtatt                                                                                                             >  1:374896/1‑149 (MQ=255)
                                     aaCACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCtatt                                                                                                             <  4:206792/149‑1 (MQ=255)
                                            ggTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgc                                                                                                      <  4:76397/149‑1 (MQ=255)
                                             gTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgcc                                                                                                     >  1:451581/1‑149 (MQ=255)
                                             gTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgcc                                                                                                     >  3:86597/1‑149 (MQ=255)
                                              tCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgccg                                                                                                    >  3:87408/1‑149 (MQ=255)
                                              tCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGAcgccg                                                                                                    >  1:744589/1‑149 (MQ=255)
                                                         gtgCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCt                                                                                         <  3:298398/149‑1 (MQ=255)
                                                           gCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTg                                                                                       <  3:619374/149‑1 (MQ=255)
                                                              gcgTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAAc                                                                                    <  1:397372/149‑1 (MQ=255)
                                                                  gtACTGGTGGTGTCGATTTTGATGCCCATCATCGTGCTGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgtt                                                                                <  1:309602/149‑1 (MQ=255)
                                                                   tACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttg                                                                               >  1:448053/1‑149 (MQ=255)
                                                                   tACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttg                                                                               <  2:249723/149‑1 (MQ=255)
                                                                              tCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGt                                                                    >  4:591872/1‑149 (MQ=255)
                                                                               cGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTa                                                                   <  4:96629/149‑1 (MQ=255)
                                                                                gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                  <  2:547276/149‑1 (MQ=255)
                                                                                gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                  <  3:237513/149‑1 (MQ=255)
                                                                                gATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAg                                                                  >  2:517748/1‑149 (MQ=255)
                                                                                 aTTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAgg                                                                 <  2:465362/149‑1 (MQ=255)
                                                                                 aTTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAgg                                                                 <  2:575007/149‑1 (MQ=255)
                                                                                  ttttGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAggg                                                                >  1:159483/1‑149 (MQ=255)
                                                                                     tgatgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGAtt                                                             <  1:492561/149‑1 (MQ=255)
                                                                                       atgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGa                                                           <  3:415432/149‑1 (MQ=255)
                                                                                        tgaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTAGCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCTGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGAc                                                          >  2:378888/1‑149 (MQ=255)
                                                                                         gaCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTAAGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACg                                                         >  2:467862/1‑149 (MQ=255)
                                                                                          aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGAc                                                                                                        <  3:605422/101‑1 (MQ=255)
                                                                                          aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGAc                                                                                                        >  4:605422/1‑101 (MQ=255)
                                                                                          aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGt                                                        <  2:451581/149‑1 (MQ=255)
                                                                                          aCCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGt                                                        <  2:204861/149‑1 (MQ=255)
                                                                                             atcatcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAActc                                                                                                                 >  2:287023/1‑89 (MQ=255)
                                                                                             atcatcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAActc                                                                                                                 <  1:287023/89‑1 (MQ=255)
                                                                                                 tcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGtt                                                 <  3:658328/149‑1 (MQ=255)
                                                                                                 tcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGtt                                                 <  3:52865/149‑1 (MQ=255)
                                                                                                 tcGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGATGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGtt                                                 >  1:214929/1‑149 (MQ=255)
                                                                                                                  ttGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                >  1:723006/1‑149 (MQ=255)
                                                                                                                  ttGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATta                                >  3:362134/1‑149 (MQ=255)
                                                                                                                    gCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATtata                              <  1:467862/149‑1 (MQ=255)
                                                                                                                      tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                             >  1:282150/1‑100 (MQ=255)
                                                                                                                      tATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                             <  2:282150/100‑1 (MQ=255)
                                                                                                                       aTTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTg                           >  2:373963/1‑149 (MQ=255)
                                                                                                                        ttAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGa                          <  1:652394/149‑1 (MQ=255)
                                                                                                                          aGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATc                        <  3:591872/149‑1 (MQ=255)
                                                                                                                             ttATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGa                     >  3:669615/1‑149 (MQ=255)
                                                                                                                                tCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcg                  >  3:658791/1‑149 (MQ=255)
                                                                                                                                  tGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcgcg                >  3:416635/1‑149 (MQ=255)
                                                                                                                                  tGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAAcgcg                >  4:99389/1‑149 (MQ=255)
                                                                                                                                      ttCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGa            >  2:380318/1‑149 (MQ=255)
                                                                                                                                      gtCCGCG‑TCAACTCGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt           >  4:700619/2‑149 (MQ=255)
                                                                                                                                       tCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGAt           <  4:508165/149‑1 (MQ=255)
                                                                                                                                       tCCGCG‑TCAACTCGTTGA‑‑GCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                  <  2:582516/108‑1 (MQ=255)
                                                                                                                                       tCCGCG‑TCAACTCGTTGA‑‑GCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGt                                                  >  1:582516/1‑108 (MQ=255)
                                                                                                                                          gtgcTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGGAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGAc        >  3:591967/3‑149 (MQ=255)
                                                                                                                                            gcTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCa                                                    >  3:588360/1‑103 (MQ=255)
                                                                                                                                            gcTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCa                                                    <  4:588360/103‑1 (MQ=255)
                                                                                                                                               cAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGc   >  2:50427/1‑149 (MQ=255)
                                                                                                                                                aaCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                             <  3:222356/74‑1 (MQ=255)
                                                                                                                                                aaCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCgttgtt                                                                             >  4:222356/1‑74 (MQ=255)
                                                                                                                                                aaCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt  >  1:51739/1‑149 (MQ=255)
                                                                                                                                                aaCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt  >  2:20042/1‑149 (MQ=255)
                                                                                                                                                aaCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt  >  3:196048/1‑149 (MQ=255)
                                                                                                                                                aaCGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt  >  3:276522/1‑149 (MQ=255)
                                                                                                                                                aaCGAGTTGA‑CGCGGACGGTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCt  >  2:287745/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCATTATTTATTGCTGTTATTCTTAACCCGCTGGTGCAACACATGGTCCGCTGGCGTGTGCCGCGTGTACTGGCGGTGTCGATTTTGATGACCATCATCGTGATGGCGATGGTGTTGCTATTAGCTTATCTGGGTTCCGCGCTCAACGAGTTGA‑CGCGGACGTTACCGCAATATCGCAACTCTATTATGACGCCGCTGCAAGCGCTTGAACCGTTGTTGCAACGCGTAGGGATTGACGTCTCAGTTGACCAGCTGGCGCATTATATTGATCCGAACGCGGCGATGACGTTGCT  >  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: ‑