Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,632,586 G→C 6.3% intergenic (‑51/+18) 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,5860GC 6.3% 180.1 / 4.9 65intergenic (‑51/+18)guaA/guaBGMP synthetase/inosine 5'‑monophosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base C (2/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

CCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTAC  >  NC_000913/2632444‑2632728
                                                                                                                                              |                                                                                                                                              
cccAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAAGGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATg                                                                                                                                                <  2:672026/143‑1 (MQ=255)
 ccAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAGCAGATAATATAAATCGCCCGACATGa                                                                                                                                               >  1:789284/1‑143 (MQ=255)
 ccAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGa                                                                                                                                               >  1:508481/1‑143 (MQ=255)
 ccAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGa                                                                                                                                               >  1:40653/1‑143 (MQ=255)
  cAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTCTTCCGAGGCAAGTGCACCCGCTAATATAAATCGCCCGCCCTGac                                                                                                                                              >  1:203271/1‑142 (MQ=255)
        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:793418/1‑143 (MQ=255)
                                                   tCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTa                                                                                             <  1:556500/143‑1 (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)
                                                                                                                                              |                                                                                                                                              
CCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTAC  >  NC_000913/2632444‑2632728

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: