Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,027,073 T→C E315E (GAA→GAG dgcQ ← putative diguanylate cyclase DgcQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,027,0730TC98.3% 991.2 / 3.3 241E315E (GAA→GAGdgcQputative diguanylate cyclase DgcQ
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (102/134);  total (103/138)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 241 total aligned reads displayed.

CTTAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGT  >  NC_000913/2026975‑2027170
                                                                                                  |                                                                                                 
cTTAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTc                                                                                                 >  1:2949493/1‑101 (MQ=255)
 ttAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTcc                                                                                                >  2:2949495/1‑101 (MQ=255)
 ttAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTcc                                                                                                <  2:2949089/101‑1 (MQ=255)
  tAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCa                                                                                               <  2:2949182/101‑1 (MQ=255)
  tAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCa                                                                                               >  1:2949496/1‑101 (MQ=255)
    aCGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCAtt                                                                                             <  2:2949172/101‑1 (MQ=255)
      gTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGc                                                                                           >  1:2949499/1‑101 (MQ=255)
      gTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCAAGCTCCATTGc                                                                                           <  2:2949221/101‑1 (MQ=255)
          ggACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTgc                                                                                       <  2:2949206/101‑1 (MQ=255)
              acGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCa                                                                                   <  1:2949038/101‑1 (MQ=255)
              acGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCa                                                                                   >  1:2949504/1‑101 (MQ=255)
               cGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAg                                                                                  <  1:2949111/101‑1 (MQ=255)
                gCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGc                                                                                 <  1:2949170/101‑1 (MQ=255)
                    cAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAgg                                                                             <  2:2949051/101‑1 (MQ=255)
                       cacaCCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCa                                                                                >  1:2949510/1‑95 (MQ=255)
                       cacaCCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCa                                                                          <  2:2949276/101‑1 (MQ=255)
                       cacaCCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCa                                                                          >  2:2949509/1‑101 (MQ=255)
                        acacCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAg                                                                         <  2:2949128/101‑1 (MQ=255)
                         cacCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGt                                                                        <  1:2949084/101‑1 (MQ=255)
                          acCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGtt                                                                       <  1:2949238/101‑1 (MQ=255)
                          acCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGtt                                                                       >  1:2949513/1‑101 (MQ=255)
                           ccGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTc                                                                      <  1:2948948/101‑1 (MQ=255)
                           ccGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTc                                                                      <  1:2949196/101‑1 (MQ=255)
                           ccGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTc                                                                      <  2:2949141/101‑1 (MQ=255)
                            cGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCa                                                                     <  1:2949218/101‑1 (MQ=255)
                              tCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACg                                                                   >  2:2949515/1‑101 (MQ=255)
                               cGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGa                                                                  <  1:2949341/101‑1 (MQ=255)
                                gAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAg                                                                 >  1:2949519/1‑101 (MQ=255)
                                gAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAg                                                                 >  1:2949517/1‑101 (MQ=255)
                                 aaaTGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAgg                                                                <  1:2949019/101‑1 (MQ=255)
                                 aaaTGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAgg                                                                >  1:2949522/1‑101 (MQ=255)
                                   aTGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGAt                                                              >  2:2949524/1‑101 (MQ=255)
                                   aTGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGAt                                                              <  2:2949134/101‑1 (MQ=255)
                                      gTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCaa                                                           <  1:2949265/101‑1 (MQ=255)
                                       tCCAGACGTTCCCAGCTAACATAGCGACTGTCTATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCaaa                                                          >  1:2949528/1‑101 (MQ=255)
                                       tCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCaaa                                                          <  1:2948887/101‑1 (MQ=255)
                                         cAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCaaaaa                                                        >  2:2949532/1‑101 (MQ=255)
                                          aGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAat                                                       >  2:2949534/1‑101 (MQ=255)
                                           gACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAata                                                      <  1:2949201/101‑1 (MQ=255)
                                             cGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATAtt                                                    <  2:2949116/101‑1 (MQ=255)
                                              gTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTc                                                   >  1:2949537/1‑101 (MQ=255)
                                               ttCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGa                                                                  <  2:2949279/85‑1 (MQ=255)
                                                tCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTccc                                                 >  1:2949539/1‑101 (MQ=255)
                                                 cccAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCt                                                >  2:2949544/1‑101 (MQ=255)
                                                 cccAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCt                                                <  2:2949208/101‑1 (MQ=255)
                                                 cccAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCt                                                >  2:2949541/1‑101 (MQ=255)
                                                  ccAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTg                                               <  2:2949271/101‑1 (MQ=255)
                                                  ccAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTg                                               >  1:2949545/1‑101 (MQ=255)
                                                  ccAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTg                                               <  1:2948985/101‑1 (MQ=255)
                                                    aGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGtt                                             <  2:2949299/101‑1 (MQ=255)
                                                      cTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTgg                                           <  2:2948881/101‑1 (MQ=255)
                                                      cTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTgg                                           <  2:2949155/101‑1 (MQ=255)
                                                      cTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTgg                                           >  1:2949548/1‑101 (MQ=255)
                                                      cTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTgg                                           >  2:2949549/1‑101 (MQ=255)
                                                         aCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATg                                        <  2:2949252/101‑1 (MQ=255)
                                                         aCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCTCGAGGATCAAAAATATTCCCTGTTGGATg                                        <  1:2949240/101‑1 (MQ=255)
                                                         aCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATg                                        >  1:2949553/1‑101 (MQ=255)
                                                         aCATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATg                                        <  2:2949005/101‑1 (MQ=255)
                                                          cATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGATCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATgg                                       <  1:2949152/101‑1 (MQ=255)
                                                           aTAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGt                                      >  2:2949555/1‑101 (MQ=255)
                                                            tAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTc                                     >  1:2949558/1‑101 (MQ=255)
                                                             aGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCa                                    <  2:2949071/101‑1 (MQ=255)
                                                               cGACTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAgg                                  >  2:2949560/1‑101 (MQ=255)
                                                                 aCTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGAt                                <  1:2949231/101‑1 (MQ=255)
                                                                 aCTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGAt                                >  2:2949563/1‑101 (MQ=255)
                                                                  cTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGAtt                               <  2:2949501/101‑1 (MQ=255)
                                                                  cTGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGAtt                               >  1:2949564/1‑101 (MQ=255)
                                                                   tGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTg                              >  2:2949567/1‑101 (MQ=255)
                                                                   tGTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTg                              <  1:2949418/101‑1 (MQ=255)
                                                                    gTCCATACGAATGCCGCCCCGCGTGTCATGCTCCATTCCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAgg                                  >  1:2949569/1‑96 (MQ=255)
                                                                      ccATACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGaa                           <  1:2949162/101‑1 (MQ=255)
                                                                        aTACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAaga                         <  1:2949096/101‑1 (MQ=255)
                                                                        aTACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAaga                         <  2:2949244/101‑1 (MQ=255)
                                                                         tACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAagag                        <  1:2949131/101‑1 (MQ=255)
                                                                         tACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAagag                        <  2:2949388/101‑1 (MQ=255)
                                                                         tACGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAagag                        <  1:2949321/101‑1 (MQ=255)
                                                                          aCGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAgg                       >  1:2949572/1‑101 (MQ=255)
                                                                          aCGAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAgg                       <  1:2949375/101‑1 (MQ=255)
                                                                            gAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTc                     >  2:2949575/1‑101 (MQ=255)
                                                                            gAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTc                     <  1:2949062/101‑1 (MQ=255)
                                                                            gAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTc                     <  1:2949147/101‑1 (MQ=255)
                                                                            gAATGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTc                     >  1:2949574/1‑101 (MQ=255)
                                                                             aaTGCCGCCCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCa                    <  1:2949357/101‑1 (MQ=255)
                                                                                gccgccCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCaaaa                 <  1:2949159/101‑1 (MQ=255)
                                                                                gccgccCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCaaaa                 >  1:2949579/1‑101 (MQ=255)
                                                                                 ccgccCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTc                                     >  2:2949582/1‑80 (MQ=255)
                                                                                  cgccCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTc                     >  1:2949583/1‑95 (MQ=255)
                                                                                   gccCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAtc              >  2:2949585/1‑101 (MQ=255)
                                                                                   gccCCGCGTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAAtc              <  2:2949188/101‑1 (MQ=255)
                                                                                       cgcgTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAg          <  1:2949444/101‑1 (MQ=255)
                                                                                         cgTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGTTTGGAAGAGGTCAAAAATCTCAGCt        >  1:2949591/1‑101 (MQ=255)
                                                                                         cgTGTCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCt        >  2:2949589/1‑101 (MQ=255)
                                                                                          gtgtCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCtt       <  2:2949292/101‑1 (MQ=255)
                                                                                          gtgtCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCtt       <  1:2949258/101‑1 (MQ=255)
                                                                                           tgtCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTa      <  2:2949334/101‑1 (MQ=255)
                                                                                           tgtCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTa      <  1:2949052/101‑1 (MQ=255)
                                                                                             tCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACt    <  1:2949330/101‑1 (MQ=255)
                                                                                             tCATGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACt    >  1:2949594/1‑101 (MQ=255)
                                                                                              cATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTg   <  1:2949295/101‑1 (MQ=255)
                                                                                               aTGCTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGt  >  2:2949595/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CTTAACGTATGGACACGCACCAGCACACCGTCGAAATGGTCCAGACGTTCCCAGCTAACATAGCGACTGTCCATACGAATGCCGCCCCGCGTGTCATGTTCCATTGCCTGCGCCAGCAAGGCCAGTTCACGAGGATCAAAAATATTCCCTGTTGGATGGTCAGGATTGGAAGAGGTCAAAAATCTCAGCTTACTGT  >  NC_000913/2026975‑2027170

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: