Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,752 G→T 7.7% intergenic (+57/+37) guaA → / ← PP_1033 glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,178,7520GT 7.7% 241.8 / 7.8 78intergenic (+57/+37)guaA/PP_1033glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (30/42);  new base T (3/3);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTT  >  NC_002947/1178611‑1178896
                                                                                                                                              |                                                                                                                                                 
tCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCgg                                                                                                                                            <  4:301451/149‑1 (MQ=255)
   tATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTAGCACCCGCAGGACGGATACACCAGAAAGCCAGCATGATTGCGGGCt                                                                                                                                         >  6:317516/1‑149 (MQ=255)
     tcatcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCtt                                                                                                                                        <  2:210477/148‑1 (MQ=255)
       atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTc                                                                                                                                      <  2:312544/148‑1 (MQ=255)
        tcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGGGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTg                                                                                                                                    <  5:337235/149‑1 (MQ=255)
        tcaACGAAATCGACGGAATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGAGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTg                                                                                                                                    <  3:234766/149‑1 (MQ=255)
          aaCGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTgg                                                                                                                                   <  2:330774/148‑1 (MQ=255)
          aaCGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGt                                                                                                                                  >  3:223473/1‑149 (MQ=255)
            cGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTc                                                                                                                                <  3:123639/149‑1 (MQ=255)
                aTCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgt                                                                                                                            <  8:114918/149‑1 (MQ=255)
                   gACGGTATCTCGCGCGGTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTCTCGAGTGGGAGGGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttc                                                                                                                          <  3:210005/148‑1 (MQ=255)
                     cGGTATCTCGCGCTTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGTGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttcgtt                                                                                                                       <  8:284160/149‑1 (MQ=255)
                       gTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttcgtt                                                                                                                       >  3:278534/1‑147 (MQ=255)
                         aTCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttcgtt                                                                                                                       >  7:265756/1‑145 (MQ=255)
                            tcGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTg                                                                                                                >  6:97261/1‑149 (MQ=255)
                             cgcgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGc                                                                                                               >  3:260740/1‑149 (MQ=255)
                              gcgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTg                                                                                                                <  6:312805/147‑1 (MQ=255)
                              gcgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGc                                                                                                               <  6:37225/148‑1 (MQ=255)
                              gcgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGc                                                                                                               <  5:300485/148‑1 (MQ=255)
                              gcgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCt                                                                                                              >  8:52367/1‑149 (MQ=255)
                               cgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGACTGCtt                                                                                                            <  4:223473/149‑1 (MQ=255)
                                gcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCGAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTa                                                                                                           <  5:101195/149‑1 (MQ=255)
                                gcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTa                                                                                                           <  7:132945/149‑1 (MQ=255)
                                gcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGAGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTa                                                                                                           <  6:319328/149‑1 (MQ=255)
                                  gTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAtt                                                                                                         <  6:92958/149‑1 (MQ=255)
                                      ccTTCGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAttcttc                                                                                                     <  3:313832/149‑1 (MQ=255)
                                      ccTACGACGTGTCGAGCTAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCGAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAttcttc                                                                                                     <  4:302696/149‑1 (MQ=255)
                                              gtgtCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTc                                                                                             >  8:312976/1‑149 (MQ=255)
                                               tgtCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCt                                                                                            <  4:260740/149‑1 (MQ=255)
                                                 tCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTgg                                                                                          >  8:234031/1‑149 (MQ=255)
                                                   gAGCAAGCCGCCTGCCACTACCGAGTGGGAGTGATCGGGTAGGCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTTTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCa                                                                                        <  7:337537/149‑1 (MQ=255)
                                                       aaGCCGCCTGCGACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCTGGCTGGCACGGGCAGGACGGATACACCAGATAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATa                                                                                    <  1:90214/149‑1 (MQ=255)
                                                       aaGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATa                                                                                    <  5:97261/149‑1 (MQ=255)
                                                         gccgccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAgc                                                                                  <  7:88696/149‑1 (MQ=255)
                                                            gccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCt                                                                               >  3:161773/1‑149 (MQ=255)
                                                             ccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACAAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTTCTTATTCTTCTGGCTTTCTTGCAAATAGCTCtt                                                                              >  8:349711/1‑149 (MQ=255)
                                                              cTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGAGCACCCGCAGGACGGATACACCAGAAAGCCCCCATGATTGCGGCCTTTCTGGTTCATGTTCGTACAGAATCCTAATTCTTCTGGATTAATGCCAAAAAGCGCtt                                                                              >  1:80111/1‑149 (MQ=255)
                                                              cTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCttt                                                                             <  7:52367/149‑1 (MQ=255)
                                                                  cACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTgg                                                                         >  2:85560/1‑149 (MQ=255)
                                                                     tATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAg                                                                      >  1:78304/1‑149 (MQ=255)
                                                                     tATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGGTCTGGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAg                                                                      >  8:221557/1‑149 (MQ=255)
                                                                       tCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCt                                                                                            <  4:112795/125‑1 (MQ=255)
                                                                          aGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGt                                                                  >  2:116307/1‑148 (MQ=255)
                                                                              ggAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAg                                                             <  2:80111/149‑1 (MQ=255)
                                                                               gAGTGATCGGGTAGCTAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttcg                                                                                                                         <  6:137118/89‑1 (MQ=255)
                                                                               gAGTGATCGGGTAGCTAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttcg                                                                                                                         >  5:137118/1‑89 (MQ=255)
                                                                                 gTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaa                                                          <  3:71478/149‑1 (MQ=255)
                                                                                  tGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                         <  1:85560/149‑1 (MQ=255)
                                                                                    aTCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAag                                                        <  4:278534/148‑1 (MQ=255)
                                                                                             ccagccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctc                                              >  8:296828/1‑149 (MQ=255)
                                                                                               agccCAGGCTAGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctcct                                            <  5:317516/149‑1 (MQ=255)
                                                                                                gccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTATATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctcct                                            <  1:116307/149‑1 (MQ=255)
                                                                                                 cccAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATTGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCa                                          >  4:20668/1‑149 (MQ=255)
                                                                                                  ccAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGAGTGCGGGGTTTCTGGGTCTTGGTCGTGCAGCTTGCTTATTCTTCTGGGGTTCTGGCTAATAGCGCGTTCGGGTAGGGGTTATAGGAAAGCTGCCTGGGCCTCAc                                         >  8:137031/1‑149 (MQ=255)
                                                                                                           gCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGa                                 <  2:13254/149‑1 (MQ=255)
                                                                                                            cACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGaa                                <  8:265756/149‑1 (MQ=255)
                                                                                                              cccGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGaa                                <  8:192301/147‑1 (MQ=255)
                                                                                                               ccGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctc                                              >  4:71347/1‑132 (MQ=255)
                                                                                                               ccGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctc                                              <  3:71347/132‑1 (MQ=255)
                                                                                                                  cAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGt                           >  4:123761/1‑148 (MQ=255)
                                                                                                                    ggACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTta                                                                <  7:261239/109‑1 (MQ=255)
                                                                                                                    ggACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTta                                                                >  8:261239/1‑109 (MQ=255)
                                                                                                                    ggACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCAc                        <  6:251002/149‑1 (MQ=255)
                                                                                                                      aCGGATACACCAGAAAGCCCGCATGATTGCGGGCTATCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCAcc                       >  8:355353/1‑148 (MQ=255)
                                                                                                                       aacaagaacCAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCt                     >  1:336563/8‑149 (MQ=38)
                                                                                                                        ggATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCg                            <  3:15141/141‑1 (MQ=255)
                                                                                                                        ggATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCg                            >  4:15141/1‑141 (MQ=255)
                                                                                                                         gATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTc                    >  3:304420/1‑148 (MQ=255)
                                                                                                                          aagaacCAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCa                  >  2:27132/5‑149 (MQ=39)
                                                                                                                          aagaacCAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCa                  >  2:33075/5‑149 (MQ=39)
                                                                                                                              acCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCAcgcg              >  7:303850/1‑149 (MQ=255)
                                                                                                                                 aGAAAGCCCGCATGATTGCGGGCTTTCTGGTTGTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGtg            >  1:59696/1‑147 (MQ=255)
                                                                                                                                  gAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgg          <  2:78304/149‑1 (MQ=255)
                                                                                                                                     aGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggggg       >  6:129854/1‑149 (MQ=255)
                                                                                                                                      gCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGt      >  3:53871/1‑149 (MQ=255)
                                                                                                                                         cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCt   <  7:234031/149‑1 (MQ=255)
                                                                                                                                         cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTgac   <  7:296828/149‑3 (MQ=255)
                                                                                                                                          gCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTgacg  <  2:59696/149‑4 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
TCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTT  >  NC_002947/1178611‑1178896

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