Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 332,895 A→G Q175Q (CAA→CAG pdeL → DNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913332,8950AG98.8% 286.4 / ‑5.1 82Q175Q (CAA→CAGpdeLDNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (44/37);  total (45/37)
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.99e-01

TCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATA  >  NC_000913/332759‑333031
                                                                                                                                        |                                                                                                                                         
tCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCtg                                                                                                                                        <  1:1488047/140‑1 (MQ=255)
          gACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGa                                                                                                                              <  2:1529281/140‑1 (MQ=255)
          gACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGa                                                                                                                              <  1:1612391/140‑1 (MQ=255)
                 gTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGa                                                                                                                       <  2:779115/140‑1 (MQ=255)
                 gTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGa                                                                                                                       <  1:871369/140‑1 (MQ=255)
                 gTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGa                                                                                                                       <  2:1551414/140‑1 (MQ=255)
                   aCTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTg                                                                                                                           >  2:210153/1‑134 (MQ=255)
                   aCTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTg                                                                                                                           <  1:210153/134‑1 (MQ=255)
                     tGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGa                                                                                                                       <  2:73542/136‑1 (MQ=255)
                     tGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGa                                                                                                                       >  1:73542/1‑136 (MQ=255)
                     tGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATAtt                                                                                                                   >  2:1445845/1‑140 (MQ=255)
                     tGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATAtt                                                                                                                   >  2:1506024/1‑140 (MQ=255)
                     tGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATAtt                                                                                                                   >  1:1092023/1‑140 (MQ=255)
                     tGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATAtt                                                                                                                   >  1:448746/1‑140 (MQ=255)
                     tGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATAtt                                                                                                                   >  1:1102666/1‑140 (MQ=255)
                         gggCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGc                                                                                                                          <  2:1092067/129‑1 (MQ=255)
                         gggCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGc                                                                                                                          >  1:1092067/1‑129 (MQ=255)
                         gggCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGAt                                                                                                              >  1:1354972/1‑140 (MQ=255)
                           gCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATg                                                                                                             <  2:837028/140‑1 (MQ=255)
                              gtgAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGAt                                                                                                              <  2:231572/136‑1 (MQ=255)
                              gtgAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGAt                                                                                                              >  1:231572/1‑136 (MQ=255)
                              gtgAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCg                                                                                                          >  2:1105583/1‑140 (MQ=255)
                               tgAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCgg                                                                                                         >  1:625247/1‑140 (MQ=255)
                                      cTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACa                                                                                                  >  1:935709/1‑140 (MQ=255)
                                          tCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTtg                                                                                              <  1:1105583/140‑1 (MQ=255)
                                           ccGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTtgc                                                                                             >  2:724986/1‑140 (MQ=255)
                                              cTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTtgctgc                                                                                          >  1:1033660/1‑140 (MQ=255)
                                                  gAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGa                                                                                      <  2:1102666/140‑1 (MQ=255)
                                                  gAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGa                                                                                      <  2:625247/140‑1 (MQ=255)
                                                    aCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGa                                                                                      <  2:789025/138‑1 (MQ=255)
                                                    aCATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGa                                                                                      >  1:789025/1‑138 (MQ=255)
                                                     cATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACaa                                                                                   <  2:1530200/140‑1 (MQ=255)
                                                          acaAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTcc                                                                              <  1:663353/140‑1 (MQ=255)
                                                           caAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCa                                                                             >  2:1440819/1‑140 (MQ=255)
                                                             aaCGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCata                                                                           >  1:1408164/1‑140 (MQ=255)
                                                             aaCGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCata                                                                           <  2:521147/140‑1 (MQ=255)
                                                             aaCGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTCCCata                                                                           <  2:1383650/140‑1 (MQ=255)
                                                                gggAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTg                                                                        >  1:614311/1‑140 (MQ=255)
                                                                  gAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGc                                                                      <  1:887975/140‑1 (MQ=255)
                                                                   aaTTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCa                                                                     >  2:819349/1‑140 (MQ=255)
                                                                    aTTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCAt                                                                    <  1:724986/140‑1 (MQ=255)
                                                                       aTCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCaa                                                                 <  2:312875/140‑1 (MQ=255)
                                                                             ccGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGtctc                                                           >  2:525145/1‑140 (MQ=255)
                                                                                gATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTcggc                                                        >  2:685978/1‑140 (MQ=255)
                                                                                    aGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGgt                                                    <  2:1146683/140‑1 (MQ=255)
                                                                                     gTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCggggt                                                   >  1:69552/1‑138 (MQ=255)
                                                                                         aTTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGttt                                               <  2:1204929/140‑1 (MQ=255)
                                                                                                tGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGgt                                                    <  2:254397/128‑1 (MQ=255)
                                                                                                tGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGgt                                                    >  1:254397/1‑128 (MQ=255)
                                                                                                tGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCa                                        >  1:1593066/1‑140 (MQ=255)
                                                                                                 ggCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAg                                       >  1:163309/1‑140 (MQ=255)
                                                                                                   cGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCg                                     <  2:1354972/140‑1 (MQ=255)
                                                                                                       gTCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGAt                                <  1:1506024/140‑1 (MQ=255)
                                                                                                        tCATCCGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGGtt                               >  2:882157/1‑140 (MQ=255)
                                                                                                            ccGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTGaa                           >  1:1075505/1‑140 (MQ=255)
                                                                                                            ccGGTCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGaa                           >  2:983470/1‑140 (MQ=255)
                                                                                                                tCTTATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGGTTTGaaaaaa                       >  2:890741/1‑139 (MQ=255)
                                                                                                                  ttATTGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGATGTTTTTTGTCAGCGGGATTTTGAAAAagg                      >  2:620013/1‑139 (MQ=255)
                                                                                                                     ttGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAGAgtgt                  >  2:1514223/1‑140 (MQ=255)
                                                                                                                     ttGTCATAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGCTGTTTTTTTGCAGCGGGGATTGCAAAagaaggt                  >  1:69232/1‑136 (MQ=255)
                                                                                                                          aTAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTTAAAAAAAGTGTCTGaa             >  2:353423/1‑140 (MQ=255)
                                                                                                                          aTAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAAAGGGTCTTaa             >  2:1096664/1‑140 (MQ=255)
                                                                                                                          aTAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGaa             >  1:1450020/1‑140 (MQ=255)
                                                                                                                          aTAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTTGAAAAGGGGGGCTGaa             >  2:291621/1‑140 (MQ=255)
                                                                                                                           tAATGACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAAc            <  1:525145/140‑1 (MQ=255)
                                                                                                                               gACCCGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTgg        <  2:163309/140‑1 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGGCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTGTTTGCAGCGGGGTTTGAAAAAGAGGGTCTGAACCAGGttt     >  2:1467026/1‑139 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGGTTTAAAAAAAAGGTCTGGAACTgggtt     >  1:829533/1‑137 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                      <  2:456776/123‑1 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                      >  2:1067806/1‑123 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                      >  1:456776/1‑123 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                      >  1:1624905/1‑123 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                      <  2:1624905/123‑1 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                      <  1:1067806/123‑1 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTta     <  1:353423/140‑1 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTta     <  1:1514223/140‑1 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTta     <  2:1412790/140‑1 (MQ=255)
                                                                                                                                  ccGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTta     >  1:1412790/1‑140 (MQ=255)
                                                                                                                                   cGCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTtaa    <  1:685978/140‑1 (MQ=255)
                                                                                                                                    gCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAAcc           >  1:897065/1‑132 (MQ=255)
                                                                                                                                    gCCAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAAcc           <  2:897065/132‑1 (MQ=255)
                                                                                                                                     ccAGCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTAGCAGCGGGATTTGTAAAAGAGCGTCTTCACCTTGTtaata  >  1:988536/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                         
TCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATA  >  NC_000913/332759‑333031

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