Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,171,863 A→G 17.2% intergenic (+30/+16) qseC → / ← ygiZ quorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,171,8630AG17.2% 333.1 / 63.6 157intergenic (+30/+16)qseC/ygiZquorum sensing sensory histidine kinase in two‑component regulatory system with QseB/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (72/58);  new base G (15/12);  total (87/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

Only 100 of 157 total aligned reads displayed.

TTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAA  >  NC_000913/3171765‑3171960
                                                                                                  |                                                                                                 
  gcatggcatgAATGTTGAATTTGGGAATGCGGAACAAAGTGGGTTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGActc                                                                                                <  1:2677306/100‑1 (MQ=255)
    atggcatgAATGTTGAATTTGGGAATGCGGAACAAGCTGGATTTGAGGCGAAGGTAAGCTGGTATCTTCTCTGTCTATGAGAGCCGTTAAAACGACtctc                                                                                              <  1:1385795/100‑1 (MQ=255)
     tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCa                                                                                             <  2:813171/100‑1 (MQ=38)
     tggcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCa                                                                                             <  1:2018128/100‑1 (MQ=255)
       gcatgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATa                                                                                           <  2:479444/100‑1 (MQ=38)
        catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAg                                                                                          <  2:772190/100‑1 (MQ=38)
        catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAg                                                                                          >  2:3101059/1‑100 (MQ=255)
         atgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGa                                                                                         <  2:606174/100‑1 (MQ=38)
          tgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGTCGTTTTAACGGCTCTCATAGac                                                                                        <  2:2294402/100‑2 (MQ=37)
                ttGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACt                                                                                  <  2:1220015/100‑1 (MQ=255)
                 tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                                 >  1:866753/1‑100 (MQ=255)
                 tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                                 >  2:75899/1‑100 (MQ=255)
                 tGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                                 >  2:432104/1‑100 (MQ=255)
                  gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                                >  1:1259122/1‑100 (MQ=255)
                  gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                                >  1:88861/1‑100 (MQ=255)
                  gAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGCAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTaa                                                                                >  1:806306/1‑100 (MQ=255)
                  aaaTTTGAGAATGCGGAACAAGATGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTa                                                                                 >  1:2476698/2‑99 (MQ=255)
                    aTTTGGGAATTCGGAACAAGGGGGTTTGGAGGCGGAGTTAAGCTGGTAACTTCTCTGTCTATGAGATCCGTTAAAACTACTCTCATAGATTTTACTAATa                                                                              >  1:3183933/1‑100 (MQ=255)
                    aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                              >  2:2718794/1‑100 (MQ=255)
                    aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATa                                                                              <  2:524706/100‑1 (MQ=255)
                    aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACGCTCATAGATTTTACTAATa                                                                              >  1:1538786/1‑100 (MQ=255)
                    aTTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTGAAAACGACTCTCATAGATTTTACTAATa                                                                              >  2:1677168/1‑100 (MQ=255)
                      ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                            <  1:336346/100‑1 (MQ=255)
                      ttGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGc                                                                            <  2:2885540/100‑1 (MQ=255)
                        gggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCa                                                                           >  1:1093152/1‑99 (MQ=255)
                         ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                         >  2:2526673/1‑100 (MQ=255)
                         ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaa                                                                         >  2:1296035/1‑100 (MQ=255)
                           aaTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAat                                                                       >  2:2758829/1‑100 (MQ=255)
                           aaTGCGGAACAACGTGGATTTGAGTCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCaaaa                                                                        <  1:3682378/99‑1 (MQ=255)
                            aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                      >  2:1402778/1‑100 (MQ=255)
                            aTGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAata                                                                      <  1:2246134/100‑1 (MQ=255)
                             tGCGGAACAATGTGGATTTGAGGCGAAGGTAAGCTTGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCTTAGATTTTACTAATAGCAAAatat                                                                     <  1:2425276/100‑1 (MQ=255)
                              gCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAAtata                                                                    <  2:3658949/100‑1 (MQ=255)
                               cGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaa                                                                   >  1:1453044/1‑100 (MQ=255)
                                ggAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATaaa                                                                  >  2:1517357/1‑100 (MQ=255)
                                 gcaccaGGTGGATTTGAGGCGAAGGTAAGCTGGTCACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCGCATACATTTTACTAATATCAAAATATAAAc                                                                 <  1:702004/95‑1 (MQ=255)
                                 gAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAAc                                                                 >  2:1447626/1‑100 (MQ=255)
                                       ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                           <  2:2981177/100‑1 (MQ=255)
                                       ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTccc                                                           >  1:564830/1‑100 (MQ=255)
                                           gATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaa                                                       >  2:23555/1‑100 (MQ=255)
                                            aTTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCaaaa                                                      >  1:1761081/1‑100 (MQ=255)
                                               tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCGGTTAAAACGACACACATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                   >  1:528613/1‑100 (MQ=255)
                                               tGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                   >  2:1417216/1‑100 (MQ=255)
                                               tGAGGCGAAGGTAAACTGGTAACTTCTCTGTCTATGAGAGCAGTAAAAACGACTATCATAGCTCTTACTACTAGCAAAAAAGAAACAGCACCCaaaaaac                                                   >  1:3539741/1‑99 (MQ=255)
                                                  ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                >  1:2259522/1‑100 (MQ=255)
                                                  ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGcca                                                >  2:3717921/1‑100 (MQ=255)
                                                   gCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGccac                                               <  2:3688748/100‑1 (MQ=255)
                                                      aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaa                                            >  2:3354339/1‑100 (MQ=255)
                                                        ggTAAGCTGGTAACTTCTCTGTCTATGCGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                          <  2:1510251/100‑1 (MQ=255)
                                                        ggTAAGCTGGTAACTAGTCTGTCTATGAGAGCCGTTAAAACGACTCACATAGATTTTACTAATAGCAAAATATAAAAAAAAAAAAAAAAAGCCaccaacc                                          <  1:782290/100‑1 (MQ=255)
                                                         gcaaGCTGGAACCTTCTCTGTCTATGAGAGCCGTTACAACGACTCTCATAGATTTTACTAATTGCAAAACATAAAACGTCCCAAAAAAAGCCaccaacc                                          <  1:316732/97‑1 (MQ=255)
                                                           aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACcaca                                       <  1:1115252/100‑1 (MQ=255)
                                                           aaGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACCCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCaccaacc                                          <  1:3427823/97‑1 (MQ=255)
                                                             gCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACaaa                                     >  1:556773/1‑100 (MQ=255)
                                                              cTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                    <  2:513480/100‑1 (MQ=255)
                                                               tGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAcc                                   <  2:1087515/100‑1 (MQ=255)
                                                                ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                  >  1:2723420/1‑100 (MQ=255)
                                                                ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                  >  1:2729946/1‑100 (MQ=255)
                                                                  tAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGa                                >  2:3786860/1‑100 (MQ=255)
                                                                   aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                               <  2:2126094/100‑1 (MQ=255)
                                                                   aaCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGAc                               <  2:3600884/100‑1 (MQ=255)
                                                                    aCTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACa                              <  1:2758815/100‑1 (MQ=255)
                                                                     cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATCGATTTTACTAATAGCAAAATATAAACAGTCCCCAAAAAAGCAACGAACCACAAACCAGACaa                             >  1:1061096/1‑100 (MQ=255)
                                                                     cTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACaa                             <  1:3788438/100‑1 (MQ=255)
                                                                         aaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc                         >  2:1883692/5‑100 (MQ=37)
                                                                         aaaatCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGCCCCCAAAAAAGCCACCAACCACAAACCAGACAAATAc                         >  1:803443/5‑100 (MQ=25)
                                                                            gTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                      <  2:1849575/100‑1 (MQ=255)
                                                                            atCTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcga                      >  2:2297272/2‑100 (MQ=37)
                                                                              cTATGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                    >  2:3184877/1‑100 (MQ=38)
                                                                              cTATGAGAGTCGTATTAACGGATCGCATAGAATTTACTAATGGCAAAATACAAAACATCCCCAAAAAAACAAACAACCACAAACCAGACCAGTACcga                      >  1:706315/1‑98 (MQ=255)
                                                                              cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                    <  2:337389/100‑1 (MQ=255)
                                                                              cTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                    >  2:3217476/1‑100 (MQ=255)
                                                                               tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                   <  2:526026/100‑1 (MQ=255)
                                                                               tATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcg                   >  1:1374741/1‑100 (MQ=255)
                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                  >  1:2751599/1‑100 (MQ=38)
                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                  >  2:883412/1‑100 (MQ=38)
                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                  >  2:2916841/1‑100 (MQ=38)
                                                                                aTGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACAAGACAAATACcgaacga                  >  1:246703/1‑100 (MQ=37)
                                                                                aTGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcga                  >  2:693625/1‑100 (MQ=255)
                                                                                 tGAGAGTCGTTTTAACGGCTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                 >  2:1966230/1‑100 (MQ=38)
                                                                                 tGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                 <  1:3082756/100‑1 (MQ=255)
                                                                                  gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCACCCACAAACCAGACAAATACCGAGCgaga                >  2:343037/1‑100 (MQ=255)
                                                                                  gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                >  1:48370/1‑100 (MQ=255)
                                                                                  gagaGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                <  2:1749451/100‑1 (MQ=255)
                                                                                  cttccgatctAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCgaga                <  1:1914228/91‑1 (MQ=255)
                                                                                   agagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                 <  1:1878527/98‑1 (MQ=255)
                                                                                   agagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaa               >  2:698646/1‑100 (MQ=255)
                                                                                    gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt              <  2:499519/100‑1 (MQ=255)
                                                                                    gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAAt              >  1:2633678/1‑100 (MQ=255)
                                                                                    gagCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCAACAACCACAAACCAGACAAATACCGAGCGAGAAt              >  1:2421826/1‑100 (MQ=255)
                                                                                     agCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATa             <  2:1761090/100‑1 (MQ=255)
                                                                                      gCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATAc            >  2:1450483/1‑100 (MQ=255)
                                                                                       ccGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACa           <  2:3296432/100‑1 (MQ=255)
                                                                                        cGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACat          <  1:1751750/100‑1 (MQ=255)
                                                                                        cGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACat          <  1:1205408/100‑1 (MQ=255)
                                                                                         gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata         >  1:113547/1‑100 (MQ=255)
                                                                                         gTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAACGAGAATACata         >  1:2533061/1‑100 (MQ=255)
                                                                                            aaaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGaagatcggaag      >  2:1914237/1‑91 (MQ=255)
                                                                                             aaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc     >  1:1693603/1‑100 (MQ=255)
                                                                                             aaaCGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATAcc     <  2:806314/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
TTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAA  >  NC_000913/3171765‑3171960

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