Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,323,955 T→G 42.5% N133T (AAC→ACC)  phnD ← phosphonate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,323,9550TG42.5% 109.7 / 274.6 205N133T (AAC→ACC) phnDphosphonate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (59/57);  new base G (47/39);  total (107/98)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

GCGAAGACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCATTGGCGCG  >  NC_000913/4323856‑4324052
                                                                                                   |                                                                                                 
gCGAAGACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGgt                                                                                                   <  2:2771467/100‑1 (MQ=255)
gCGAAGACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGgt                                                                                                   >  2:1120060/1‑100 (MQ=255)
 cGAAGACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGgtt                                                                                                  >  2:2570609/1‑100 (MQ=255)
 cGAAGACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGgtt                                                                                                  >  2:243653/1‑100 (MQ=255)
  gAAGACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGgttg                                                                                                 <  2:2477701/100‑1 (MQ=255)
  gAAGACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGggtg                                                                                                 >  2:4074047/1‑100 (MQ=255)
   aaGACGTAGTAGCCGGGGACGAGGACGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGgttgt                                                                                                <  2:1070910/100‑1 (MQ=255)
    aGACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGgttgtt                                                                                               >  1:2769509/1‑100 (MQ=255)
    aGACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGgttgtt                                                                                               <  1:2580128/100‑1 (MQ=255)
     gACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGggtgtt                                                                                               <  1:1901593/99‑1 (MQ=255)
     gACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGggtgt                                                                                                <  2:3900191/98‑1 (MQ=255)
       cGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGAt                                                                                            >  2:3459498/1‑100 (MQ=255)
       cGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGAt                                                                                            >  2:3156382/1‑100 (MQ=255)
        gtagtaGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATc                                                                                           >  2:642876/1‑100 (MQ=255)
        gtagtaGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATc                                                                                           <  1:2867665/100‑1 (MQ=255)
            tagCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGAc                                                                                       <  1:1239071/100‑1 (MQ=255)
               ccGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGt                                                                                    >  1:116351/1‑100 (MQ=255)
                cGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTc                                                                                   >  1:19284/1‑100 (MQ=255)
                cGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTc                                                                                   >  2:542141/1‑100 (MQ=255)
                 ggggACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCt                                                                                  <  2:2495775/100‑1 (MQ=255)
                    gACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTttg                                                                               <  1:2032281/100‑1 (MQ=255)
                     aCGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTttgt                                                                              <  1:1774681/100‑1 (MQ=255)
                      cGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTttgtt                                                                             >  2:1896070/1‑100 (MQ=255)
                       gAGGAAGCCAGAGGTGGAGTTATGCTCGCCAGTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTttgt                                                                              <  1:3550529/98‑1 (MQ=255)
                       gAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGGTCGTTCAGGGCGTTGATCGGACTGTCTttgttg                                                                            <  1:2487122/100‑1 (MQ=255)
                       gAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTttgttg                                                                            <  2:1529670/100‑1 (MQ=255)
                         ngAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGAc                                                                          >  1:3745990/2‑100 (MQ=255)
                         ggAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGAc                                                                          <  2:4113878/100‑1 (MQ=255)
                         ggAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGAc                                                                          <  2:1097192/100‑1 (MQ=255)
                         ggAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGAc                                                                          >  2:2285436/1‑100 (MQ=255)
                         ggAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGAc                                                                          <  2:4167470/100‑1 (MQ=255)
                           aaGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGa                                                                        >  1:3785781/1‑100 (MQ=255)
                            aGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGAt                                                                       <  2:84556/100‑1 (MQ=255)
                            aGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGAt                                                                       <  1:2174974/100‑1 (MQ=255)
                             gCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATc                                                                      <  2:1342530/100‑1 (MQ=255)
                             gCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATc                                                                      <  2:3884506/100‑1 (MQ=255)
                              ctagagGTCGAGTTACGATCGCCATTGCCAAAGTTGAGATCTTTCCGCTTCTCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGAt                                                                       <  1:3968052/96‑1 (MQ=255)
                              ccAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCa                                                                     >  1:3476715/1‑100 (MQ=255)
                              ccAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCa                                                                     <  1:3075946/100‑1 (MQ=255)
                              ccAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCa                                                                     >  2:2351983/1‑100 (MQ=255)
                              ccAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCa                                                                     <  1:2096707/100‑1 (MQ=255)
                               cAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCaa                                                                    <  2:3115800/100‑1 (MQ=255)
                               cAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCaa                                                                    >  1:2420338/1‑100 (MQ=255)
                               cAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCaa                                                                    >  2:1064149/1‑100 (MQ=255)
                               cAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCaa                                                                    >  2:1202263/1‑100 (MQ=255)
                                 gagGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAaca                                                                  <  2:1687894/100‑1 (MQ=255)
                                  agGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCTGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAacac                                                                 <  1:2362533/100‑1 (MQ=255)
                                  agGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAacac                                                                 <  1:1742585/100‑1 (MQ=255)
                                   ggTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACg                                                                >  1:2940405/1‑100 (MQ=255)
                                   ggTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACg                                                                >  2:1247874/1‑100 (MQ=255)
                                       gAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCGACACGCTcc                                                            <  2:1301386/100‑1 (MQ=255)
                                        aGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCa                                                           >  2:3202373/1‑100 (MQ=255)
                                        aGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCa                                                           >  2:604582/1‑100 (MQ=255)
                                        aGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCa                                                           >  2:235299/1‑100 (MQ=255)
                                        aGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCa                                                           >  2:1366231/1‑100 (MQ=255)
                                        aGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTcc                                                            <  1:602653/99‑1 (MQ=255)
                                         gTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAg                                                          <  2:3567751/100‑1 (MQ=255)
                                           tAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTa                                                        >  2:3072627/1‑100 (MQ=255)
                                            aGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCACCAAGCTCCAGaaa                                                       >  2:3388133/1‑100 (MQ=255)
                                              gATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAAcc                                                     <  1:4081376/100‑1 (MQ=255)
                                               aTCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAAccc                                                    >  2:1608982/1‑100 (MQ=255)
                                                    cATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGa                                               >  2:3288680/1‑100 (MQ=255)
                                                    cATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTCACGATCAACACGCTCCAGTAACCCGGCg                                                <  1:336342/99‑1 (MQ=255)
                                                      atctcAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGatc                                             <  2:2839428/96‑1 (MQ=255)
                                                        gCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTTATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTACCCCGGCGatcca                                           >  1:324186/1‑100 (MQ=255)
                                                        gCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGatcca                                           >  1:2679354/1‑100 (MQ=255)
                                                        gCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGatcca                                           >  2:4023347/1‑100 (MQ=255)
                                                        gCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGatcca                                           >  1:1221599/1‑100 (MQ=255)
                                                          cAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGatccatc                                         >  2:744606/1‑100 (MQ=255)
                                                          cAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGatcagat                                         >  1:2839417/1‑96 (MQ=255)
                                                          cAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGatccatc                                         >  2:3938234/1‑100 (MQ=255)
                                                          cAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGatccatc                                         >  1:1261645/1‑100 (MQ=255)
                                                            aaGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATccg                                       <  1:579587/100‑1 (MQ=255)
                                                               gTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTGGTTGAGGGGACTGTCTTTGTTGACGATCAACATGGTCCAGTAACCAGGCGATCCATccgccg                                    >  1:2293813/1‑100 (MQ=255)
                                                                tGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCcgccgc                                   <  1:544721/100‑1 (MQ=255)
                                                                   gaTCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGAc                                <  1:1026364/100‑1 (MQ=255)
                                                                      cTTTCCGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGt                             >  1:14449/1‑100 (MQ=255)
                                                                       tttCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTc                            <  1:1881848/100‑1 (MQ=255)
                                                                        ttCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCt                           <  1:795002/100‑1 (MQ=255)
                                                                          ccGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTgg                         >  2:2917174/1‑100 (MQ=255)
                                                                           cGCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTggg                        <  1:2351977/100‑1 (MQ=255)
                                                                            gCTTCGCCAGCAGATCGTTCAGGGTGTTGATCGGACTGCCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGGCGCGACCGTCTGGGc                       >  1:2638307/1‑100 (MQ=255)
                                                                                agcggaaagATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAg                   <  2:2979637/93‑1 (MQ=255)
                                                                                 gCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGa                  >  2:157290/1‑100 (MQ=255)
                                                                                 gCCAGCAGATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGa                  >  2:1268174/1‑100 (MQ=255)
                                                                                   cagcagATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGAcc                <  2:2050950/100‑1 (MQ=255)
                                                                                     gcagATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTg              <  2:2225591/100‑1 (MQ=255)
                                                                                      cagATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTgg             >  1:1498743/1‑100 (MQ=255)
                                                                                       agATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGc            >  2:2784992/1‑100 (MQ=255)
                                                                                        gATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGcc           >  2:3869128/1‑100 (MQ=255)
                                                                                        gATCGTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGcc           >  1:3284220/1‑100 (MQ=255)
                                                                                            gTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCATTg       >  1:3156051/1‑100 (MQ=255)
                                                                                            gTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCATTg       >  1:1526073/1‑100 (MQ=255)
                                                                                            gTTCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCATTg       >  2:3563744/1‑100 (MQ=255)
                                                                                             ttCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCATTgg      <  1:3302127/100‑1 (MQ=255)
                                                                                              tCAGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCATTGgc     >  2:2123678/1‑100 (MQ=255)
                                                                                                aGGTTGTTGATCGGACCGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCATTGgcgc   <  2:3659367/100‑1 (MQ=255)
                                                                                                aGGGTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCATTGgcgc   <  1:4166165/100‑1 (MQ=255)
                                                                                                 ggTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCATTGGcgcg  >  2:221251/1‑100 (MQ=255)
                                                                                                 ggTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCATTGGcgcg  <  2:3809411/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GCGAAGACGTAGTAGCCGGGGACGAGGAAGCCAGAGGTGGAGTTAGGATCGCCATTGCCAAAGGTGAGATCTTTCCGCTTCGCCAGCAGATCGTTCAGGTTGTTGATCGGACTGTCTTTGTTGACGATCAACACGCTCCAGTAACCCGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCATTGGCGCG  >  NC_000913/4323856‑4324052

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