Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,590,077 G→C 8.2% A42P (GCG→CCG)  PP_2267 → single‑stranded DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,590,0770GC 8.2% 179.9 / 7.3 61A42P (GCG→CCG) PP_2267single‑stranded DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base C (3/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

GCCCAAAGAACGTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAG  >  NC_002947/2589937‑2590220
                                                                                                                                            |                                                                                                                                               
gCCCAAAGAACGTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCtctc                                                                                                                                         <  2:355822/149‑1 (MQ=255)
  ccAAAGAACGTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCt                                                                                                                                       >  8:23868/1‑149 (MQ=255)
  ccAAAGAACGTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGCCGCCtcttct                                                                                                                                       >  6:228432/1‑149 (MQ=255)
   cAAAGAACGTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCct                                                                                                                                            >  8:70235/1‑143 (MQ=255)
   cAAAGAACGTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCct                                                                                                                                            <  7:70235/143‑1 (MQ=255)
   cAAAGAACGTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGCCGCCtcttctt                                                                                                                                      <  5:228432/149‑2 (MQ=255)
        aaCGTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCg                                                                                                                                   >  6:181833/1‑147 (MQ=255)
        aaCGTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCg                                                                                                                                   >  6:175996/1‑147 (MQ=255)
         aCGTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGaa                                                                                                                                 <  2:66210/148‑1 (MQ=255)
           gTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCt                                                                                                                              >  3:180612/1‑149 (MQ=255)
           gTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGCCACTCACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCATGACCCAAGAATAGGCGGCGCCTCTCCTGGCGAAGCt                                                                                                                              >  6:129482/1‑149 (MQ=255)
                     ccTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCC‑AGGAGAGGCGCCGCCtcttct                                                                                                                                       >  3:141675/1‑129 (MQ=255)
                     ccTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCC‑AGGAGAGGCGCCGCCtcttct                                                                                                                                       <  4:141675/129‑1 (MQ=255)
                      cTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTg                                                                                                                    >  6:4684/1‑148 (MQ=255)
                           cgcgTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGaa                                                                                                              <  5:4684/149‑1 (MQ=255)
                            gcgTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAg                                                                                                             >  1:359829/1‑149 (MQ=255)
                              gTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAgcc                                                                                                           <  3:139499/149‑1 (MQ=255)
                               tGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCTATGAGTACCCGAAGGCGGACTGTTAGTACAAAGTCACCCTGATCCTGACCCAAGAAGACGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAgccg                                                                                                          <  2:318332/149‑1 (MQ=255)
                               tGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCATGAAgccg                                                                                                          >  2:122655/1‑149 (MQ=255)
                                ggCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGATGCACGAAgccgc                                                                                                         >  7:330703/1‑149 (MQ=255)
                                 gCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTcc                                                                                                                                        >  1:139503/1‑117 (MQ=255)
                                 gCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTcc                                                                                                                                        <  2:139503/117‑1 (MQ=255)
                                 gCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAgccgc                                                                                                         >  5:258159/1‑148 (MQ=255)
                                   aTTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGc                                                                                                                               >  6:19937/1‑124 (MQ=255)
                                   aTTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGc                                                                                                                               <  5:19937/124‑1 (MQ=255)
                                              cccGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGcc                                                                                           <  4:180612/149‑1 (MQ=255)
                                                   cACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGaa                                                                                      <  6:323594/149‑1 (MQ=255)
                                                        aCCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGtt                                                                                 >  2:48006/1‑149 (MQ=255)
                                                          cAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTc                                                                                >  4:27764/1‑148 (MQ=255)
                                                          cAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAACCCGCCATCGAGAATGGCCGTGAAGAGTTc                                                                                >  6:278432/1‑148 (MQ=255)
                                                             aCCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCaaga                                                                            <  5:175996/149‑1 (MQ=255)
                                                             aCCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCaaga                                                                            <  5:181833/149‑1 (MQ=255)
                                                               cTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCaagaag                                                                          <  3:27764/149‑1 (MQ=255)
                                                                tGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGc                                                                         >  1:128727/1‑149 (MQ=255)
                                                                 gACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCt                                                                        >  6:329507/1‑149 (MQ=255)
                                                                   attCGGGACCGTTGAGTACCCGTAGCCGGACGGTGAGTACAAAGTCACCCTTATCCTGACTCAAGTAGAGGCGGCGCCTCTCCTGGCTAAGCTGGTACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCa                                                                      <  5:129482/148‑1 (MQ=255)
                                                                      cGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAg                                                                    <  2:180461/148‑1 (MQ=255)
                                                                       ggCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGa                                                                   <  2:359829/148‑1 (MQ=255)
                                                                           ccGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGaa                                                              <  7:15670/149‑1 (MQ=255)
                                                                            cGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGaa                                                                                      >  4:279506/1‑124 (MQ=255)
                                                                            cGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGaa                                                                                      <  3:279506/124‑1 (MQ=255)
                                                                             gATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGc                                                            >  4:155231/1‑149 (MQ=255)
                                                                                      ccGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTAAAGATCGAAGCCCGCaag                                                     >  2:341470/1‑147 (MQ=255)
                                                                                        gAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGc                                                 >  8:37770/1‑149 (MQ=255)
                                                                                        gAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGc                                                 <  6:258159/149‑1 (MQ=255)
                                                                                        gAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGc                                                 >  8:289277/1‑149 (MQ=255)
                                                                                                aCGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGATGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGAc                                         <  8:330703/149‑1 (MQ=255)
                                                                                                    tGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTaa                                     >  1:251764/1‑149 (MQ=255)
                                                                                                       gTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCg                                  >  5:13236/1‑149 (MQ=255)
                                                                                                          cAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTg                                >  7:29444/1‑148 (MQ=255)
                                                                                                           aaaGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTg                                >  4:326132/1‑147 (MQ=255)
                                                                                                                  cccTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGctc                       <  1:48006/149‑1 (MQ=255)
                                                                                                                      gATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCa                   >  4:230618/1‑149 (MQ=255)
                                                                                                                       aTCCTGACCCTAGAAGAGGCGGCGGCTCTCCTGGCGAAGCTGGAAGCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGATAGTCGTAAGCGTGCAGCCGCTCTTc                    <  1:216861/147‑1 (MQ=255)
                                                                                                                             aCCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCgag              >  3:133043/1‑147 (MQ=255)
                                                                                                                             aCCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCgag              <  7:37770/147‑1 (MQ=255)
                                                                                                                               ccAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCgaggag           >  5:279224/1‑148 (MQ=255)
                                                                                                                               ccAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCgaggag           >  4:364988/1‑148 (MQ=255)
                                                                                                                                 caggagagGCGCCGCCTCTCCTGGCGAAGCTGGAACCGCGGCACGAGGCCGCCATCGAGAAGGGCCGTGAAGAGTTCAAGAAGCTCAAGAGCGAAGCCCGCCAGAAGCTGAAAGAGGTATGCGTGAAGCGGCTCGTCAGAGAGGAGTAc        >  8:313171/5‑149 (MQ=255)
                                                                                                                                       agGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGTAGCTGAAGATCGAAGACCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACaac  <  5:278432/149‑2 (MQ=255)
                                                                                                                                       agGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAg  <  8:29444/149‑1 (MQ=255)
                                                                                                                                         gcggcgCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGACAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAg  >  5:39364/1‑147 (MQ=255)
                                                                                                                                            |                                                                                                                                               
GCCCAAAGAACGTAACCATGACCTCGCCGCGTGGCATTTTCCGCTACCCGGCACTGACCAAACCTGACTTCGGCACCGATGAGTACCCGAAGCCGGACGGTGAGTACAAAGTCACCCTGATCCTGACCCAAGAAGAGGCGGCGCCTCTCCTGGCGAAGCTGGAACCGCTGCACGAAGCCGCCATCGAGAATGGCCGTGAAGAGTTCAAGAAGCTCAAGATCGAAGCCCGCAAGAAGCTGAAAGACGTAACCGTGCAGCCGCTCTTCACCGAGGAGTACGACAAG  >  NC_002947/2589937‑2590220

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: ‑