Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,575,528 T→A 6.5% pseudogene (1394/2721 nt) yjiV → putative uncharacterized protein YjiV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,575,5280TA 6.5% 242.6 / 5.5 77pseudogene (1394/2721 nt)yjiVputative uncharacterized protein YjiV
Reads supporting (aligned to +/- strand):  ref base T (34/38);  new base A (2/3);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

TTCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAA  >  NC_000913/4575380‑4575675
                                                                                                                                                    |                                                                                                                                                   
ttCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGAtatt                                                                                                                                                    <  2:517267/150‑1 (MQ=255)
   gTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATc                                                                                                                                                 <  2:245565/150‑1 (MQ=255)
   gTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATc                                                                                                                                                 <  2:415982/150‑1 (MQ=255)
   gTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATc                                                                                                                                                 <  2:439440/150‑1 (MQ=255)
       gTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATa                                                                                                                                             <  2:1476007/150‑1 (MQ=255)
       gTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATa                                                                                                                                             >  2:1222708/1‑150 (MQ=255)
       gTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATa                                                                                                                                             <  1:848778/150‑1 (MQ=255)
           cccGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATg                                                                                                                                         <  2:398158/150‑1 (MQ=255)
           cccGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCAACGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATg                                                                                                                                         >  1:750649/1‑150 (MQ=255)
            ccGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATAATATCAATagttgc                                                                                                                                        <  2:279643/150‑6 (MQ=255)
             cGGAATCAAAAGATCCTGAAGATTATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGc                                                                                                                                       <  1:700465/150‑1 (MQ=255)
             cGGAATCAAAAGATCCTGAAGATAATCATAAGAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTTGGTGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATATATggg                                                                                                                                       <  2:1535858/150‑2 (MQ=255)
              ggAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATAATATCAATagttgccg                                                                                                                                      <  1:789982/150‑8 (MQ=255)
                  tCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATTCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGaaa                                                                                                                                  >  1:205653/1‑150 (MQ=255)
                  tCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGaaa                                                                                                                                  <  2:559319/150‑1 (MQ=255)
                   cAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAaag                                                                                                                                 <  2:930325/150‑1 (MQ=255)
                    aaaaGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAaaga                                                                                                                                >  1:36409/1‑150 (MQ=255)
                    aaaaGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATAATATCAATagttg                                                                                                                                         >  2:404715/1‑137 (MQ=255)
                    aaaaGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATAATATCAATagttg                                                                                                                                         <  1:404715/141‑5 (MQ=255)
                      aaGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAaagaag                                                                                                                              <  1:1214078/150‑1 (MQ=255)
                           ccTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATgg                                                                                                                         <  1:898385/150‑1 (MQ=255)
                            cTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGatat                                                                                                                                                     <  2:582703/121‑1 (MQ=255)
                            cTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGatat                                                                                                                                                     >  1:582703/1‑121 (MQ=255)
                            cTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATggg                                                                                                                        <  2:1018552/150‑1 (MQ=255)
                            cTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATggg                                                                                                                        <  2:79234/150‑1 (MQ=255)
                                aGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCa                                                                                                                    >  1:1430091/1‑150 (MQ=255)
                                aGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCa                                                                                                                    >  1:1191551/1‑150 (MQ=255)
                                 gATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCTTCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCaa                                                                                                                   >  1:626184/1‑150 (MQ=255)
                                  aTAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAAc                                                                                                                  <  2:440190/150‑1 (MQ=255)
                                  aTAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAAc                                                                                                                  <  2:1221583/150‑1 (MQ=255)
                                      tCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTAc                                                                                                              >  1:1397995/1‑150 (MQ=255)
                                        aGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCg                                                                                                            >  2:715496/1‑150 (MQ=255)
                                           aCAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGAc                                                                                                         <  2:394892/150‑1 (MQ=255)
                                              gACGCTGGCCTCAAATTCTTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTAc                                                                                                              <  1:1343171/142‑1 (MQ=255)
                                              gACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTAc                                                                                                              >  2:1343171/1‑142 (MQ=255)
                                              gACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGc                                                                                                      >  2:1276766/1‑150 (MQ=255)
                                                cGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATggg                                                                                                                        <  1:1304996/130‑1 (MQ=255)
                                                cGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATggg                                                                                                                        >  2:1304996/1‑130 (MQ=255)
                                                 gCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTg                                                                                                   <  1:690087/150‑1 (MQ=255)
                                                  cTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGc                                                                                                  >  1:870286/1‑150 (MQ=255)
                                                  cTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGc                                                                                                  <  2:319027/150‑1 (MQ=255)
                                                      ccTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTa                                                                                              >  1:123940/1‑150 (MQ=255)
                                                        tCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACt                                                                                            <  1:938985/150‑1 (MQ=255)
                                                          aaaTTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGa                                                                                          <  1:551215/150‑1 (MQ=255)
                                                                   aTGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATg                                                                                 <  1:413881/150‑1 (MQ=255)
                                                                   aTGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATg                                                                                 >  1:157729/1‑150 (MQ=255)
                                                                      gCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCa                                                                              >  2:21415/1‑150 (MQ=255)
                                                                          cGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCg                                                                          >  2:1448646/1‑150 (MQ=255)
                                                                                ccATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTttta                                                                    <  2:260046/150‑1 (MQ=255)
                                                                                 cATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTtttat                                                                   <  2:333517/150‑1 (MQ=255)
                                                                                      tttAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGa                                                                                          >  2:104074/1‑122 (MQ=255)
                                                                                      tttAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGa                                                                                          <  1:104074/122‑1 (MQ=255)
                                                                                          gCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGaaa                                                          >  2:1430290/1‑150 (MQ=255)
                                                                                           cGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAag                                                         <  2:1191551/150‑1 (MQ=255)
                                                                                             aaGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAagag                                                       <  1:1276766/150‑1 (MQ=255)
                                                                                                ttGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATGCATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCat                                                    <  2:1430091/150‑1 (MQ=255)
                                                                                                 tGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCata                                                   >  1:131818/1‑150 (MQ=255)
                                                                                                 tGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTCCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCata                                                   >  1:838470/1‑150 (MQ=255)
                                                                                                   cTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATg                                                 >  1:975716/1‑150 (MQ=255)
                                                                                                    tGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGa                                                <  1:406802/150‑1 (MQ=255)
                                                                                                         cTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACtttt                                           >  2:676071/1‑150 (MQ=255)
                                                                                                                     ttCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACg                               <  1:406625/150‑1 (MQ=255)
                                                                                                                       cAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAACCCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGa                             >  1:182192/1‑150 (MQ=255)
                                                                                                                       cAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCatat                                                  >  2:316292/1‑129 (MQ=255)
                                                                                                                       cAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCatat                                                  <  1:316292/129‑1 (MQ=255)
                                                                                                                       cAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGa                             >  1:1087604/1‑150 (MQ=255)
                                                                                                                       cAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGa                             >  2:1188465/1‑150 (MQ=255)
                                                                                                                       cAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGa                             >  1:1437923/1‑150 (MQ=255)
                                                                                                                           aCCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGcc                         <  2:1397995/150‑1 (MQ=255)
                                                                                                                             cGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTg                       <  2:123940/150‑1 (MQ=255)
                                                                                                                             cGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTg                       <  2:157729/150‑1 (MQ=255)
                                                                                                                             cGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTg                       >  1:22912/1‑150 (MQ=255)
                                                                                                                                 aaTGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGAtt                   <  1:1448646/150‑1 (MQ=255)
                                                                                                                                    gCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGc                >  1:404314/1‑150 (MQ=255)
                                                                                                                                       gCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGCCCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGccag             >  1:517527/1‑148 (MQ=255)
                                                                                                                                        catgtATTGATAATATCAATACATGGCTGAAGGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCgg            >  1:1442856/5‑150 (MQ=255)
                                                                                                                                             aTTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATg                           <  1:833477/130‑1 (MQ=255)
                                                                                                                                             aTTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATg                           >  2:833477/1‑130 (MQ=255)
                                                                                                                                               tGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCaa     >  2:1153970/1‑150 (MQ=255)
                                                                                                                                                gATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGTACCGCTTGCAGTACTGAAACCCGATGGCAACCCTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACCGATGCCTGGATTTGCCCGGTAACCAAt    >  2:1461970/1‑150 (MQ=255)
                                                                                                                                                  tattatCAATACATGGCAGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATaa  <  1:1272808/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCGTCCGTAACCCGGAATCAAAAGATCCTGAAGATAATCAGAACAGACGCTGGCCTCAAATTCGTAATGGCAACGTATCCCATCGTTTAGCGAAGTTGCTGATGCTGGGGGCTGGATTCAAAACCGTCAATGCGGCAACTATTGATATTATCAATACATGGCTGAAAGAAGCATGGGCCCAACTTACCGGACCGCTTGCAGTACTGAAACCCGATGGCAACCGTTTTTATTTACCGAAAGAGCATATGACTTTTTCTTTAATCACGGATGCCTGGATTTGCCCGGTAACCAATAA  >  NC_000913/4575380‑4575675

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: