Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,445,885 A→G I115V (ATT→GTT)  feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,445,8850AG98.9% 284.3 / ‑4.4 92I115V (ATT→GTT) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (52/39);  total (52/40)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

CGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTAAAACGCTGG  >  NC_000913/1445795‑1445972
                                                                                          |                                                                                       
cGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGtt                                                                                       <  1:2375787/93‑1 (MQ=255)
    ttCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTccc                                                                                   >  2:4263228/1‑93 (MQ=255)
    ttCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTccc                                                                                   >  1:2288079/1‑93 (MQ=255)
    ttCTGCTACCTTTTGCGGATTTGGTGGAACAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTccc                                                                                   <  2:3624693/93‑1 (MQ=255)
     tCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCg                                                                                  <  2:1460003/93‑1 (MQ=255)
     tCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCg                                                                                  <  2:4105230/93‑1 (MQ=255)
      cTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGt                                                                                 >  2:2249004/1‑93 (MQ=255)
      cTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGt                                                                                 >  1:1436462/1‑93 (MQ=255)
        gCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGc                                                                               >  1:2620683/1‑93 (MQ=255)
         cTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCt                                                                              >  2:2300138/1‑93 (MQ=255)
          tACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCtt                                                                             >  2:3720881/1‑93 (MQ=255)
          tACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCtt                                                                             >  1:702845/1‑93 (MQ=255)
          tACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCtt                                                                             >  2:3787980/1‑93 (MQ=255)
          tACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGATGGCGCAACTGGCAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCtt                                                                             >  2:2016485/1‑93 (MQ=255)
           aCGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCttt                                                                            <  1:1146744/93‑1 (MQ=255)
           aCGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCttt                                                                            >  2:477949/1‑93 (MQ=255)
           aCGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCttt                                                                            >  1:1312144/1‑93 (MQ=255)
             gTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCCGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTg                                                                          >  1:2228967/1‑93 (MQ=255)
                ttGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAg                                                                       <  2:2817386/93‑1 (MQ=255)
                 tGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAACCCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGt                                                                      <  1:2094554/93‑1 (MQ=255)
                 tGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGt                                                                      >  1:1125713/1‑93 (MQ=255)
                 tGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGt                                                                      <  1:2562918/93‑1 (MQ=255)
                   cTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTgg                                                                    <  2:592508/93‑1 (MQ=255)
                    cgATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTggg                                                                   <  2:3527986/92‑1 (MQ=255)
                     gATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGc                                                                  <  1:3882950/93‑1 (MQ=255)
                     gATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGc                                                                  <  1:2918842/93‑1 (MQ=255)
                     gATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGc                                                                  <  1:2113327/93‑1 (MQ=255)
                       tCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTg                                                                >  2:2428732/1‑93 (MQ=255)
                       tCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTg                                                                >  2:2593626/1‑93 (MQ=255)
                        cTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGt                                                               <  2:1912221/93‑1 (MQ=255)
                        cTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGt                                                               <  1:1521702/93‑1 (MQ=255)
                           gtggAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACg                                                            >  1:53028/1‑93 (MQ=255)
                            tggAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGc                                                           >  2:602251/1‑93 (MQ=255)
                               agcagcACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGa                                                        <  2:3291218/93‑1 (MQ=255)
                                gcagcaCAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGaa                                                       <  2:3981522/93‑1 (MQ=255)
                                   gcaCAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTg                                                    >  1:1419236/1‑93 (MQ=255)
                                    cacaGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTgg                                                   <  2:4627046/93‑1 (MQ=255)
                                     acaGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCAGTTTCCCGTGCTTTTGAAGTGGGCTGGACGCTGAACTGGa                                                  >  1:357562/1‑93 (MQ=255)
                                      caGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGAt                                                 >  1:2478566/1‑93 (MQ=255)
                                        gTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGc                                               >  2:4442659/1‑93 (MQ=255)
                                        gTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGc                                               >  2:2943731/1‑93 (MQ=255)
                                        gTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGc                                               >  2:2639481/1‑93 (MQ=255)
                                          gaggagCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGt                                             <  2:706701/93‑1 (MQ=255)
                                            ggagCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTta                                           <  2:1307349/93‑1 (MQ=255)
                                            ggagCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTta                                           <  2:2568948/93‑1 (MQ=255)
                                               gCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATAc                                        <  1:2374779/93‑1 (MQ=255)
                                               gCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATAc                                        >  1:246044/1‑93 (MQ=255)
                                                cTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATAcc                                       >  2:3421348/1‑93 (MQ=255)
                                                cTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATAcc                                       >  1:371384/1‑93 (MQ=255)
                                                 tGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATAccg                                      >  2:1889699/1‑93 (MQ=255)
                                                    cgcAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATAccgccg                                   >  1:4521406/1‑93 (MQ=255)
                                                     gcAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCgg                                  <  2:316302/93‑1 (MQ=255)
                                                     gcAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCgg                                  >  1:2232401/1‑93 (MQ=255)
                                                      cAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCggg                                 <  2:658568/93‑1 (MQ=255)
                                                       aaCTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGt                                <  1:519081/93‑1 (MQ=255)
                                                       aaCTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGt                                <  2:4622193/93‑1 (MQ=255)
                                                        aCTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTTTACGCTGAACTGGATGCGTTATACCGCCGGGtt                               >  1:1421035/1‑93 (MQ=255)
                                                         cTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTa                              <  1:681158/93‑1 (MQ=255)
                                                         cTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGCTa                              >  2:1869573/1‑93 (MQ=255)
                                                          tGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTaa                             >  2:829549/1‑93 (MQ=255)
                                                          tGGAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTaa                             >  2:702559/1‑93 (MQ=255)
                                                           ggAAACCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAAc                            <  2:1332753/93‑1 (MQ=255)
                                                               aCCCTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGAcc                        >  1:2322829/1‑93 (MQ=255)
                                                                 ccTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCaa                      <  2:1604943/93‑1 (MQ=255)
                                                                 ccTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCaa                      >  2:1825107/1‑93 (MQ=255)
                                                                 ccTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAAAGACCaa                      >  2:969609/1‑93 (MQ=255)
                                                                  cTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCaaa                     >  1:817636/1‑93 (MQ=255)
                                                                  cTGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCaaa                     >  1:4673321/1‑93 (MQ=255)
                                                                   tGGAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCaaaa                    <  2:1285900/93‑1 (MQ=255)
                                                                    ggAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAAt                   >  1:1736801/1‑93 (MQ=255)
                                                                    ggAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAAt                   >  1:2800881/1‑93 (MQ=255)
                                                                     gAGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATc                  <  1:1570418/93‑1 (MQ=255)
                                                                      aGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATcg                 <  2:4497962/93‑1 (MQ=255)
                                                                      aGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATcg                 >  1:2439553/1‑93 (MQ=255)
                                                                      aGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATcg                 >  1:2823501/1‑93 (MQ=255)
                                                                      aGCAAGGCAAGTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATcg                 >  2:4493691/1‑93 (MQ=255)
                                                                         aaggcaagTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGTCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCgg              <  2:501442/93‑1 (MQ=255)
                                                                         aaggcaagTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCgg              >  1:3807507/1‑93 (MQ=255)
                                                                         aaggcaagTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCgg              >  2:762264/1‑93 (MQ=255)
                                                                          aggcaagTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCggg             >  2:4480385/1‑93 (MQ=255)
                                                                          aggcaagTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCggg             <  2:346446/93‑1 (MQ=255)
                                                                          aggcaagTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCggg             >  2:2584092/1‑93 (MQ=255)
                                                                             caagTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTaa          <  1:3635050/93‑1 (MQ=255)
                                                                             caagTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTaa          >  1:1954829/1‑93 (MQ=255)
                                                                              aagTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTaaa         >  1:2457611/1‑93 (MQ=255)
                                                                                gTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTAAAAc       <  1:3045679/93‑1 (MQ=255)
                                                                                gTCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTAAAAc       >  1:2845020/1‑93 (MQ=255)
                                                                                 tCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTAAAACg      <  1:3847986/93‑1 (MQ=255)
                                                                                 tCAATTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTAAAACg      >  2:3678743/1‑93 (MQ=255)
                                                                                  cAATTGCTGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTAAAACGc     <  2:1555943/93‑1 (MQ=255)
                                                                                   aaTTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCTAAATCGCGGGTAAAACGCt    <  2:4525812/93‑1 (MQ=255)
                                                                                   aaTTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTAAAACGCt    <  2:553043/93‑1 (MQ=255)
                                                                                   aaTTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGGAAAACGCt    >  2:4171767/1‑93 (MQ=255)
                                                                                    aTTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTAAAACGCTg   <  1:487713/93‑1 (MQ=255)
                                                                                    aTTGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGGTAACGACCAAAATCGCGGGTAAAACGCTg   >  2:919943/1‑93 (MQ=255)
                                                                                     ttGCCGTTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTAAAACGCTgg  >  2:2102594/1‑93 (MQ=255)
                                                                                          |                                                                                       
CGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCCGGGTTAACGACCAAAATCGCGGGTAAAACGCTGG  >  NC_000913/1445795‑1445972

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