Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,632,588 A→T 6.4% intergenic (‑53/+16) guaA ← / ← guaB GMP synthetase/inosine 5'‑monophosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,632,5880AT 6.4% 171.6 / 3.7 64intergenic (‑53/+16)guaA/guaBGMP synthetase/inosine 5'‑monophosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (25/34);  new base T (2/2);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-02

ACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACC  >  NC_000913/2632452‑2632729
                                                                                                                                        |                                                                                                                                             
acgcacgcGGCGCGCAACCAGTTGAGTGTACTGCGAACCGAAGTCCAGAATGAGGATGTGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCCAGTGAAACAGCTAATCTAACTCGCCCGACTTCATGTCggg                                                                                                                                         <  2:680078/143‑1 (MQ=255)
acgcacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCggg                                                                                                                                         >  1:352783/1‑143 (MQ=255)
   cacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCGa                                                                                                                                      <  1:542290/143‑1 (MQ=255)
      gcggcgCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGCCGCGCGAaga                                                                                                                                   >  2:424593/1‑143 (MQ=255)
      gccgcgcgcAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAaga                                                                                                                                   <  1:850542/140‑1 (MQ=255)
       gggcgCGCTGCCAGTGGTGTGTACTTTGACCCGATGTCCAGTATGATGGTGCGATGCTGATGTAGGTTTGCCGTCACTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAagag                                                                                                                                  <  2:230568/142‑1 (MQ=255)
           gcgcAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGACAATATAAATCGCCCGACATGAACTCGGGCGAAGAGAATc                                                                                                                              >  2:672987/1‑143 (MQ=255)
            cgcAACCCGTTGAGTGTACTGAGAACCGAAGGCGAGAATGAGGATGCGATGCTTATGAAGGTTTTCCGTCATGGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCa                                                                                                                             <  1:307160/143‑1 (MQ=255)
             gcAACCAGTGGAGTGTACTGAGAACCGAAGTCCAGAATGAGCATGCGAGGCTTATGAATGTTTATCGTCATTGACGCTTATTTCGAGGCAAATGAAATAGATATTATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAg                                                                                                                            <  1:340399/143‑1 (MQ=255)
              cAACCAGTTGAGTGTACTGAGAACCGAAGTCCAAGATGAGTATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGCCGCAGAGAATCAgg                                                                                                                           >  1:919500/1‑143 (MQ=255)
               aagcAGGTGAGGGGACTGAGAACCGAAGTCCAGAAGGGGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAACATAAATCGCCGGTCACGAAGTCGGGCGAAGAGAATCAGGa                                                                                                                          <  1:702041/140‑1 (MQ=255)
                 ccAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGaac                                                                                                                        >  1:589760/1‑141 (MQ=255)
                    gTTGAGTGTACTGAGAACCGACGTCCAGACTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGTCATGAAGTCGGGCGAAGAGAATCAGGAGCCCa                                                                                                                     <  2:477993/143‑1 (MQ=255)
                      tGAGTGTACTGAGATCCGACGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGa                                                                                                                   <  2:425500/143‑1 (MQ=255)
                         gtgtACTGAGAACCGAAGTCCAGAATGAGGATGCTATGCTTATGAAGGTTTTCCGTCATTGACGCTTATCCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACgg                                                                                                                <  2:537254/143‑1 (MQ=255)
                               tGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAAAGAATCAGGAGCCCAGACGGTAGTTc                                                                                                          >  2:718528/1‑143 (MQ=255)
                                         aGTCCCGAATGAGGATGCGATGCTTATGAATGTTTTCCGGCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCtt                                                                                                <  2:347545/143‑1 (MQ=255)
                                           tCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTa                                                                                              <  1:556500/143‑1 (MQ=255)
                                           tCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTa                                                                                              >  1:793418/1‑143 (MQ=255)
                                           tCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGAGAATAGAATCAGGAGCCCAGACGGTAGTCCGGGGACTCTTTa                                                                                              >  2:790868/1‑143 (MQ=255)
                                             cAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGGAGAGAATCAGGAGCCCAGTCGGGAGTTCGGGGACTCTTTAgg                                                                                            >  1:860680/1‑142 (MQ=255)
                                                 aTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATg                                                                                        <  1:538304/143‑1 (MQ=255)
                                                  tGAGGAGTCGAGGCTTATGGATGTTTTCCGTCAGTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATGAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATgg                                                                                       <  2:182412/143‑1 (MQ=255)
                                                   gAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGt                                                                                      <  1:580142/143‑1 (MQ=255)
                                                            tggtgtATGAATGTTCTCCGTCATGGGCGCTTATTCCGAGGCAAGTGAAACAGATAATATAAAGCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTg                                                                             <  2:266429/138‑1 (MQ=255)
                                                            atgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTg                                                                             <  1:849051/143‑1 (MQ=255)
                                                              gcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTACCATGCACGTCGTGaa                                                                           >  2:789939/1‑143 (MQ=255)
                                                               cTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAAc                                                                          <  2:768703/143‑1 (MQ=255)
                                                                ttATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTCTAGTAATGGTCACGTCGTGAACg                                                                         >  2:490860/1‑143 (MQ=255)
                                                                   tGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGCCTCTTTAGTAATGGTCTCGCCGTCCACGTgg                                                                      >  1:90524/1‑143 (MQ=255)
                                                                   tGAACGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTgg                                                                      <  1:559472/143‑1 (MQ=255)
                                                                       tGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGCCATAAAGTCGGGCGCAGAGAATCAGGAGCCCAGACGGTAGTTCCGGGACTCTTTAGTAATGGTCACGTCGCGACCGTGGCttt                                                                  >  2:852859/1‑143 (MQ=255)
                                                                        gTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTc                                                                 <  1:753224/143‑1 (MQ=255)
                                                                          tttGCGTCATTGACGCTTATTCCGATGTAGTTGAAACAGATAATATAAATCGCCGGACATGTAGTCGTGCGAAGAGAATCAGGAGCCCAGACGGTATTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCt                                                               <  2:848827/143‑1 (MQ=255)
                                                                          gtcccGCCACTGACGCTCAGTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCt                                                               <  1:232490/140‑1 (MQ=255)
                                                                               gTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATg                                                          <  2:352783/143‑1 (MQ=255)
                                                                                   ttGACGCGTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg                                                      <  2:685683/143‑1 (MQ=255)
                                                                                   ttGACGCGTATCCCGAGGCAAGCGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg                                                      <  1:317652/143‑1 (MQ=255)
                                                                                         cTTATCCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGc                                                <  1:317414/143‑1 (MQ=255)
                                                                                          ttATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGATCCTAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGGCCGCACCGCt                                               >  2:692661/1‑143 (MQ=255)
                                                                                          gtagtCCGAGGCAAGTGAAGCAGAGAATATAAGTCGCGGGACATGAGGTCGGGCGGAGAGAAGCAGGAGCCCAGACGGTAGTTCGGGGCCTCTTTAGTACAAGTATCGGTCTCACCGCGTCTTTCCTGAATGCCCGCACCGCt                                               <  2:506215/139‑1 (MQ=255)
                                                                                                cGAGGCAAGTGAAACAGATAATATAAATCGCACGACATGAAGGAGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGTACTCTTTAGTAATGTTCACGTCGTGAAAATGGCTTTCCTGAATGCCAGAACCGCTGATaag                                         >  2:61755/1‑141 (MQ=255)
                                                                                                  aGGCAAGTGACACAGATAATATAAACCGCCCGCCATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTa                                       <  1:269737/143‑1 (MQ=255)
                                                                                                      aaGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGGGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACaaa                                   >  1:641481/1‑143 (MQ=255)
                                                                                                         tGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGACCGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACta                                >  2:309760/1‑142 (MQ=255)
                                                                                                           aggcAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCg                              <  1:66445/140‑1 (MQ=255)
                                                                                                           aaaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCg                              <  1:639033/143‑1 (MQ=255)
                                                                                                             aCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCGCCGCTGATACGTACAAACTCCGCt                            >  1:337199/1‑143 (MQ=255)
                                                                                                              cAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCtt                           <  1:218154/143‑1 (MQ=255)
                                                                                                               aGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGCATCAGGAGCCCAGACGGTAGGTCGTGGGCGCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCttt                          >  1:577858/1‑143 (MQ=255)
                                                                                                                gctaataTGAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTa                         <  1:810992/141‑1 (MQ=255)
                                                                                                                gATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTa                         <  1:220978/143‑1 (MQ=255)
                                                                                                                 ataataTAAATCGGCCGGCGTGAAGTCGGGCGAAGAGAATCACGAGCCCAGACGGTAGTGCGGGGACTCTTTCGTCATGGTCACGTCGTGAACGCGGCTTTCCTGAATGCCCGCACCGCTGACACGTACAAACTCCGCTTTAg                        <  2:769040/143‑1 (MQ=255)
                                                                                                                  tgacataAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGt                       <  1:481786/139‑1 (MQ=255)
                                                                                                                  taataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGACCGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACATACTCCGCttttct                       >  1:478695/1‑140 (MQ=255)
                                                                                                                  taataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGCATGCCCGCACCGCTGATACGTACAAACTCCGCtataat                       >  1:788734/1‑141 (MQ=255)
                                                                                                                  taataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGCAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACGATCTCCGCTTTAGt                       >  2:220170/1‑143 (MQ=255)
                                                                                                                           tCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTtc              >  1:62361/1‑143 (MQ=255)
                                                                                                                             gCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGAACGCGGTtcgt            >  2:213410/1‑143 (MQ=255)
                                                                                                                             gCCCGACTTCATGTCCGGCGAAGAGAATCAGGAGCCCACACGGTAGTTCCGCCACTCTGTAGTAACGGTCACCTCCTCAACGCCCCTCTCCTCAATGCCCGCACCGCTGATACGTACAACCCCCGCTCTAGTCCGCCGtccct            >  2:876689/1‑139 (MQ=255)
                                                                                                                                cGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCCCGTCATGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGCCGa         >  1:575542/1‑143 (MQ=255)
                                                                                                                                     ggAAGTCCGGCGGAGAGAATCAGGTGCCCAGACGGTAGGTCGGGGGCTCTTGAGTAAGGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTa    <  2:605205/142‑1 (MQ=255)
                                                                                                                                      gAAGTCGGTCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTAc   <  1:807038/143‑1 (MQ=255)
                                                                                                                                      gAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCCGTTCGTCGATAGTAc   >  1:349504/1‑143 (MQ=255)
                                                                                                                                       aaGTCGGGCGAAGAGAATCAGGAGCCCAGACGTTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTAcc  <  2:40653/143‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                             
ACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACC  >  NC_000913/2632452‑2632729

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: