Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,485,380 C→G 5.8% intergenic (+41/‑29) spuA → / → spuC‑I glutamine amidotransferase/polyamine:pyruvate transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,485,3800CG 5.8% 220.3 / 4.1 69intergenic (+41/‑29)spuA/spuC‑Iglutamine amidotransferase/polyamine:pyruvate transaminase
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base G (2/2);  total (35/34)
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.32e-01

GTTCCAGGTCATGGATAACCCCCATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTGGCCGAGAAAGGCCCACGCATCATCA  >  NC_002947/2485234‑2485520
                                                                                                                                                  |                                                                                                                                             
gTTCCAGGTCATGGATAACCCCCATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCgg                                                                                                                                             >  7:32821/1‑149 (MQ=255)
    cAGGTCATGGATAACCCCCATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCGGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTgc                                                                                                                                         >  7:241939/1‑149 (MQ=255)
       gTCATGGATAACCCCCATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTgcgc                                                                                                                                       >  2:40111/1‑148 (MQ=255)
        tCATGGATAACCCCCATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCaa                                                                                                                                     <  7:155151/149‑1 (MQ=255)
             gATAACCCCCATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCCGCGCTgcga                                                                                                                                       <  8:221257/142‑2 (MQ=255)
             gATAACCCCCATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCCGCGCTgcga                                                                                                                                       >  7:221257/1‑141 (MQ=255)
                  cccccATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAg                                                                                                                           >  5:128565/1‑149 (MQ=255)
                   ccccATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGt                                                                                                                          >  4:171536/1‑149 (MQ=255)
                    cccATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGt                                                                                                                          >  6:134625/1‑148 (MQ=255)
                     ccATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAg                                                                                                                           <  8:144631/146‑1 (MQ=255)
                     ccATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAg                                                                                                                           >  7:144631/1‑146 (MQ=255)
                     ccATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTg                                                                                                                         >  4:306584/1‑148 (MQ=255)
                            cTGGCGATCTTCCAGTCCTTTGGGAAGGCGTGCCGGCACCGCTCGGTTTTGCGCCAGAGGCTGTTGCTGTCTGGCTGAGCAGCGGCCTTTCCCACCTGTAAAAATTACCGTCGCAGCGCGGGTGCGCAACTTGATGGTGTGAAGTCTAt                                                                                                                 <  1:79208/149‑1 (MQ=255)
                            cTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAAt                                                                                                                 <  8:7259/149‑1 (MQ=255)
                            cTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAAt                                                                                                                 <  5:134625/149‑1 (MQ=255)
                               gCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGa                                                                                                               <  3:157906/148‑1 (MQ=255)
                                         aGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCCGCGCTgcgagg                                                                                                                                     <  7:78020/116‑4 (MQ=255)
                                         aGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCCGCGCTgcgagg                                                                                                                                     >  8:78020/1‑113 (MQ=255)
                                          gTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGAGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACCGACtt                                                                                                   >  4:147392/1‑149 (MQ=255)
                                           tCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTc                                                                                                  <  3:306584/149‑1 (MQ=255)
                                             cTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAAt                                                                                                                 >  1:178679/1‑132 (MQ=255)
                                             cTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAAt                                                                                                                 <  2:178679/132‑1 (MQ=255)
                                              tttGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATg                                                                                                                              >  3:126503/1‑118 (MQ=255)
                                              tttGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATg                                                                                                                              <  4:126503/118‑1 (MQ=255)
                                                  gCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGAcc                                                                                          >  3:139995/1‑149 (MQ=255)
                                                  gCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGAcc                                                                                          >  8:181326/1‑149 (MQ=255)
                                                    aaGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCt                                                                                        >  8:5929/1‑149 (MQ=255)
                                                            ccGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTgg                                                                                 <  7:5929/148‑1 (MQ=255)
                                                             cGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGGAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGc                                                                                >  1:219209/1‑148 (MQ=255)
                                                               gCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAAt                                                                                                                 >  8:256643/1‑114 (MQ=255)
                                                               gCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAAt                                                                                                                 <  7:256643/114‑1 (MQ=255)
                                                               gCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCaaa                                                                             <  3:63806/149‑1 (MQ=255)
                                                                  ccGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCa                                                                          <  7:181326/149‑1 (MQ=255)
                                                                     ctctGTATTGCGCCATAGTCTGTATCTGTTTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTAGCAGCGCTGCTGCTCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGa                                                                       <  3:48021/149‑1 (MQ=255)
                                                                            ttGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGa                                                                >  2:90741/1‑149 (MQ=255)
                                                                            ttGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGa                                                                >  6:30512/1‑149 (MQ=255)
                                                                               cgcCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGAcc                                                              <  8:32821/148‑1 (MQ=255)
                                                                                gcCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGt                                                                                                             <  3:277827/101‑1 (MQ=255)
                                                                                gcCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGt                                                                                                             >  4:277827/1‑101 (MQ=255)
                                                                                 ccAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTgg                                                                                 >  3:183750/1‑127 (MQ=255)
                                                                                 ccAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTgg                                                                                 <  4:183750/127‑1 (MQ=255)
                                                                                  cAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCa                                                          >  1:216890/1‑149 (MQ=255)
                                                                                   agagGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCAc                                                         >  2:46578/1‑149 (MQ=255)
                                                                                    gagGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGTAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCAcc                                                        >  5:312825/1‑149 (MQ=255)
                                                                                         tGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGcc                                                   >  7:274501/1‑149 (MQ=255)
                                                                                         tGTTGCTGTCTGCCTGAGCAGCGCCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGcc                                                   >  1:275198/1‑149 (MQ=255)
                                                                                                     ccTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGc                                          <  1:40111/146‑1 (MQ=255)
                                                                                                         agcagcGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGt                                     >  5:80661/1‑147 (MQ=255)
                                                                                                           cagcGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGGAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAgg                                 <  2:219209/149‑1 (MQ=255)
                                                                                                            agcGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGc                                <  5:30512/149‑1 (MQ=255)
                                                                                                            agcGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGc                                <  8:274501/149‑1 (MQ=255)
                                                                                                             gcGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACTCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGc                                <  1:90741/149‑1 (MQ=255)
                                                                                                                 ccTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTg                           <  3:171536/149‑1 (MQ=255)
                                                                                                                    ttCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCa                                   <  1:316360/138‑1 (MQ=255)
                                                                                                                    ttCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCa                                   >  2:316360/1‑138 (MQ=255)
                                                                                                                           cTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTGGCCGAGAAAg                 <  8:181811/149‑1 (MQ=255)
                                                                                                                             gTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGc                                          <  5:108293/122‑1 (MQ=255)
                                                                                                                             gTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGc                                          >  6:108293/1‑122 (MQ=255)
                                                                                                                              tAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGa                                                                       <  7:125395/92‑1 (MQ=255)
                                                                                                                              tAAAAAATACCCTCGCAGCGCGGCAGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGa                                                                       >  8:125395/1‑92 (MQ=255)
                                                                                                                                aaaaaTACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTGGCCGAGAAAGGCCCa            <  4:139995/149‑1 (MQ=255)
                                                                                                                                 aaaaTACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCTATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTGGc                         >  1:222794/1‑135 (MQ=255)
                                                                                                                                 aaaaTACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTGGc                         <  2:222794/135‑1 (MQ=255)
                                                                                                                                 aaaaTACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTGGCCGAGAAAGGCCCAc           <  2:216890/149‑1 (MQ=255)
                                                                                                                                  aaaTACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTGGCCGAGAAAGGCCCACg          >  5:45366/1‑149 (MQ=255)
                                                                                                                                   aaTACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGGTCAGCGATGTCAGGCAACTGGCCGAGAAAGGCCCACGc         <  6:80661/149‑1 (MQ=255)
                                                                                                                                    aTACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTGGCCGAGAAAGGCCCACGca        <  6:128565/149‑1 (MQ=255)
                                                                                                                                          tCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTGGCCGAGAAAGGCCCACGCAtcatca  >  5:3778/1‑149 (MQ=255)
                                                                                                                                          tCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTGGCCGAGAAAGGCCCACGCAtcatca  >  8:51721/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GTTCCAGGTCATGGATAACCCCCATTACCTGGCCATCTTCCAGTCCTTTGGCAAGGCCTGCCGGCACCGCTCGGTATTGCGCCAGAGGCTGTTGCTGTCTGCCTGAGCAGCGGCCTTTCCCACCTGTAAAAAATACCCTCGCAGCGCGGCTGCGCAACTTGATGGAGTGAAGTCAATGAGTGAACAGAATTCGCAGACCCTTGCCTGGCAAACCATGAGCCGCGACCATCACCTGGCCCCGTTCAGCGATGTCAGGCAACTGGCCGAGAAAGGCCCACGCATCATCA  >  NC_002947/2485234‑2485520

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: