Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,460,484 T→A 6.9% intergenic (‑120/+38) nrdG ← / ← nrdD anaerobic ribonucleoside‑triphosphate reductase activating protein/anaerobic ribonucleoside‑triphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,460,4840TA 6.9% 601.8 / 21.8 188intergenic (‑120/+38)nrdG/nrdDanaerobic ribonucleoside‑triphosphate reductase activating protein/anaerobic ribonucleoside‑triphosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (96/79);  new base A (4/9);  total (100/88)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 192 total aligned reads displayed.

GTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACACGCGGC  >  NC_000913/4460335‑4460631
                                                                                                                                                     |                                                                                                                                                   
 tCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTa                                                                                                                                                    <  1:875015/150‑1 (MQ=255)
 tCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTa                                                                                                                                                    <  2:31289/150‑1 (MQ=255)
     cAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                <  2:888204/150‑1 (MQ=255)
       agATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATTCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                              >  1:1218761/2‑150 (MQ=255)
           aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                          <  1:515609/150‑1 (MQ=38)
           aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                          >  2:790356/1‑150 (MQ=38)
           aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                          <  1:520548/150‑1 (MQ=38)
           aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                          >  1:148323/1‑150 (MQ=255)
            gTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                           <  1:334970/148‑6 (MQ=38)
            gTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                           >  2:334970/1‑143 (MQ=38)
             tACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                        <  1:790356/150‑1 (MQ=38)
             tACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                        <  1:159994/150‑1 (MQ=38)
             tACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                        <  1:604638/150‑1 (MQ=255)
              aCTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                       <  2:636792/150‑1 (MQ=38)
              aCTGATGATAATTCATTATTTCACCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                       >  1:605263/1‑150 (MQ=255)
                tgatgaTAATTCATTATTTCGCCTGGCTTATCAGTCCGTGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcct                                                                                                                                     <  2:1271622/150‑3 (MQ=37)
                tgatgaTAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAt                                                                                                                                     >  1:1170748/1‑150 (MQ=255)
                  atgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAttt                                                                                                                                   <  2:945748/150‑1 (MQ=255)
                   tgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTg                                                                                                                                  >  1:174546/1‑150 (MQ=255)
                    gatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGc                                                                                                                                 <  1:1238021/150‑1 (MQ=255)
                     atAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCt                                                                                                                                >  2:1211338/1‑150 (MQ=255)
                     atAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCt                                                                                                                                <  1:289813/150‑1 (MQ=255)
                     atAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCt                                                                                                                                <  1:722869/150‑1 (MQ=255)
                       aaTTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTa                                                                                                                              >  1:198620/1‑150 (MQ=255)
                       aaTTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTa                                                                                                                              <  2:1093838/150‑1 (MQ=255)
                        aTTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAg                                                                                                                             <  2:1159235/150‑1 (MQ=255)
                        aTTCATTATTTCGCCTGGCGTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTCGCTTAg                                                                                                                             >  2:402757/1‑150 (MQ=255)
                           cATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTc                                                                                                                          >  1:1265966/1‑150 (MQ=255)
                            attattTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCt                                                                                                                                <  1:403423/143‑1 (MQ=255)
                                      cTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGAt                                                                                                               >  2:706298/1‑150 (MQ=255)
                                        ggCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGc                                                                                                                                 >  2:1170828/1‑130 (MQ=255)
                                           ttATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAcc                                                                                                          <  2:263114/150‑1 (MQ=255)
                                           ttATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAcc                                                                                                          <  2:146308/150‑1 (MQ=255)
                                            tATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCt                                                                                                         >  2:303271/1‑150 (MQ=255)
                                             aTCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTa                                                                                                        >  1:675442/1‑150 (MQ=255)
                                                        tgtAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTc                                                                                             >  2:1375562/1‑150 (MQ=255)
                                                        tgtAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTc                                                                                             <  2:1265966/150‑1 (MQ=255)
                                                          tAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATAAGACCTGGAGCAAATGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTccc                                                                                           >  2:151029/1‑150 (MQ=255)
                                                            ggCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCa                                                                                         <  2:594159/150‑1 (MQ=255)
                                                             gCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCaa                                                                                        >  2:853137/1‑150 (MQ=255)
                                                               tGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAAt                                                                                      >  1:1123059/1‑150 (MQ=255)
                                                               tGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAAt                                                                                      >  1:740254/1‑150 (MQ=255)
                                                               tGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAAt                                                                                      >  2:841972/1‑150 (MQ=255)
                                                                  aTCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGtt                                                                                   <  1:314323/150‑1 (MQ=255)
                                                                   tCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGttt                                                                                  <  2:369267/150‑1 (MQ=255)
                                                                    ccGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTa                                                                                 >  1:46232/1‑150 (MQ=255)
                                                                           tCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCTGATAAGGCGTTTACGGCGCATCAGGCATTTTCTTAGCTCAGAGTGAAGATTAATCTATCTGGCCATTCCCTAAATGTTTAAggcggc                                                                          <  2:1155269/150‑1 (MQ=38)
                                                                            cAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACgcggcg                                                                         >  1:1069926/1‑150 (MQ=255)
                                                                              aTGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCt                                                                       <  1:614707/150‑1 (MQ=255)
                                                                                   tGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAAct                                                                  <  2:1218761/150‑1 (MQ=255)
                                                                                     aTGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAActtc                                                                <  2:1123059/150‑1 (MQ=255)
                                                                                         gACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTcc                                                            >  1:1400578/1‑150 (MQ=255)
                                                                                           cGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTg                                                          >  1:232497/1‑150 (MQ=255)
                                                                                                   cgTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAg                                                  <  2:864873/150‑1 (MQ=255)
                                                                                                     tCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTa                                                                                                             >  2:327937/1‑89 (MQ=255)
                                                                                                        tATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTa                                             <  1:1426584/150‑1 (MQ=255)
                                                                                                          tCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaaa                                           >  1:1380477/1‑150 (MQ=255)
                                                                                                           cAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                          >  2:508347/1‑150 (MQ=255)
                                                                                                           cAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                          >  1:423944/1‑150 (MQ=255)
                                                                                                           cAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                          <  1:637852/150‑1 (MQ=255)
                                                                                                           cAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                          <  1:774959/150‑1 (MQ=255)
                                                                                                            aGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGGCTTAACTTCTTCCTGCTTACCAGCGTTAAACgg                                        <  1:429784/150‑1 (MQ=255)
                                                                                                               ccTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGAc                                      >  2:958464/1‑150 (MQ=255)
                                                                                                                cTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACg                                     <  1:1312730/150‑1 (MQ=255)
                                                                                                                cTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACg                                     <  2:675442/150‑1 (MQ=255)
                                                                                                                  ggAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                                   <  1:1268715/150‑1 (MQ=255)
                                                                                                                  ggAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                                   >  1:48258/1‑150 (MQ=255)
                                                                                                                       aaTTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATcc                              >  1:351362/1‑150 (MQ=255)
                                                                                                                       aaTTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATcc                              >  1:66329/1‑150 (MQ=255)
                                                                                                                         ttGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCgg                            <  1:843948/150‑1 (MQ=255)
                                                                                                                         ttGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCgg                            <  1:706298/150‑1 (MQ=255)
                                                                                                                          tGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCTGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCggg                           <  2:740254/150‑1 (MQ=255)
                                                                                                                          tGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCggg                           <  1:750826/150‑1 (MQ=255)
                                                                                                                           gCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGc                          >  2:1297823/1‑150 (MQ=255)
                                                                                                                           cccAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGc                          >  1:325887/2‑150 (MQ=255)
                                                                                                                            ccAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCt                         <  2:610767/150‑1 (MQ=255)
                                                                                                                             cAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTa                        <  2:1029951/150‑1 (MQ=255)
                                                                                                                              aaaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTAc                       <  2:198620/150‑1 (MQ=255)
                                                                                                                              aaaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTAc                       >  1:1143230/1‑150 (MQ=255)
                                                                                                                                 tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTa                    <  1:563377/150‑1 (MQ=255)
                                                                                                                                   tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaaa                  <  2:516976/150‑1 (MQ=255)
                                                                                                                                   tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaaa                  >  2:31298/1‑150 (MQ=255)
                                                                                                                                   aatgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaaa                  >  1:682212/4‑150 (MQ=255)
                                                                                                                                    aGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAat                 >  1:1005572/1‑150 (MQ=255)
                                                                                                                                     ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAata                <  2:446886/150‑1 (MQ=255)
                                                                                                                                     ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAata                >  2:502866/1‑150 (MQ=255)
                                                                                                                                     ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCCCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAata                <  1:394805/150‑1 (MQ=255)
                                                                                                                                      gCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAatat               >  2:778728/1‑150 (MQ=255)
                                                                                                                                       ccGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTCCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATc              >  2:1020796/1‑150 (MQ=255)
                                                                                                                                           aTAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGca          <  2:414494/150‑1 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCTATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTAcc                                                    >  1:91824/1‑104 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACAcg      >  1:856844/1‑150 (MQ=255)
                                                                                                                                                gCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACAcgc     >  1:1318526/1‑150 (MQ=255)
                                                                                                                                                 cGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACGTATCTGACCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGAGCATCCGGGCTAGCTAAATATCCGCACAcgcg    >  2:1322078/1‑150 (MQ=255)
                                                                                                                                                 cGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACAcgcg    <  1:499172/150‑1 (MQ=255)
                                                                                                                                                 cGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACAcgcg    >  1:1149920/1‑150 (MQ=255)
                                                                                                                                                  gTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACACgcgg   <  1:1157101/150‑1 (MQ=255)
                                                                                                                                                   tttACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACACgcggc  <  2:228124/150‑1 (MQ=255)
                                                                                                                                                     tACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATc              <  2:1011768/136‑1 (MQ=255)
                                                                                                                                                     tACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATc              >  1:1011768/1‑136 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACACGCGGC  >  NC_000913/4460335‑4460631

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: