Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,217,796 Δ1 bp coding (469/1524 nt) pdeG → putative c‑di‑GMP phosphodiesterase PdeG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,217,7930G.99.7% 1536.5 / ‑2.6 294coding (466/1524 nt)pdeGputative c‑di‑GMP phosphodiesterase PdeG
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (155/138);  total (156/138)
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.92e-01

Only 100 of 323 total aligned reads displayed.

TCCAGCAATCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACGCTTAATGGAC  >  NC_000913/1217693‑1217889
                                                                                                    |                                                                                                
 ccAGCAATCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTgg                                                                                                 >  1:2049443/1‑101 (MQ=255)
 ccAGCAATCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTgg                                                                                                 <  1:2049128/101‑1 (MQ=255)
  cAGCAATCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTggg                                                                                                >  2:2049445/1‑101 (MQ=255)
  cAGCAATCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTg                                                                                                  >  1:2049447/1‑99 (MQ=255)
     cAATCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTgggcac                                                                                             <  2:2048977/101‑4 (MQ=255)
     cAATCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTgggcac                                                                                             <  2:2048922/101‑4 (MQ=255)
      aaTCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTgggcact                                                                                            >  2:2049451/1‑97 (MQ=255)
      aaTCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTgggcact                                                                                            >  1:2049454/1‑97 (MQ=255)
       aTCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACtt                                                                                          >  2:2049458/1‑101 (MQ=255)
        tCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTa                                                                                         >  1:2049459/1‑101 (MQ=255)
         cGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTaa                                                                                        <  2:2049103/101‑1 (MQ=255)
          gACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAAt                                                                                       <  2:2048844/101‑1 (MQ=255)
          gACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAAt                                                                                       >  1:2049462/1‑101 (MQ=255)
           acacCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATg                                                                                      >  1:2049466/1‑101 (MQ=255)
             acCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTa                                                                                    >  1:2049472/1‑101 (MQ=255)
               cGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATg                                                                                      <  1:2049343/97‑1 (MQ=255)
                gTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTAccc                                                                                 <  1:2048742/101‑1 (MQ=255)
                gTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTAccc                                                                                 >  2:2049477/1‑101 (MQ=255)
                    aTAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGa                                                                             >  1:2049480/1‑101 (MQ=255)
                     tAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGaa                                                                            <  1:2049070/101‑1 (MQ=255)
                       gATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAg                                                                          >  1:2049484/1‑101 (MQ=255)
                       gATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAg                                                                          >  2:2049485/1‑101 (MQ=255)
                        aTTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAgg                                                                         <  1:2048557/101‑1 (MQ=255)
                          tACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAgggg                                                                       >  2:2049489/1‑101 (MQ=255)
                           aCCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGt                                                                      <  1:2048889/101‑1 (MQ=255)
                             cTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTaa                                                                    <  2:2049078/101‑1 (MQ=255)
                             cTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTaa                                                                    >  1:2049493/1‑101 (MQ=255)
                               aTCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAgg                                                                  <  2:2048948/101‑1 (MQ=255)
                               aTCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAgg                                                                  >  1:2049497/1‑101 (MQ=255)
                                tCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGgt                                                                 >  2:2049499/1‑101 (MQ=255)
                                   tGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGgtatg                                                              >  1:2049502/1‑101 (MQ=255)
                                      tctATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTg                                                                                                  >  1:2049508/1‑63 (MQ=255)
                                      tctATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGgtatgtat                                                           <  1:2049327/101‑1 (MQ=255)
                                      tctATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGgtatgtat                                                           >  2:2049507/1‑101 (MQ=255)
                                        tATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTa                                                         <  1:2049235/101‑1 (MQ=255)
                                        tATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTa                                                         >  1:2049515/1‑101 (MQ=255)
                                        tATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTa                                                         >  1:2049512/1‑101 (MQ=255)
                                         aTCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTAc                                                        <  1:2048933/101‑1 (MQ=255)
                                          tCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTAcg                                                       <  1:2048997/101‑1 (MQ=255)
                                           cAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTAcgc                                                      >  1:2049519/1‑101 (MQ=255)
                                           cAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTAcgc                                                      <  1:2049019/101‑1 (MQ=255)
                                            aGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTAcgcg                                                     >  2:2049524/1‑101 (MQ=255)
                                            aGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTAcgcg                                                     >  1:2049523/1‑101 (MQ=255)
                                             gAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGt                                                    <  1:2049213/101‑1 (MQ=255)
                                              aaCCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTAcg                                                       >  1:2049528/1‑97 (MQ=255)
                                              aaCCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTg                                                   <  1:2049048/101‑1 (MQ=255)
                                                ccAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGc                                                 >  1:2049531/1‑101 (MQ=255)
                                                ccAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGc                                                 <  2:2049117/101‑1 (MQ=255)
                                                 cAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCg                                                <  2:2049086/101‑1 (MQ=255)
                                                 cAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCg                                                >  1:2049536/1‑101 (MQ=255)
                                                  aaTTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACtt                                                                                          <  2:2048871/58‑1 (MQ=255)
                                                  aaTTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCgg                                               >  2:2049537/1‑101 (MQ=255)
                                                  aaTTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCgg                                               <  1:2049106/101‑1 (MQ=255)
                                                   aTTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGa                                              <  2:2049125/101‑1 (MQ=255)
                                                   aTTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGa                                              >  1:2049540/1‑101 (MQ=255)
                                                   aTTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGa                                              <  2:2049046/101‑1 (MQ=255)
                                                    tttGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGAt                                             <  2:2049266/101‑1 (MQ=255)
                                                     ttGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCgg                                               >  2:2049546/1‑98 (MQ=255)
                                                     ttGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGAt                                             <  1:2048981/100‑1 (MQ=255)
                                                     ttGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATg                                            <  2:2049000/101‑1 (MQ=255)
                                                     ttGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATg                                            >  1:2049545/1‑101 (MQ=255)
                                                        cAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACa                                         >  1:2049553/1‑101 (MQ=255)
                                                        cAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACa                                         >  2:2049551/1‑101 (MQ=255)
                                                         aGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGaa                                                                            >  1:2049559/1‑65 (MQ=255)
                                                           aTACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACatta                                      >  2:2049561/1‑101 (MQ=255)
                                                            tACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACattat                                     <  2:2049063/101‑1 (MQ=255)
                                                              cgcgCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTg                                   >  1:2049563/1‑101 (MQ=255)
                                                              cgcgCATTTTGGTTACGATAAGAGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTg                                   >  1:2049566/1‑101 (MQ=255)
                                                                cgcATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGa                                           <  2:2049183/91‑1 (MQ=255)
                                                                cgcATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGa                                 >  2:2049570/1‑101 (MQ=255)
                                                                 gcATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGAc                                <  1:2049140/101‑1 (MQ=255)
                                                                  cATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTa                                                                     >  2:2049574/1‑63 (MQ=255)
                                                                   aTTTTGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCa                              >  1:2049575/1‑101 (MQ=255)
                                                                    ttttGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCaa                             >  2:2049579/1‑101 (MQ=255)
                                                                     tttGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAAc                            >  1:2049581/1‑101 (MQ=255)
                                                                     tttGGTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAAc                            <  1:2049114/101‑1 (MQ=255)
                                                                         gTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGgt                        <  1:2049037/101‑1 (MQ=255)
                                                                         gTTACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGgt                        <  1:2049053/101‑1 (MQ=255)
                                                                          ttACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGgtg                       <  1:2049014/101‑1 (MQ=255)
                                                                          ttACGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGgtg                       >  2:2049588/1‑101 (MQ=255)
                                                                            aCGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGgtgat                     >  1:2049590/1‑101 (MQ=255)
                                                                             cGATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGgtgatg                    >  1:2049593/1‑101 (MQ=255)
                                                                              gATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAgg                                                                         >  2:2049596/1‑47 (MQ=255)
                                                                              gATAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGgtgatgt                   <  2:2049196/101‑1 (MQ=255)
                                                                               aTAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGgtgatgtg                  >  1:2049598/1‑101 (MQ=255)
                                                                                tAAGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGgtgatgtga                 <  2:2049469/101‑1 (MQ=255)
                                                                                 aaGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGgtgatgtgat                <  2:2049281/101‑1 (MQ=255)
                                                                                 aaGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGgtgatgtgat                <  2:2049147/101‑1 (MQ=255)
                                                                                  aGCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGtgatgtgatg               <  2:2049030/101‑1 (MQ=255)
                                                                                   gCGATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACggg                         >  2:2049602/1‑90 (MQ=255)
                                                                                     gATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACg            >  1:2049604/1‑101 (MQ=255)
                                                                                     gATCAGCATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACg            <  2:2048741/101‑1 (MQ=255)
                                                                                           cATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACGCTTAAt      >  2:2049609/1‑101 (MQ=255)
                                                                                           cATATTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACGCTTAAt      <  1:2049179/101‑1 (MQ=255)
                                                                                            atatTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACGCTTAATg     <  1:2049084/101‑1 (MQ=255)
                                                                                            atatTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACGCTTAATg     <  2:2049333/101‑1 (MQ=255)
                                                                                            atatTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACGCTTAATg     >  1:2049613/1‑101 (MQ=255)
                                                                                             tatTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACGCTTAATgg    <  1:2049304/101‑1 (MQ=255)
                                                                                              atTCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACGCTTAATGGa   <  1:2048955/101‑1 (MQ=255)
                                                                                               ttCGTGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACGCTTAATGGAc  <  1:2049248/101‑1 (MQ=255)
                                                                                                    |                                                                                                
TCCAGCAATCGACACCGTTAATAGATTACCTATCCTGCTCTATCAGAACCAATTTGCAGATACGCGCATTTTGGTTACGATAAGCGATCAGCATATTCGTGGGGCACTTAATGTACCCTTGAAAGGGGTAAGGTATGTATTACGCGTGGCGGATGACATTATTGGACCAACGGGTGATGTGATGACGCTTAATGGAC  >  NC_000913/1217693‑1217889

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