Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 446,789 A→G 16.8% intergenic (‑123/+26) yajR ← / ← cyoE putative major facilitator superfamily transport protein YajR/heme O synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913446,7890AG16.8% 200.3 / 33.2 89intergenic (‑123/+26)yajR/cyoEputative major facilitator superfamily transport protein YajR/heme O synthase
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base G (5/10);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

ACCTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGT  >  NC_000913/446640‑446936
                                                                                                                                                     |                                                                                                                                                   
aCCTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTa                                                                                                                                                     >  2:1035045/1‑150 (MQ=255)
  cTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGc                                                                                                                                                   >  1:154453/1‑150 (MQ=255)
     gCGTCATTTTATAATTGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCAcc                                                                                                                                                >  1:846095/1‑150 (MQ=255)
       gTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCtt                                                                                                                                              >  1:313694/1‑150 (MQ=255)
       gTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCtt                                                                                                                                              >  1:1500336/1‑150 (MQ=255)
        tCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttt                                                                                                                                             <  1:629223/150‑1 (MQ=255)
         cATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATTCTAAGAAAAAAGGTGCCGTAGCACCtttt                                                                                                                                            <  2:830292/150‑1 (MQ=255)
         cATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCtttt                                                                                                                                            <  1:665615/150‑1 (MQ=255)
         cATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCtttt                                                                                                                                            <  1:41760/150‑1 (MQ=255)
         cATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCtttt                                                                                                                                            <  1:868479/150‑1 (MQ=255)
          aTTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCttttt                                                                                                                                           <  2:501149/150‑1 (MQ=255)
          aTTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCttttt                                                                                                                                           >  2:26470/1‑150 (MQ=255)
           ttttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCtttttt                                                                                                                                          <  1:788499/150‑2 (MQ=255)
            tttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttttttc                                                                                                                                         <  2:1368960/150‑3 (MQ=255)
            tttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttttttc                                                                                                                                         <  1:1629331/150‑3 (MQ=255)
            tttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttttt                                                                                                                                           >  2:56751/1‑148 (MQ=255)
            tttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttttt                                                                                                                                           <  1:56751/148‑1 (MQ=255)
             ttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttttttct                                                                                                                                        <  2:959645/150‑4 (MQ=255)
             ttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttttttct                                                                                                                                        >  1:959645/1‑147 (MQ=255)
             ttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAAt                                                                                                                                        <  2:550566/150‑1 (MQ=255)
              tataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATa                                                                                                                                       <  2:704132/150‑1 (MQ=255)
              tataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATa                                                                                                                                       <  2:21295/150‑1 (MQ=255)
               ataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCGAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATag                                                                                                                                      <  2:1045565/150‑1 (MQ=255)
               ataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttttttctta                                                                                                                                      <  2:553969/150‑6 (MQ=255)
               ataAGCGTTCATTACCCCCTCTCGTTTAACATCGCGACTAGTGTAAAGCGGCAAACGACTGGGGTAAACATGCCGATTTTTTAGTAAATATTTCGAGAGGGCTATGTTTCTAATGCTAAGAAAAAACGGGCCGTAGCACCTTTTTAATag                                                                                                                                      <  2:853545/150‑1 (MQ=255)
                taATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttttttcttag                                                                                                                                     <  2:601758/150‑7 (MQ=255)
                taATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATaga                                                                                                                                     >  2:1425006/1‑150 (MQ=255)
                   tCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATagagag                                                                                                                                  <  1:1012168/150‑1 (MQ=255)
                    cGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAgg                                                                                                                                 <  2:498354/150‑1 (MQ=255)
                    cGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAgg                                                                                                                                 <  1:216597/150‑1 (MQ=255)
                      ttCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGtt                                                                                                                               >  1:992170/1‑150 (MQ=255)
                           tACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGtt                                                                                                                          <  2:1350183/150‑1 (MQ=255)
                            accaccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTa                                                                                                                         <  1:1491360/150‑1 (MQ=255)
                              caccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTAcc                                                                                                                       <  2:313694/150‑1 (MQ=255)
                                ccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACcac                                                                                                                     >  1:1041449/1‑150 (MQ=255)
                                              tCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagccagcag                                                                                                       <  1:1031342/150‑1 (MQ=255)
                                                gcgACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCg                                                                                                     >  1:1227157/1‑150 (MQ=255)
                                                    caaGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATg                                                                                                 >  2:1361012/3‑150 (MQ=255)
                                                    cTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATg                                                                                                 >  1:1585000/1‑150 (MQ=255)
                                                    cTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATg                                                                                                 >  1:1089042/1‑150 (MQ=255)
                                                      aGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCg                                                                                               <  1:48526/150‑1 (MQ=255)
                                                                 tcAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCATTATCACCTTTTTAATAGAGATGTTTTGTTACCAAACAGCAGCCAGCAGCGTATGCGAGTCCGGTAcc                                                                                    >  1:639177/2‑150 (MQ=255)
                                                                  cAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCa                                                                                   >  1:1081215/1‑150 (MQ=255)
                                                                    aaCGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATa                                                                                 <  1:1380450/150‑1 (MQ=255)
                                                                        aaTGGGGTAAAAAAGCGGATTTTTTATTAAATATTTGGAGAGGGTTATGTTTCTAATTCTAAGAAAAAACGTGCCGTAGCAACTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAAt                                                                             <  2:672255/150‑1 (MQ=255)
                                                                           ggggTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCaa                                                                          >  2:756521/1‑150 (MQ=255)
                                                                             ggTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACg                                                                        <  2:992170/150‑1 (MQ=255)
                                                                             ggTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACg                                                                        <  2:1500336/150‑1 (MQ=255)
                                                                             ggTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACg                                                                        >  1:507664/1‑150 (MQ=255)
                                                                              gTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACgg                                                                       <  2:1275504/150‑1 (MQ=255)
                                                                                aaaCATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGAc                                                                     <  2:879884/150‑1 (MQ=255)
                                                                                   cATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGAcatc                                                                  >  2:582375/1‑150 (MQ=255)
                                                                                        ggATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCAc                                                             >  1:212375/1‑150 (MQ=255)
                                                                                              tttAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGagag                                                       <  2:513205/150‑1 (MQ=255)
                                                                                                  gTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCa                                                   <  1:1425006/150‑1 (MQ=255)
                                                                                                   tAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAg                                                  <  2:1089042/150‑1 (MQ=255)
                                                                                                    aaaTATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGt                                                 >  1:858163/1‑150 (MQ=255)
                                                                                                    aaaTATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGt                                                 >  2:1208222/1‑150 (MQ=255)
                                                                                                      atatTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGa                                               <  1:34094/150‑1 (MQ=255)
                                                                                                           tCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCg                                          >  1:1324402/1‑150 (MQ=255)
                                                                                                               gagGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCgatga                                      <  1:1418715/150‑1 (MQ=255)
                                                                                                                agGGATATGTTTCTAATGCTAAGAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATa                                    <  1:1569307/150‑1 (MQ=255)
                                                                                                                 gggATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATa                                    >  2:471603/1‑150 (MQ=255)
                                                                                                                     tatTTTTCTAATGCTAAGAAAATAGGGTGCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCATCAGAGTATGCGTGTCCGGTACCATAAAATCAATGGACATCATCACCGAGAGGGCTGTTATGGAGATGATAgaga                                <  2:639177/150‑1 (MQ=255)
                                                                                                                       tGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCAAACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATCGCGATGATAGAGAAg                              >  2:1214234/1‑150 (MQ=255)
                                                                                                                         tttCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGcc                            <  1:582375/150‑1 (MQ=255)
                                                                                                                             tAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAAc                        <  1:381074/150‑1 (MQ=255)
                                                                                                                             aaaTGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAAc                        <  2:1081215/149‑1 (MQ=255)
                                                                                                                               aTGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAg                      >  1:733024/1‑150 (MQ=255)
                                                                                                                                tGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGc                     >  2:962597/1‑150 (MQ=255)
                                                                                                                                  cTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCtt                   >  1:1117958/1‑150 (MQ=255)
                                                                                                                                  cTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCtt                   >  1:1263933/1‑150 (MQ=255)
                                                                                                                                  cTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCtt                   >  2:1129540/1‑150 (MQ=255)
                                                                                                                                   tAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCg                           <  1:1001508/141‑1 (MQ=255)
                                                                                                                                   tAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCg                           >  2:1001508/1‑141 (MQ=255)
                                                                                                                                     attaaaaaGGTGCTCCGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcg                >  1:1489344/4‑150 (MQ=255)
                                                                                                                                     aGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcg                >  1:306708/1‑150 (MQ=255)
                                                                                                                                     aGAAAAAAGGTGCCGTAGCAACTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcg                >  2:938282/1‑150 (MQ=255)
                                                                                                                                      ttaaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgc               >  1:1416037/3‑150 (MQ=255)
                                                                                                                                      gAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATagag                                 <  2:876526/132‑1 (MQ=255)
                                                                                                                                      gAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATagag                                 >  1:876526/1‑132 (MQ=255)
                                                                                                                                      gAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgc               >  2:669805/1‑150 (MQ=255)
                                                                                                                                        aaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGc                                           >  2:512171/1‑120 (MQ=255)
                                                                                                                                        aaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGc                                           <  1:512171/120‑1 (MQ=255)
                                                                                                                                        aaaaaGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcgc             >  1:195356/1‑150 (MQ=255)
                                                                                                                                            aGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAg         <  1:26470/150‑1 (MQ=255)
                                                                                                                                                gCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTc     <  1:1035048/150‑1 (MQ=255)
                                                                                                                                                  cGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGGGCCCAGATTCTg   >  2:1584730/1‑150 (MQ=255)
                                                                                                                                                  cGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTg   >  1:1353245/1‑150 (MQ=255)
                                                                                                                                                   gTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGACTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGAGATGATAGAGAAGCCGAACAGCATGCGCGCCAAGATTCTgg  >  2:1043960/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACCTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGT  >  NC_000913/446640‑446936

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