Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,873,917 G→C 7.7% H26D (CAT→GAT) ‡ cysC ← adenosine 5'‑phosphosulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,873,9170GC 7.7% 303.2 / 12.4 93H26D (CAT→GAT) ‡cysCadenosine 5'‑phosphosulfate kinase
Reads supporting (aligned to +/- strand):  ref base G (46/39);  new base C (3/3);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

GCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGT  >  NC_000913/2873774‑2874064
                                                                                                                                               |                                                                                                                                                   
gCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCg                                                                                                                                                 <  3:338340/148‑1 (MQ=255)
    cATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACTCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGATGCCGGTAAACCACAGCACTACACCACGATGACCgtggtg                                                                                                                                            >  2:251079/1‑149 (MQ=255)
    cATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtggt                                                                                                                                             >  8:39109/1‑148 (MQ=255)
    cATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtggt                                                                                                                                             <  7:39109/148‑1 (MQ=255)
      tAATCCGTGGTTAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTa                                                                                                                                          <  4:102961/149‑1 (MQ=255)
      tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTa                                                                                                                                          <  1:29879/149‑1 (MQ=255)
      tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTa                                                                                                                                          <  2:360784/149‑1 (MQ=255)
      tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtggtg                                                                                                                                            >  7:419175/1‑147 (MQ=255)
        aTCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGATGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAg                                                                                                                                         <  1:251079/148‑1 (MQ=255)
           cGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCCCGATGACCGTGGTGTAGCTcg                                                                                                                                     >  3:264910/1‑149 (MQ=255)
            gTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcgc                                                                                                                                    <  2:166033/149‑1 (MQ=255)
            gTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcgc                                                                                                                                    <  8:249571/149‑1 (MQ=255)
             tGGCGAACATTGTCGCCATCCAGCTGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcgcg                                                                                                                                   <  3:240649/149‑1 (MQ=255)
               gCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATAACCGTGGTGTAGCTCGCgtt                                                                                                                                 >  1:365531/1‑149 (MQ=255)
                 gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAgtgctgtg                                                                                                                                  <  7:270547/146‑8 (MQ=255)
                 gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAgtgctgtg                                                                                                                                  >  8:270547/1‑139 (MQ=255)
                      ttGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGca                                                                                                                           >  8:286794/1‑148 (MQ=255)
                        gTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                         <  2:9194/148‑1 (MQ=255)
                         tCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGAAGACCGTGGTGTAGCTCGCGTTGTTGCACAGt                                                                                                                       >  3:38868/1‑149 (MQ=255)
                         tCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCAAGCTCACCGGCGACCGTTGATTTACCGGACCCGGGGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGt                                                                                                                       <  5:291350/149‑1 (MQ=255)
                            ccATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCCGCGCCCCGTCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGAc                                                                                                                    <  6:315042/149‑1 (MQ=255)
                             cATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGAcc                                                                                                                   >  1:235258/1‑149 (MQ=255)
                                    cagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGc                                                                                                            >  8:138704/1‑149 (MQ=255)
                                       aTACGTACTGACGCCGGGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTAt                                                                                                         >  2:266536/1‑149 (MQ=255)
                                        tACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATg                                                                                                        <  5:55030/149‑1 (MQ=255)
                                         aCGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtg                                                                                                                                               >  2:88052/1‑109 (MQ=255)
                                         aCGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtg                                                                                                                                               <  1:88052/109‑1 (MQ=255)
                                          cGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATAACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGc                                                                                                       <  2:365531/148‑1 (MQ=255)
                                              cTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAgac                                                                                                  <  2:63944/149‑1 (MQ=255)
                                                   gCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGCCCGGATGGCTATGCCAGACGACGt                                                                                             >  3:13473/1‑149 (MQ=255)
                                                   gCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGt                                                                                             >  6:136117/1‑149 (MQ=255)
                                                     cGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGttt                                                                                           >  6:282555/1‑149 (MQ=255)
                                                     cGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGttt                                                                                           <  4:239763/149‑1 (MQ=255)
                                                      gAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGtttt                                                                                          >  7:143074/1‑149 (MQ=255)
                                                         tttATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCAGGAGAGGCCGGTAAACCACAGCACTACACCACGATGAACGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGt                                                                                       >  7:395420/1‑149 (MQ=255)
                                                           tATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTc                                                                                      <  1:211541/148‑1 (MQ=255)
                                                           tATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTc                                                                                      <  7:106190/148‑1 (MQ=255)
                                                           tATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTc                                                                                      <  5:232680/148‑1 (MQ=255)
                                                            aTGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCa                                                                                     <  8:419175/148‑1 (MQ=255)
                                                               tAACGCCTCCTCAAGCGCCCCGCCGACCGTTGATTTACCGGACCCGGAGAGGCCGCTAAACCACAGCACTACACCACGATGACCGTCGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCa                                                                                 >  8:330452/1‑149 (MQ=255)
                                                                aaCGCCTCCTCCAGCGCCCCAGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                         <  3:121713/108‑1 (MQ=255)
                                                                aaCGCCTCCTCCAGCGCCCCAGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                         >  4:121713/1‑108 (MQ=255)
                                                                   gCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgcgc                                                                             <  8:396659/149‑1 (MQ=255)
                                                                      tcctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCat                                                                          >  6:112206/1‑149 (MQ=255)
                                                                      tcctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACCTTTTCGTCATGCAGCGCCat                                                                          >  4:172638/1‑149 (MQ=255)
                                                                        ctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatta                                                                        >  3:111405/1‑149 (MQ=255)
                                                                           cAGCGCCCCGGCGACCGTTGCTTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCattatt                                                                      <  2:384129/148‑1 (MQ=255)
                                                                           cAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCattatt                                                                      >  3:95583/1‑148 (MQ=255)
                                                                            aGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttattta                                                                    <  5:136117/149‑1 (MQ=255)
                                                                                 cccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcccc                                                               <  4:135041/149‑1 (MQ=255)
                                                                                   cGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGCCGTTTTCGTCATGCAGCGCCATTATTTATcccc                                                               >  5:139083/1‑147 (MQ=255)
                                                                                   cGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGt                                                                                       >  7:315144/1‑123 (MQ=255)
                                                                                   cGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGt                                                                                       <  8:315144/123‑1 (MQ=255)
                                                                                   cGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCa                                                             <  2:120203/149‑1 (MQ=255)
                                                                                   cGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCa                                                             <  5:112206/149‑1 (MQ=255)
                                                                                          cGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttattt                                                                     <  1:64677/134‑1 (MQ=255)
                                                                                          cGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttattt                                                                     >  2:64677/1‑134 (MQ=255)
                                                                                          cGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATc                                                      <  8:345232/149‑1 (MQ=255)
                                                                                          cGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATc                                                      >  8:135276/1‑149 (MQ=255)
                                                                                           gTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcg                                                     >  5:83865/1‑149 (MQ=255)
                                                                                              gATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgcgc                                                  >  8:359939/1‑149 (MQ=255)
                                                                                              gATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgcg                                                   <  7:135276/148‑1 (MQ=255)
                                                                                              gATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTATCCCCCAGCAAATcgcgc                                                  <  7:286794/148‑1 (MQ=255)
                                                                                               aTTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgcgcg                                                 >  8:316611/1‑149 (MQ=255)
                                                                                               aTTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgcgc                                                  <  4:188738/148‑1 (MQ=255)
                                                                                                tttACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                         >  5:103693/1‑76 (MQ=255)
                                                                                                tttACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                         <  6:103693/76‑1 (MQ=255)
                                                                                                     cGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGccc                                            >  8:249149/1‑148 (MQ=255)
                                                                                                      ggACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTGGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCa                                          <  8:143074/149‑1 (MQ=255)
                                                                                                           cGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCg                                     <  7:316611/149‑1 (MQ=255)
                                                                                                                       tAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcgac                         >  7:21689/1‑149 (MQ=255)
                                                                                                                            cacaGCACTACACCACGATGCCCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAAc                     >  2:131172/1‑148 (MQ=255)
                                                                                                                            cacaGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTAt                                                                                                         <  7:432734/64‑1 (MQ=255)
                                                                                                                            cacaGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTAt                                                                                                         >  8:432734/1‑64 (MQ=255)
                                                                                                                            acgcgagcTACACCACGGTCATCGTGGTGTAGTTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAAc                     >  6:259620/8‑148 (MQ=255)
                                                                                                                              caGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATc                                                                  <  6:50475/101‑1 (MQ=255)
                                                                                                                              caGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATc                                                                  >  5:50475/1‑101 (MQ=255)
                                                                                                                                gCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCACGCAGCGCCATTATTTATGCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCagag                >  7:111548/1‑149 (MQ=255)
                                                                                                                                    tACACCACGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAtt            >  1:210005/1‑149 (MQ=255)
                                                                                                                                    tACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatt                                                                         <  3:350794/88‑1 (MQ=255)
                                                                                                                                    tACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatt                                                                         >  4:350794/1‑88 (MQ=255)
                                                                                                                                    tACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAtt            >  7:399638/1‑149 (MQ=255)
                                                                                                                                     acacCACGCTGACCGTGGTGTAGCTCGCGTAGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAtt            >  3:394845/1‑148 (MQ=255)
                                                                                                                                     acacCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAtt            >  1:333459/1‑148 (MQ=255)
                                                                                                                                       accacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCaa         >  5:10461/1‑149 (MQ=255)
                                                                                                                                         cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGCCCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAtt       >  6:114992/1‑149 (MQ=255)
                                                                                                                                         cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTTGCGACGAACCAGAGCATTCAAtt       >  4:54928/1‑149 (MQ=255)
                                                                                                                                         cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAtt       >  8:152016/1‑149 (MQ=255)
                                                                                                                                         cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAt        <  5:329498/148‑1 (MQ=255)
                                                                                                                                         cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATATATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAt        >  6:329498/1‑148 (MQ=255)
                                                                                                                                          acGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTc      >  8:246574/1‑149 (MQ=255)
                                                                                                                                            gATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCa    <  7:138704/149‑1 (MQ=255)
                                                                                                                                              tGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGt  >  7:270799/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGT  >  NC_000913/2873774‑2874064

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