Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,448,947 A→C 26.4% T477P (ACG→CCG)  feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,448,9470AC26.4% 138.0 / 60.5 87T477P (ACG→CCG) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base C (12/11);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGC  >  NC_000913/1448850‑1449046
                                                                                                 |                                                                                                   
ggAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACg                                                                                                   >  1:3704232/1‑100 (MQ=255)
ggAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACg                                                                                                   >  2:1700882/1‑100 (MQ=255)
   atatAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGc                                                                                                <  2:1542280/100‑1 (MQ=255)
   atatAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGc                                                                                                <  2:277418/100‑1 (MQ=255)
   atatAGCAATCACTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGc                                                                                                <  2:3196819/100‑1 (MQ=255)
    tataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGcc                                                                                               >  1:398414/1‑100 (MQ=255)
    tataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTAAATTGACGCTAACTGACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGc                                                                                                <  2:3564633/99‑1 (MQ=255)
     ataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAACCATCCATACCTTCATTGACGCTAACTTACCGTTTGGTAGGATGAAGCAGTCAGGACCGGGCCg                                                                                              <  2:545664/100‑1 (MQ=255)
      taGCGCTCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCCTACCTTCATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGcgg                                                                                              >  2:3904742/1‑97 (MQ=255)
      taGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGt                                                                                             >  1:2644584/1‑100 (MQ=255)
         cGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATAGCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAt                                                                                          >  2:2260709/1‑100 (MQ=255)
            tCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGAtttt                                                                                       <  1:3321997/100‑1 (MQ=255)
            tCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAtttt                                                                                       <  2:2963922/100‑1 (MQ=255)
             cGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACATTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGGCAGGAACGGGCCGTGATTTTg                                                                                      <  2:1165427/100‑1 (MQ=255)
              gCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTgg                                                                                     <  1:2890763/100‑1 (MQ=255)
               cTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGc                                                                                    <  1:1317392/100‑1 (MQ=255)
               cTTACAGGCAGGGACGGTGTGGGCAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGc                                                                                    >  1:2327635/1‑100 (MQ=255)
                ttACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGcc                                                                                   <  1:2315785/100‑1 (MQ=255)
                ttACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGcc                                                                                   <  2:2986782/100‑1 (MQ=255)
                 tACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGccc                                                                                  >  1:357550/1‑100 (MQ=255)
                 tACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGGTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGccc                                                                                  >  2:3383339/1‑100 (MQ=255)
                   caggcaggGACGGTGTGGGTAAACAGCCATACCTTACTTGACGCTAACTTACCGGTTGGGGGGATGAAGCAGTCAGGACCGGGCCGTGATGTTGGCCCCg                                                                                >  1:3259125/1‑100 (MQ=255)
                   caggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGcc                                                                                   >  2:3438949/1‑97 (MQ=255)
                   caggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCg                                                                                >  2:150310/1‑100 (MQ=255)
                   caggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGAAGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGGTTTTGGCCCCg                                                                                >  2:2307356/1‑100 (MQ=255)
                     ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGAc                                                                              <  2:1746119/100‑1 (MQ=255)
                     ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGGGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGAc                                                                              >  2:3759078/1‑100 (MQ=255)
                        aggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActgg                                                                           >  2:289013/1‑100 (MQ=255)
                         gggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActggc                                                                          <  1:3389286/100‑1 (MQ=255)
                          ggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActggct                                                                         <  2:3324602/100‑1 (MQ=255)
                             cGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGa                                                                      <  2:1410383/100‑1 (MQ=255)
                             cGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGa                                                                      <  2:2108479/100‑1 (MQ=255)
                              ggTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActggctggtc                                                                     <  1:2791979/100‑3 (MQ=255)
                              ggTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGAc                                                                     <  1:1440859/100‑1 (MQ=255)
                              ggTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGAc                                                                     <  1:1576938/100‑1 (MQ=255)
                                 gtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt                                                                  >  2:1720729/1‑100 (MQ=255)
                                 gtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt                                                                  >  1:3989522/1‑100 (MQ=255)
                                  tgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTCGACGGtt                                                                 <  1:2161223/100‑1 (MQ=255)
                                   gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTg                                                                >  2:1877566/1‑100 (MQ=255)
                                   gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCTGTGATTTTAGCCCCGACTGGCTGGACGGTTg                                                                >  2:2781008/1‑100 (MQ=255)
                                   gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGtt                                                                 >  1:425257/1‑99 (MQ=255)
                                   gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTCTGGCGGGATGAAGCAGTCAGGAACGGGCCGTGATTATGGCCGCGACTGGCTGGACGGTTg                                                                >  2:864281/1‑100 (MQ=255)
                                    ggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg                                                               <  2:3621510/100‑1 (MQ=255)
                                     gTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgt                                                              >  2:1347836/1‑100 (MQ=255)
                                     gTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgt                                                              >  2:386003/1‑100 (MQ=255)
                                         aCAGCCATAGCTTAATTGACGCTAACTTGCCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGa                                                          <  2:3630667/100‑1 (MQ=255)
                                           aGCCATACCTTAATTGACGCGAACTTACCGTTTGTTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaa                                                         <  2:1916272/99‑1 (MQ=255)
                                            gCCATACGTTAATTGACGCTAACGTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaaa                                                        <  2:915110/99‑1 (MQ=255)
                                            gCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAAc                                                       <  1:2317904/100‑1 (MQ=255)
                                            gCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAAc                                                       >  2:3278468/1‑100 (MQ=255)
                                             ccATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACt                                                      >  2:1054535/1‑100 (MQ=255)
                                               aTACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTaa                                                    <  2:3109013/100‑1 (MQ=255)
                                                tACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAg                                                   <  1:268742/100‑1 (MQ=255)
                                                  ccTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                                 >  2:4111337/1‑100 (MQ=255)
                                                    ttAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCg                                                <  1:2157400/99‑1 (MQ=255)
                                                        ttGGCGCTAACTTACCGTGTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgt                                            <  2:19860/99‑1 (MQ=255)
                                                          gACGCTAACTTACCGTTTGGTGGGAGGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCGGGACGGTTGGTGTGAAACTAAGTCGgtgtgtg                                         <  2:580834/100‑1 (MQ=255)
                                                            cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTa                                       <  1:1168778/100‑1 (MQ=255)
                                                            cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTa                                       >  2:2975704/1‑100 (MQ=255)
                                                            cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTa                                       >  2:519649/1‑100 (MQ=255)
                                                            cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTa                                       >  1:1638317/1‑100 (MQ=255)
                                                               tAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACgg                                    <  1:3064034/100‑1 (MQ=255)
                                                                 aCTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTa                                  <  1:2548041/100‑1 (MQ=255)
                                                                  cTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAt                                 >  1:1054801/1‑100 (MQ=255)
                                                                   ttACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAtt                                >  1:1984902/1‑100 (MQ=255)
                                                                    tACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTa                               <  2:3596691/100‑1 (MQ=255)
                                                                     aCCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTa                               >  2:82842/1‑99 (MQ=255)
                                                                        gTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCt                           >  2:3100262/1‑100 (MQ=255)
                                                                          ttGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTgg                         <  2:3225781/100‑1 (MQ=255)
                                                                             gtggGGTGAAGGAGTCAGGACCGGGCCGTGATTTTGGGCGCGTCTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTc                      <  2:3177505/100‑1 (MQ=255)
                                                                             gtggGATGAAGCAGTGGGGACCGGGCCGTGATTTTGGGCCCGACTGGCTGGCCGGTGGGTGTGAAACTAAGTCGGGGTGTGTACGGTATTAATCTGGTTc                      <  2:3177364/100‑1 (MQ=255)
                                                                             gtggGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTc                      <  1:2216045/100‑1 (MQ=255)
                                                                              tggGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCg                     >  1:1151881/1‑100 (MQ=255)
                                                                              tggGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGATACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCg                     <  1:879665/100‑1 (MQ=255)
                                                                              tggGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGATGGTGTGAAACTAAGTCGGTGTGTGTACGGAATTAATCTGGTTCg                     >  2:592879/1‑100 (MQ=255)
                                                                                  aTGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCa                 <  1:3240863/100‑1 (MQ=255)
                                                                                  aTGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCa                 >  1:1169821/1‑100 (MQ=255)
                                                                                   tGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCAt                >  1:3068036/1‑100 (MQ=255)
                                                                                    gAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa               <  1:3126420/100‑1 (MQ=255)
                                                                                    gAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa               >  1:494989/1‑100 (MQ=255)
                                                                                       gCAGTCAGGACCGGGCCGTTGACCTAGACCCGACTGGCTGGACGGTTGGAGTGAAACTGAGACGGTGTGTGGACGGACTTACTCTGGTTCGCTCATAAGt            >  2:728949/1‑100 (MQ=255)
                                                                                       gCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGt            >  2:3302382/1‑100 (MQ=255)
                                                                                          gTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTTGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTa           <  2:1582751/98‑1 (MQ=255)
                                                                                            cAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaa       <  2:1951265/100‑1 (MQ=255)
                                                                                              ggACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAAc     >  2:590283/1‑100 (MQ=255)
                                                                                              ggACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAAc     >  1:3430555/1‑100 (MQ=255)
                                                                                              ggAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAAc     >  1:537007/1‑100 (MQ=255)
                                                                                                 aCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGc  <  1:3896594/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
GGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGC  >  NC_000913/1448850‑1449046

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