Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,051 C→G 6.7% intergenic (+34/‑164) lldD → / → trmL L‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,0510CG 6.7% 164.7 / 3.2 60intergenic (+34/‑164)lldD/trmLL‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/28);  new base G (3/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

ACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCC  >  NC_000913/3780910‑3781191
                                                                                                                                             |                                                                                                                                             
aCTGGCGCGAAATCCATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                                              <  1:216699/143‑1 (MQ=255)
  tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGGGGGGCCAACATCCCCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                            <  1:491814/143‑1 (MQ=255)
  tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                            <  2:580950/143‑1 (MQ=255)
       cgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGTACTGGCTCCCCTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCATCCTAACCCTCTCCCTGAGGGagag                                                                                                                                       <  1:763343/143‑1 (MQ=255)
        gAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGagagc                                                                                                                                      >  2:534060/1‑142 (MQ=255)
           aTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGAAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                                                                   <  2:528489/143‑1 (MQ=255)
           aTCGATCAGCGAAATTACGCAAGAGTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                                                                   <  1:62463/143‑1 (MQ=255)
            tCGATCAGCGAACGTACGCAAGATTCGGTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGAc                                                                                                                                  <  2:101525/143‑1 (MQ=255)
              gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCCTGGCGAAAGGGAATGCGGCATAGTCGTTTGCTCCCCCTCACCCTAACGCTCTCCCTGGGGGAGGGGGGACCg                                                                                                                                >  1:397237/1‑143 (MQ=255)
                   gCGAACTTACGCAACATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCgg                                                                                                                           <  1:323319/143‑1 (MQ=255)
                   gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCgg                                                                                                                           <  1:1014310/143‑1 (MQ=255)
                    cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                                                                          >  2:968659/1‑143 (MQ=255)
                     gAAATTACGCAAGCTTCGCTGGTGCAGGGGCTGGGTACAGACTGGCCTGCGCCACCGGCTCCCATGGCGAGAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                                                                        <  2:336659/143‑1 (MQ=255)
                        aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                                                                     <  1:920819/143‑1 (MQ=255)
                           aCGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGATGCCTGCGGCACTGGCTCCCATGGCCAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtat                                                                                                                  <  1:541066/143‑1 (MQ=255)
                             gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCATCAGGGAGAGGGGACCGTTCGCCGCcgtatgt                                                                                                                >  2:598078/1‑143 (MQ=255)
                              cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                                               <  1:967433/143‑1 (MQ=255)
                                aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCGTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGGCGTTTGCGGCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACt                                                                                                             <  2:1046351/143‑1 (MQ=255)
                                   tgcGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGTGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                          <  2:152643/141‑1 (MQ=255)
                                    tCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCACTCAGGGAGTGGGTTTCTTTCGGCGCTGTATGTACTCCCt                                                                                                         <  1:1022541/143‑1 (MQ=255)
                                              aGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCAAACCCTGaaa                                                                                               >  2:458548/1‑143 (MQ=255)
                                               ggggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCGCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAAc                                                                                              <  2:30697/143‑1 (MQ=255)
                                                    tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACACACTCACTCTGAAACGACAcc                                                                                        >  1:520364/1‑143 (MQ=255)
                                                     gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTCTGTACTCCCTCACTCTGAAACGCCAcc                                                                                        >  2:23753/1‑143 (MQ=255)
                                                        tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCa                                                                                     >  1:889757/1‑143 (MQ=255)
                                                        tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCCGCGCGGTATGAACTCCCTCACTCTGAAGCGACACCGcc                                                                                     >  2:472446/1‑142 (MQ=255)
                                                          aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGGCCGTTCGGCGCTGTCTGGACTCCCTGCCTCTGAACCGAGCCCGCAct                                                                                   >  1:669084/1‑143 (MQ=255)
                                                           agcggtGCCTGCGGCACGGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCGCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCActc                                                                                  <  2:293329/138‑1 (MQ=255)
                                                           agagTGGCCGGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCActc                                                                                  <  1:967784/143‑1 (MQ=255)
                                                              gTTGCGTGCGGCGCTGGCGCCCAGGGGGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCTCACCCTACCCCTCTCCCCCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                               <  2:491772/143‑1 (MQ=255)
                                                                  ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                           <  1:945314/143‑1 (MQ=255)
                                                                   cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                                                                          >  1:886332/1‑143 (MQ=255)
                                                                    tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACCCTGAAACGACACCGCAATCTTTTTTtct                                                                         >  1:717237/1‑143 (MQ=255)
                                                                     gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACACCGCACTCTTTTtctctc                                                                        >  2:992176/1‑143 (MQ=255)
                                                                       ggCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGCGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAGCGACACCGCACTCTTTTTTTCTccc                                                                      >  2:497306/1‑143 (MQ=255)
                                                                           cTGGCTCCCATGGCCAAAGCGAATGCTGCATATTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                  <  1:633619/143‑1 (MQ=255)
                                                                                tCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGCGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTc                                                             <  1:534060/143‑1 (MQ=255)
                                                                                tCCCAAGTCGGAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACACTAACCCTCTCCCTCAGGGAGAGGGGGCCGTTCGGCGCGGTATGTACTCCCTCCATCTGAAACGACACCGCACCCTTTTTTTCTCCCTCGcccccc                                                             >  2:435067/1‑141 (MQ=255)
                                                                                  ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCg                                                           <  1:1049465/143‑1 (MQ=255)
                                                                                  ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGCAGAGGCGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGCAACGACCCCGCACTCTTTTTTTCTCCCTCGCACCTcag                                                           >  2:669053/1‑141 (MQ=255)
                                                                                   cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGCAGAGGCGACCGTTCGGCGCTGTATGTACTCCCTCACTCTCAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                          >  2:153943/1‑143 (MQ=255)
                                                                                    gtcgCGAAAGGGGATGCGGCACGGGCGTTTGCGCCGCCGCACCCTACCCCCCGCCCGGAGGGAGAGGGGCCCGTTCGGCGCCGCATGTACTCCCTCACCCTGAAACGACACCGCACTCTTTTTTTCTCCCGCGCCCCTCCggg                                                         <  2:860143/140‑1 (MQ=255)
                                                                                    aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                         <  2:551316/143‑1 (MQ=255)
                                                                                        cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGTGgaga                                                     >  1:909058/1‑143 (MQ=255)
                                                                                         gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACAGCCCTCTTTTTTTCTCCCTCGCCCCTCCTGGGagag                                                    >  2:203909/1‑143 (MQ=255)
                                                                                            aGGGACTGCGGCAAAGGCGGTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGGACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGc                                                 <  1:948992/143‑1 (MQ=255)
                                                                                              ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCGGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCCCCCGGGACAGGCCCg                                               >  2:226215/1‑143 (MQ=39)
                                                                                              ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTATCTGAGATTGAACCCTTCGACAAACGTGTTGTGTCCGCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                               >  1:1070772/1‑143 (MQ=255)
                                                                                                 aTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggg                                            <  1:919827/143‑1 (MQ=255)
                                                                                                    cGGCATAGTCGTTTGCGCCCCCTCACCCTGACGCGCTCCCTCAGGGAGGGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCGCTCTTTTTTTCTCCCCCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                         <  2:118865/143‑1 (MQ=255)
                                                                                                      gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGCCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGAGCCCGCCCTCTTTTTTGCCCCCTCGCCTCTCCGGGGACGGGGCCGGGGTGAgg                                       >  2:584422/1‑143 (MQ=255)
                                                                                                         tAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGGCGGGGTGAggggg                                    >  1:234288/1‑142 (MQ=255)
                                                                                                            tCGTTTGCGCCCCCTCACCCTACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                 <  1:858755/143‑1 (MQ=255)
                                                                                                            tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCGGGTTCGTGAACGTCAGAGTTAGCTCTTGTGCTCCCTCGCCCCTCCGGGGGGAGGGCCGGGGGGAGGGGaaga                                 >  2:1044374/1‑141 (MQ=19)
                                                                                                                       cccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGc                      >  2:1031454/1‑143 (MQ=255)
                                                                                                                                   cccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGCGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc          >  2:442481/1‑143 (MQ=255)
                                                                                                                                    ccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt         >  2:336122/1‑143 (MQ=255)
                                                                                                                                    ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGAAGAGGCCCGGGGTGAGGGGAACAGACCGAAATGCTGACTAATTTTCt         >  1:661865/1‑143 (MQ=255)
                                                                                                                                     ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGCGAGAGGGCCGGCGTGAGGGGAGAAGGCCGCACTGCTCACCAATCTtgtg        >  2:882280/1‑143 (MQ=255)
                                                                                                                                           cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAATCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGCGGGGGCCGGGGTGAGGGGAAAAGGCCGCCCTGCTCACCAATTTTCTGCTATcc  >  1:228289/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                             
ACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCC  >  NC_000913/3780910‑3781191

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