Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,004,515 G→C 100% G61G (GGC→GGG yigI ← 4HBT thioesterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,5150GC51.5% 23.6 / 278.8 174G61G (GGC→GGGyigI4HBT thioesterase family protein
Reads supporting (aligned to +/- strand):  ref base G (41/43);  new base C (41/48);  total (83/91)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 177 total aligned reads displayed.

GTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCC  >  NC_000913/4004421‑4004614
                                                                                              |                                                                                                   
gTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTg                                                                                                <  2:587522/100‑1 (MQ=255)
      cATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCTATTACCCCCCCGtcc                                                                                               >  1:3940026/1‑93 (MQ=255)
      cATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAAt                                                                                          <  2:1252999/100‑1 (MQ=255)
      cATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAAt                                                                                          >  1:3455162/1‑100 (MQ=255)
       aTCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATg                                                                                         <  2:3879401/100‑1 (MQ=255)
         cTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCt                                                                                       >  2:2590921/1‑100 (MQ=255)
          ttCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCtt                                                                                      >  2:3598578/1‑100 (MQ=255)
           tCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGCCATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCttt                                                                                     >  2:2571716/1‑100 (MQ=255)
            cACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTg                                                                                    <  2:3522433/100‑1 (MQ=255)
            cACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTg                                                                                    >  1:1197427/1‑100 (MQ=255)
             aCTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTgc                                                                                   >  1:3457613/1‑100 (MQ=255)
             aCTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTgc                                                                                   >  2:3353305/1‑100 (MQ=255)
              cTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTgcg                                                                                  <  2:420585/100‑1 (MQ=255)
              cTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTgcg                                                                                  >  1:2680727/1‑100 (MQ=255)
              cTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTgcg                                                                                  <  1:711955/100‑1 (MQ=255)
               tGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTgcgc                                                                                 >  1:2144920/1‑100 (MQ=255)
                gATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGAACGGCGGCGACATCCCGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGcc                                                                                >  1:31265/1‑100 (MQ=255)
                 aTGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGccc                                                                               <  2:2749704/100‑1 (MQ=255)
                 aTGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGccc                                                                               <  1:450046/100‑1 (MQ=255)
                 aTGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGccc                                                                               <  2:241354/100‑1 (MQ=255)
                   ggTTTCGTGGCGGGTTAAGGTACTTCCCCCGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAg                                                                             <  2:191906/100‑1 (MQ=255)
                    gTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGt                                                                            >  2:230032/1‑100 (MQ=255)
                    gTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGt                                                                            <  2:538980/100‑1 (MQ=255)
                    gTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGt                                                                            <  1:2826627/100‑1 (MQ=255)
                    gTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGt                                                                            >  2:2781906/1‑100 (MQ=255)
                    gTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGt                                                                            <  2:3715875/100‑1 (MQ=255)
                    gTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGt                                                                            <  2:4121728/100‑1 (MQ=255)
                    gTTTCGTGGCCGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAg                                                                             <  2:3167546/99‑1 (MQ=255)
                      ttCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTg                                                                          <  1:2159432/100‑1 (MQ=255)
                      ttCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTg                                                                          <  1:2922295/100‑1 (MQ=255)
                        cGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGcc                                                                        >  2:1959490/1‑100 (MQ=255)
                         gTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGccc                                                                       <  2:198016/100‑1 (MQ=255)
                         gTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGccc                                                                       <  1:1030918/100‑1 (MQ=255)
                         gTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGccc                                                                       >  1:3224520/1‑100 (MQ=255)
                          tggctcgTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCcca                                                                      <  1:1855847/94‑1 (MQ=255)
                           ggCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccac                                                                     <  2:703532/100‑1 (MQ=255)
                                gTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACcatca                                                                >  2:2104224/1‑100 (MQ=255)
                                gTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACcatca                                                                >  1:33593/1‑100 (MQ=255)
                                 ttAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACcatcat                                                               >  2:1346605/1‑100 (MQ=255)
                                   aaGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTg                                                             <  2:2587432/100‑1 (MQ=255)
                                   aaGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTg                                                             <  2:2795393/100‑1 (MQ=255)
                                   aaGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTg                                                             >  1:577024/1‑100 (MQ=255)
                                    aGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTgg                                                            >  1:1079737/1‑100 (MQ=255)
                                     ggTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGc                                                           >  2:886656/1‑100 (MQ=255)
                                     ggTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGc                                                           <  1:2206307/100‑1 (MQ=255)
                                        aCTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGa                                                        >  2:4072186/1‑100 (MQ=255)
                                         cTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAt                                                       >  2:155601/1‑100 (MQ=255)
                                         cTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAt                                                       <  1:447427/100‑1 (MQ=255)
                                         attCTCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGAGCCCGCCGTGCAAAAGGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAt                                                       >  1:787972/2‑100 (MQ=255)
                                          ttCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAtt                                                      >  1:2335053/1‑100 (MQ=255)
                                           tCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttt                                                     <  1:3155114/100‑1 (MQ=255)
                                            ccgtggCACACCAGCCCGGCGGCGACAGGCAGGGCCGACGCAATGAACCCCCCGTGGAAAATGCGTTGCGCCCAGTTGCCCACCATCATTGGc                                                           <  1:4094780/88‑1 (MQ=255)
                                              cccGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTa                                                  <  1:117327/98‑1 (MQ=255)
                                              cACGCACACCAGACCGGCGGCGACATCTAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTa                                                  <  1:1977407/100‑1 (MQ=255)
                                                 gCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAg                                               >  2:3981213/1‑100 (MQ=255)
                                                   acacCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGc                                             >  1:600446/1‑100 (MQ=255)
                                                    cacCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGcc                                            <  1:123839/100‑1 (MQ=255)
                                                      ccAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAg                                          >  2:844597/1‑100 (MQ=255)
                                                      ccAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAg                                          >  2:338328/1‑100 (MQ=255)
                                                       cAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGc                                         <  1:595670/100‑1 (MQ=255)
                                                       cAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGc                                         <  1:285117/100‑1 (MQ=255)
                                                         gACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtg                                       >  1:3603918/1‑100 (MQ=255)
                                                          aCCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgt                                      >  1:2424660/1‑100 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                     >  1:1830403/1‑100 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                     <  2:4258159/100‑1 (MQ=255)
                                                            cggcggCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGc                                    <  2:199759/100‑1 (MQ=255)
                                                             ggcggcGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCg                                   >  1:3908351/1‑100 (MQ=255)
                                                             ggcggcGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCg                                   >  2:2363537/1‑100 (MQ=255)
                                                             ggcggcGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCg                                   >  2:2637379/1‑100 (MQ=255)
                                                                ggcgACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTTCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGaa                                 >  1:3007584/1‑99 (MQ=255)
                                                                 gcgACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAAct                               <  1:173624/100‑1 (MQ=255)
                                                                       tCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCttttt                         >  2:1495417/1‑100 (MQ=255)
                                                                       tCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCttttt                         <  2:3455175/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                        >  1:399899/1‑100 (MQ=255)
                                                                        ccAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                        <  2:2907658/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                        >  2:2941723/1‑100 (MQ=255)
                                                                         cAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCg                       <  2:600448/100‑1 (MQ=255)
                                                                         cAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCg                       <  2:1497792/100‑1 (MQ=255)
                                                                          aGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGt                      <  1:2740645/100‑1 (MQ=255)
                                                                          aGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGt                      <  2:3914897/100‑1 (MQ=255)
                                                                            cgcCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTaa                    >  2:41735/1‑100 (MQ=255)
                                                                            cgcCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGATCTCTTTTTCGtta                    <  2:2898231/100‑3 (MQ=255)
                                                                             gcCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAAc                   <  1:1255740/100‑1 (MQ=255)
                                                                              ccGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAAc                   >  1:981034/1‑99 (MQ=255)
                                                                              ccGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACg                  <  2:1198385/100‑1 (MQ=255)
                                                                              ccGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATCGGCTGATTTTTAATGTCTTTTTGTGCGAACTCTTTTTc                        >  1:241760/1‑94 (MQ=255)
                                                                               cGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTaa                    <  1:4222037/97‑1 (MQ=255)
                                                                               cGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGc                 >  2:3424431/1‑100 (MQ=255)
                                                                                  cGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTcc              <  1:1645652/100‑1 (MQ=255)
                                                                                   gCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTctttttttt                >  1:550327/1‑89 (MQ=255)
                                                                                    cAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAg            <  1:978303/100‑1 (MQ=255)
                                                                                    cAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAg            >  1:1404970/1‑100 (MQ=255)
                                                                                    cAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAg            <  1:2071757/100‑1 (MQ=255)
                                                                                       tGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTc         <  2:3975113/100‑1 (MQ=255)
                                                                                        gACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTcc        >  2:3051490/1‑100 (MQ=255)
                                                                                        gACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTcc        >  1:2863284/1‑100 (MQ=255)
                                                                                         aCCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCa       <  2:2088082/100‑1 (MQ=255)
                                                                                          ccccccCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCAt      >  2:481009/1‑100 (MQ=255)
                                                                                             cgccgTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATccc   <  1:3557198/100‑1 (MQ=255)
                                                                                              gccgTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATcccc  <  2:3749294/100‑1 (MQ=255)
                                                                                              |                                                                                                   
GTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCC  >  NC_000913/4004421‑4004614

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: