Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,205,616 C→A 7.2% V273L (GTG→TTG)  PP_4586 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,205,6160CA 7.2% 198.3 / 4.3 64V273L (GTG→TTG) PP_4586hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/24);  new base A (2/3);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

TTGCTCAGACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGTTACGCGGGAA  >  NC_002947/5205471‑5205763
                                                                                                                                                 |                                                                                                                                                   
ttgctCAGACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCaa                                                                                                                                                  >  5:118691/1‑149 (MQ=255)
    tCAGACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTcg                                                                                                                                              >  5:250973/1‑149 (MQ=255)
      aGACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTcgc                                                                                                                                             <  8:72123/148‑1 (MQ=255)
       gACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAATTGGTGCAACTGTTCgcgc                                                                                                                                           <  6:274870/149‑1 (MQ=255)
        aCGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAATTGGTGCAACTGTTCGCgca                                                                                                                                          <  5:208813/149‑1 (MQ=255)
          gcgcCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAATTGGTGCAACTGTTCGCgcaga                                                                                                                                        <  5:64834/149‑2 (MQ=255)
                               ggTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCg                                                                                                                   >  3:168762/1‑149 (MQ=255)
                                   atcaCTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCg                                                                                                                   >  8:50613/1‑145 (MQ=255)
                                   atcaCTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCg                                                                                                                   <  7:50613/145‑1 (MQ=255)
                                    tcaCTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGAc                                                                                                              <  6:122747/149‑1 (MQ=255)
                                      aCTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCCGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGCCCa                                                                                                            >  6:311293/1‑149 (MQ=255)
                                        tCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATg                                                                                                          >  5:88307/1‑149 (MQ=255)
                                        tCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATg                                                                                                          >  8:132928/1‑149 (MQ=255)
                                        tCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGAGACCAGTTGAACCAATTCGCGCAGCAACGGGCGATTGCCGTTTTCGCTGACCATg                                                                                                          >  5:154083/1‑149 (MQ=255)
                                          aTCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCCCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATggg                                                                                                        >  8:95675/1‑149 (MQ=255)
                                              cgcccgATCGGCGTAACAACTAGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACTCCACTCAATCTGCGCGAACAGTTGCACCAATTAGCGCAGCAACGGGCTATTGCCGTTTTAGATGACAATGGGCCgg                                                                                                    >  6:158583/1‑149 (MQ=255)
                                              cgcccgAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCgg                                                                                                    >  5:83686/1‑149 (MQ=255)
                                                         cgTAACAACTCGTTGGCGCAATAATCGGTCAGGCTCTCTACCAGCCATTCACTGCCGTCACGCCCGGCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGTTCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAg                                                                                         <  6:229686/149‑1 (MQ=255)
                                                               aacTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACaa                                                                                   <  6:83686/149‑1 (MQ=255)
                                                                    gTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGc                                                                              >  6:22288/1‑149 (MQ=255)
                                                                                  cGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCa                                                                  >  3:189225/1‑147 (MQ=255)
                                                                                         gCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTc                                                                                                                    <  4:215817/90‑1 (MQ=255)
                                                                                         gCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTc                                                                                                                    >  3:215817/1‑90 (MQ=255)
                                                                                         gCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTc                                                         >  7:23020/1‑149 (MQ=255)
                                                                                         gCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCCCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGAACTCTTCACGTACAACGAGCCATGCGAAGCCATAGCCCCTc                                                         >  5:38503/1‑149 (MQ=255)
                                                                                             tctACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGGGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGcc                                                      <  6:250973/148‑1 (MQ=255)
                                                                                             tctACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGc                                                       <  7:132928/147‑1 (MQ=255)
                                                                                                aCCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACAt                                                  >  2:167638/1‑149 (MQ=255)
                                                                                                 ccagccagTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACAt                                                  <  4:189225/148‑1 (MQ=255)
                                                                                                 ccagccagTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACAt                                                  <  1:223213/148‑1 (MQ=255)
                                                                                                  cagccagTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGc                                                >  6:308016/1‑149 (MQ=255)
                                                                                                    gccagTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCa                                                                      >  2:320933/1‑125 (MQ=255)
                                                                                                    gccagTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCa                                                                      <  1:320933/125‑1 (MQ=255)
                                                                                                    gccagTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCg                                              <  6:88307/149‑1 (MQ=255)
                                                                                                     ccagTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGt                                             >  6:46494/1‑149 (MQ=255)
                                                                                                     ccagTCACTGCCGTCACGCCCCACAATCTCCGAGAACAGTTGCAACAATTCACGCAGCAACCGGCTATTGACGTATTCGCTGACCATGGGACTGGCACTATTCAGGTAAAACGAGCCCTGCGCAGCCATAGCCCCTCGACACATGCCGt                                             >  6:105830/1‑149 (MQ=255)
                                                                                                           aCTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCAcgcg                                       >  3:80055/1‑149 (MQ=255)
                                                                                                             tGCCGCCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGATGACCATGGGCCGGGCACAGTTCAGGTACAACGAGCCCTCCGAAGCCATACCCCCTCGCCACATGCCGTCACGCGcc                                     >  2:121835/1‑149 (MQ=255)
                                                                                                                cGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCTGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGcccc                                   >  8:113872/1‑148 (MQ=255)
                                                                                                                cGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGcccc                                   <  5:105830/148‑1 (MQ=255)
                                                                                                                    aCGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGAcc                               <  7:95675/148‑1 (MQ=255)
                                                                                                                      gCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGAcc                               >  6:101236/1‑146 (MQ=255)
                                                                                                                       ccccTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCa                              <  8:281836/146‑1 (MQ=255)
                                                                                                                       ccccTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCa                              >  7:281836/1‑146 (MQ=255)
                                                                                                                        cccTCAATCTGCGCTAACAGTTGCACCAATTAGCTCAGCAACGGGGTATTGCCTTTTTCCCTGACCATGGGCCGGTCACAGTTCAGGTACAACTAGCCCTGCGCAGCCATTTCCCCTCTCAACATGCCGTCACGCGCCCCGACCAagag                          <  6:154083/149‑1 (MQ=255)
                                                                                                                        cccTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGa                                 <  7:106598/142‑1 (MQ=255)
                                                                                                                        cccTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGa                                 >  8:106598/1‑142 (MQ=255)
                                                                                                                        cccTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAagag                          >  4:172591/1‑149 (MQ=255)
                                                                                                                         ccTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGCGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGagc                         <  8:157978/149‑1 (MQ=255)
                                                                                                                          cTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAagag                          >  2:29314/1‑147 (MQ=255)
                                                                                                                                 tGCGCGAATTGGTGCAACTGTTCGCGCGGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTggg                 >  5:276116/1‑149 (MQ=255)
                                                                                                                                  gcgcGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGc                >  3:79490/1‑149 (MQ=255)
                                                                                                                                     cgAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCggg             >  3:72732/1‑149 (MQ=255)
                                                                                                                                     cctaCAGTTTGACCAATTCGCGCTGCAAGGGGCTATTGCCGTTTTCGCTGACGATGGGCCGGGCACTGTTGAGGTAGAACGAGCCCTGCGCTGCCTTATCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCggg             <  6:38503/146‑1 (MQ=255)
                                                                                                                                      gAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGt            >  1:304517/1‑149 (MQ=255)
                                                                                                                                        aCAGTTGCACCTATTCGCTCAGCAAAGGGCTATTGCAGTTTTCGATGACCATGTGCCGGGCACTGTTCATGTACAACGATCCCTGCGCTGCCATAGCCCCTCGCCACATGCAGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGt            <  5:158583/147‑1 (MQ=255)
                                                                                                                                          aGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCcac                                                    >  5:315527/1‑105 (MQ=255)
                                                                                                                                          aGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCcac                                                    <  6:315527/105‑1 (MQ=255)
                                                                                                                                           gTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGTTAcgc       <  3:172591/149‑1 (MQ=255)
                                                                                                                                           gTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGTTAcgc       <  6:118691/149‑1 (MQ=255)
                                                                                                                                               cACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGTTACGCggg    >  8:43051/1‑148 (MQ=255)
                                                                                                                                               cACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGTTACGCggg    >  3:86955/1‑148 (MQ=255)
                                                                                                                                               cACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGTTACGCggg    <  5:22288/148‑1 (MQ=255)
                                                                                                                                                aCCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGTTACGCGGGaa  >  4:113719/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTGCTCAGACGCGCCTGCCAGACCTGGTAGCGGTCATCACTCATCCCGCCCGAGCGGCGTAACAACTCGTTGGCGTAATAATCGGTCAGGCTCTCTACCAGCCAGTCACTGCCGTCACGCCCCTCAATCTGCGCGAACAGTTGCACCAATTCGCGCAGCAACGGGCTATTGCCGTTTTCGCTGACCATGGGCCGGGCACTGTTCAGGTACAACGAGCCCTGCGCAGCCATAGCCCCTCGCCACATGCCGTCACGCGCCCCGACCAAGAGCAGCTTGGGCGGGTTACGCGGGAA  >  NC_002947/5205471‑5205763

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