Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,901,761 A→T 6.3% D771V (GAC→GTC)  PP_1703 → assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,901,7610AT 6.3% 180.8 / 3.0 64D771V (GAC→GTC) PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base T (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

CCCTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCATGGCATCAG  >  NC_002947/1901617‑1901909
                                                                                                                                                |                                                                                                                                                    
cccTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAg                                                                                                                                                  <  3:202308/149‑1 (MQ=255)
 ccTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAgg                                                                                                                                                 >  8:219417/1‑149 (MQ=255)
   ttGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGCGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGcc                                                                                                                                               >  4:180537/1‑149 (MQ=255)
     gATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGccc                                                                                                                                              <  3:45624/148‑1 (MQ=255)
         tGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGTCCAGGcccagcag                                                                                                                                         <  8:81008/149‑6 (MQ=255)
          gAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGcc                                                                                                                                         >  7:22842/1‑148 (MQ=255)
           aaCACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCgcgc                                                                                                                                          <  4:79673/146‑1 (MQ=255)
               cagCCACTGTAGAGCCTCTCCGCATCCCGACAGCCACCCTGTACCGTCTGCTAGCCCTGCACACCTTGCCGGTTCCGGTGCTGGCAGAAGAAGAGCTTCATTACCTCAAAGGATTCCTGCTGGGCCTGGACCAGGCACTCGCCGAggg                                                                                                                                    >  3:220097/3‑148 (MQ=255)
                   cACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGACATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTAc                                                                                                                               >  7:218762/1‑149 (MQ=255)
                    aCTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTAcc                                                                                                                              >  5:237977/1‑149 (MQ=255)
                          gagCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCt                                                                                                                        >  1:322289/1‑149 (MQ=255)
                              cTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCg                                                                                                                    <  7:219417/149‑1 (MQ=255)
                              cTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCg                                                                                                                    <  8:22842/149‑1 (MQ=255)
                               tGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCgg                                                                                                                   >  2:145355/1‑149 (MQ=255)
                                gTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGc                                                                                                                  <  7:78975/149‑1 (MQ=255)
                                   cGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAg                                                                                                               >  5:286755/1‑149 (MQ=255)
                                      aTGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTgc                                                                                                            <  6:286755/149‑1 (MQ=255)
                                            aCCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGACATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCgct                                                                                                      <  8:218762/149‑1 (MQ=255)
                                              cGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCGTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCgctgg                                                                                                    <  1:50736/149‑1 (MQ=255)
                                                       gTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGTCAGcc                                                                                           >  8:66534/1‑149 (MQ=255)
                                                               tGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCa                                                                                              <  4:114051/138‑1 (MQ=255)
                                                               tGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCa                                                                                              >  3:114051/1‑138 (MQ=255)
                                                                 cTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgtttgt                                                                                  <  4:15780/148‑1 (MQ=255)
                                                                  ttGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgtttgttt                                                                                >  3:64273/1‑149 (MQ=255)
                                                                     gCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTc                                                                             >  3:106194/1‑149 (MQ=255)
                                                                     gCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTc                                                                             >  6:276309/1‑149 (MQ=255)
                                                                       cggCACACCTTGCCGGGTCCGGTGCGGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGa                                                                           <  8:238408/147‑1 (MQ=255)
                                                                       cTGCACACCTTGCCGGTTCCGGTGCTGGCCGAATAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGa                                                                           <  3:132479/149‑1 (MQ=255)
                                                                          cacaCCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgttt                                                                                    <  7:193999/137‑1 (MQ=255)
                                                                          cacaCCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgttt                                                                                    >  8:193999/1‑137 (MQ=255)
                                                                          cacaCCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACgg                                                                        >  2:244374/1‑149 (MQ=255)
                                                                             acCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGTCCAGGccc                                                                                                                                              <  6:276877/76‑1 (MQ=255)
                                                                             acCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGTCCAGGccc                                                                                                                                              >  5:276877/1‑76 (MQ=255)
                                                                             acCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAg                                                                                             >  3:6696/1‑125 (MQ=255)
                                                                             acCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAg                                                                                             <  4:6696/125‑1 (MQ=255)
                                                                             acCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCCCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGc                                                                       <  5:44043/147‑1 (MQ=255)
                                                                             acCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCCCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGc                                                                       >  6:44043/1‑147 (MQ=255)
                                                                               cTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGa                                                                           >  5:139958/1‑141 (MQ=255)
                                                                               cTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGa                                                                           <  6:139958/141‑1 (MQ=255)
                                                                                        tCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTgcg                                                                                                           >  7:162452/1‑100 (MQ=255)
                                                                                        tCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTgcg                                                                                                           <  8:162452/100‑1 (MQ=255)
                                                                                             ggCTGGCCGAAGAAGAGCGTCATTACCTCCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTc                                                     <  2:281050/148‑1 (MQ=255)
                                                                                                 ggCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGccc                                                 <  5:276309/149‑1 (MQ=255)
                                                                                                    cGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGTCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGc                                              <  7:66534/149‑1 (MQ=255)
                                                                                                      aagaagAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCt                                            >  6:222079/1‑149 (MQ=255)
                                                                                                        gaagaGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGc                                          <  3:180537/149‑1 (MQ=255)
                                                                                                           gagCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAAc                                       >  3:319409/1‑149 (MQ=255)
                                                                                                             gCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTg                                     >  5:63686/1‑149 (MQ=255)
                                                                                                               gTCATTACCTGCAAGGTTTCCTGCTGGGCCGGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCa                                   <  1:244374/149‑1 (MQ=255)
                                                                                                                tCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGGCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCAc                                  >  6:204085/1‑149 (MQ=255)
                                                                                                                  aTTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCAccc                                <  4:64273/149‑1 (MQ=255)
                                                                                                                   ttACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCa                               >  6:161415/1‑149 (MQ=255)
                                                                                                                   ttACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCa                               >  3:237983/1‑149 (MQ=255)
                                                                                                                       cTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATcgc                           >  6:60124/1‑149 (MQ=255)
                                                                                                                        tGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATcgcg                          >  6:211663/1‑149 (MQ=255)
                                                                                                                         gCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCgg                         >  1:36068/1‑149 (MQ=255)
                                                                                                                          cAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCgg                         >  7:202559/1‑148 (MQ=255)
                                                                                                                          cAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCgg                         <  2:36068/148‑1 (MQ=255)
                                                                                                                           aaGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGca                       >  2:26709/1‑149 (MQ=255)
                                                                                                                            tggTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGGGCCGCTGTCTGCCAGCCGTCGTTTGTTTGGCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGTcac                      <  6:237977/148‑1 (MQ=255)
                                                                                                                                      ctgGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCgcgc             <  4:319409/148‑1 (MQ=255)
                                                                                                                                       cgggCCTGGTCCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCa           >  7:15060/2‑149 (MQ=255)
                                                                                                                                          gCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCATgg        <  5:161415/149‑1 (MQ=255)
                                                                                                                                               gACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCATGGCATCa   <  4:106194/149‑1 (MQ=255)
                                                                                                                                                aCCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCATGGCATCAg  <  5:211663/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCCTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAAGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCATGGCATCAG  >  NC_002947/1901617‑1901909

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