Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,045,986 Δ1 bp intergenic (‑85/+182) ygfF ← / ← gcvP putative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/glycine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,045,9830T.98.9% 408.5 / ‑1.8 94intergenic (‑82/+185)ygfF/gcvPputative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/glycine decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (34/59);  total (34/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

GTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCAT  >  NC_000913/3045852‑3046117
                                                                                                                                   |                                                                                                                                       
gTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAATAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCcccct                                                                                                                                 <  1:976967/140‑2 (MQ=255)
 ttGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTc                                                                                                                               >  1:1156073/1‑140 (MQ=255)
 ttGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTc                                                                                                                               >  2:681002/1‑140 (MQ=255)
  tGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCg                                                                                                                              <  2:1567033/140‑1 (MQ=255)
  tGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCg                                                                                                                              <  1:899520/140‑1 (MQ=255)
  tGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCg                                                                                                                              <  1:1136815/140‑1 (MQ=255)
  tGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCg                                                                                                                              <  2:812747/140‑1 (MQ=255)
  tGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCg                                                                                                                              <  2:114340/140‑1 (MQ=255)
  tGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCg                                                                                                                              <  2:454048/140‑1 (MQ=255)
   gcccgcccGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGc                                                                                                                             >  2:897341/1‑140 (MQ=255)
   gcccgcccGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGc                                                                                                                             >  2:1538563/1‑140 (MQ=255)
   gcccgcccGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGc                                                                                                                             >  1:353134/1‑140 (MQ=255)
   gcccgcccGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGc                                                                                                                             >  2:148871/1‑140 (MQ=255)
   gcccgcccGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGc                                                                                                                             >  2:1443039/1‑140 (MQ=255)
   gcccgcccGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGc                                                                                                                             >  1:747062/1‑140 (MQ=255)
   gcccgcccGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGc                                                                                                                             >  1:1199782/1‑140 (MQ=255)
   gcccgcccGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGCTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGc                                                                                                                             >  2:28792/1‑140 (MQ=255)
    cccgcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGcc                                                                                                                            <  2:1391467/140‑1 (MQ=255)
    cccgcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGcc                                                                                                                            <  1:2403/140‑1 (MQ=255)
    cccgcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGcc                                                                                                                            <  1:270564/140‑1 (MQ=255)
    cccgcccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGACTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGcc                                                                                                                            <  2:402783/140‑1 (MQ=255)
      cgcccgCTGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCt                                                                                                                          <  2:1146082/140‑1 (MQ=255)
         ccgATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTc                                                                                                                       <  1:428899/140‑1 (MQ=255)
              gCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCTCGCTTTCCCCCTCGCCCTTTCAGGga                                                                                                                  <  1:828958/140‑1 (MQ=255)
              gCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGga                                                                                                                  <  1:518522/140‑1 (MQ=255)
                    aCTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGc                                                                                                            <  1:88145/140‑1 (MQ=255)
                     cTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGcc                                                                                                           <  1:1187000/140‑1 (MQ=255)
                       gCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCgg                                                                                                         >  2:1330983/1‑140 (MQ=255)
                       gCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCgg                                                                                                         >  1:770780/1‑140 (MQ=255)
                        ccaccaGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCggg                                                                                                        <  1:795399/140‑1 (MQ=255)
                          cccaGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGt                                                                                                      <  2:54515/139‑1 (MQ=255)
                               tCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAggg                                                                                                 <  1:85789/140‑1 (MQ=255)
                                 acaAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTa                                                                                               >  2:258343/1‑140 (MQ=255)
                                   aaGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGACACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTaaa                                                                                             >  1:817592/1‑140 (MQ=255)
                                       gCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATAtt                                                                                         <  1:1438158/140‑1 (MQ=255)
                                        cTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTc                                                                                        >  1:1531564/1‑140 (MQ=255)
                                        cTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTc                                                                                        >  2:499757/1‑140 (MQ=255)
                                         tataGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTcg                                                                                       >  2:163079/1‑140 (MQ=255)
                                           taGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTcgcg                                                                                     >  1:1408685/1‑140 (MQ=255)
                                            aGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCgcgc                                                                                    <  1:1354399/140‑1 (MQ=255)
                                             gCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGcc                                                                                   >  2:199200/1‑140 (MQ=255)
                                                aTTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGt                                                                                >  1:759559/1‑140 (MQ=255)
                                                aTTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGt                                                                                >  1:1197229/1‑140 (MQ=255)
                                                aTTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGt                                                                                >  2:1448694/1‑140 (MQ=255)
                                                   tGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGcc                                                                             >  1:1131874/1‑140 (MQ=255)
                                                    gTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCg                                                                            <  1:258343/140‑1 (MQ=255)
                                                      ttttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGc                                                                          <  1:1069127/140‑1 (MQ=255)
                                                         tCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTg                                                                       <  1:597786/140‑1 (MQ=255)
                                                          ccGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGt                                                                      <  1:1538563/140‑1 (MQ=255)
                                                            gagaAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTc                                                                    >  2:1352827/1‑140 (MQ=255)
                                                               aaaaaTTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTcccc                                                                 <  1:681002/140‑1 (MQ=255)
                                                                aaaaTTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCt                                                                <  1:163079/140‑1 (MQ=255)
                                                                aaaaTTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCt                                                                <  1:28792/140‑1 (MQ=255)
                                                                aaaaTTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCt                                                                <  1:55908/140‑1 (MQ=255)
                                                                aaaaTTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCt                                                                <  1:698181/140‑1 (MQ=255)
                                                                    ttaattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                            <  1:695241/140‑1 (MQ=255)
                                                                          cGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAcc                                                     <  1:870093/140‑1 (MQ=255)
                                                                          cGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAAGGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAcc                                                     <  2:466040/140‑1 (MQ=255)
                                                                            ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCGCGCTCTCCCTCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCTCCCTAAcccg                                                   <  2:1250144/140‑2 (MQ=255)
                                                                            ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCct                                                   <  2:775073/140‑1 (MQ=255)
                                                                            ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCct                                                   <  1:1443039/140‑1 (MQ=255)
                                                                            ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCct                                                   <  2:17420/140‑1 (MQ=255)
                                                                            ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCct                                                   <  1:1352827/140‑1 (MQ=255)
                                                                                   aTCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCa                                            <  2:1156073/140‑1 (MQ=255)
                                                                                     cTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                          <  2:219835/140‑1 (MQ=255)
                                                                                     cTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                          <  1:873679/140‑1 (MQ=255)
                                                                                       ttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAgg                                        <  1:199200/140‑1 (MQ=255)
                                                                                       ttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAgg                                        <  1:148871/140‑1 (MQ=255)
                                                                                         gagaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAgggg                                      >  2:460558/1‑140 (MQ=255)
                                                                                           gaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCg                                    <  2:1131874/140‑1 (MQ=255)
                                                                                           gaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCg                                    <  2:1408685/140‑1 (MQ=255)
                                                                                            aTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGa                                   <  1:885883/140‑1 (MQ=255)
                                                                                             ttGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                  <  2:1328498/140‑1 (MQ=255)
                                                                                              tGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                 >  1:581017/1‑140 (MQ=255)
                                                                                               gCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                >  1:1391780/1‑140 (MQ=255)
                                                                                               gCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                >  2:1140956/1‑140 (MQ=255)
                                                                                                 tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                              <  1:1448694/140‑1 (MQ=255)
                                                                                                 tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                              <  2:1181272/140‑1 (MQ=255)
                                                                                                  aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGAc                             <  2:666520/140‑1 (MQ=255)
                                                                                                          ggCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                     <  2:1199782/140‑1 (MQ=255)
                                                                                                          ggCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                     <  1:991630/140‑1 (MQ=255)
                                                                                                          ggCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                     <  1:971143/140‑1 (MQ=255)
                                                                                                          ggCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                     >  1:965924/1‑140 (MQ=255)
                                                                                                          ggCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                     <  1:919493/140‑1 (MQ=255)
                                                                                                          ggCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                     <  2:747062/140‑1 (MQ=255)
                                                                                                          ggCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                     >  2:567461/1‑140 (MQ=255)
                                                                                                          ggCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                     <  1:255306/140‑1 (MQ=255)
                                                                                                          ggCTCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATACTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                     <  1:328108/140‑1 (MQ=255)
                                                                                                             tCAATCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCAc                  >  2:496441/1‑140 (MQ=255)
                                                                                                             ccAATCTTCAATGGTGGCCCCGTTTCCCCCCCCCCCTTTCCGGGGGGGGGGCGGGGGGGGGGTAAAATTTCGCCCCAGTGCCCGCCTGGTCCCCTCCCCCTTACCCCCCCCCCCAAGGGGGCGAGGGGACTGTCTGAGCa                   <  2:353134/139‑1 (MQ=255)
                                                                                                               aaTCTTCAATGTTGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAATGGGCGAGGGGACTGTCTGAGCACtt                >  1:1039920/1‑140 (MQ=255)
                                                                                                                           tGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCTCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTc    >  2:1239630/1‑140 (MQ=255)
                                                                                                                           tGGCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTc    <  1:823580/140‑1 (MQ=255)
                                                                                                                            ggCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCa   >  2:599071/1‑140 (MQ=255)
                                                                                                                            ggCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCa   <  1:386657/140‑1 (MQ=255)
                                                                                                                             gCACGCTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCAt  <  1:1416415/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
GTTGCCCGCCCGATGCCGCGACTGCCACCAGTCACAAGTGCTATAGCCATTTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCAT  >  NC_000913/3045852‑3046117

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