Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 368,813 A→T 6.3% Q68L (CAG→CTG)  mhpA → 3‑(3‑hydroxyphenyl)propionate hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913368,8130AT 6.3% 233.2 / 4.0 73Q68L (CAG→CTG) mhpA3‑(3‑hydroxyphenyl)propionate hydroxylase
Reads supporting (aligned to +/- strand):  ref base A (37/31);  new base T (2/3);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCGTAACGCCTTTATTCAGC  >  NC_000913/368665‑368959
                                                                                                                                                    |                                                                                                                                                  
gCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCa                                                                                                                                                    <  6:138863/149‑1 (MQ=255)
  gagcGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAAGTCGATATGTTGATCGTCTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGGAgt                                                                                                                                                  <  4:321850/146‑1 (MQ=255)
   aTCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGCCGTGCTGGTGGTGGAGAAACTCGATATGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCAGGCAgtc                                                                                                                                                 >  5:244772/1‑149 (MQ=255)
   aTCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACCTGCTGGTGGTGGAGATACTCGATAAGTTGATAGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAgtc                                                                                                                                                 >  2:276471/1‑149 (MQ=255)
          gTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGc                                                                                                                                          <  8:181194/149‑1 (MQ=255)
            gCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTAGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGTCGCTGCGCACCATGCAGTCGGTCGGCCt                                                                                                                                        >  5:189344/1‑149 (MQ=255)
                    ggTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCgat                                                                                                                                 >  3:42171/1‑148 (MQ=255)
                     gTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCgatga                                                                                                                               >  7:200422/1‑149 (MQ=255)
                     gTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCgatga                                                                                                                               >  5:18752/1‑149 (MQ=255)
                        gggCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGGGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGatgatg                                                                                                                             <  5:278213/148‑1 (MQ=255)
                          gCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATAACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTc                                                                                                                          <  2:181719/149‑1 (MQ=255)
                             gatgatGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAGCTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGc                                                                                                                       >  6:179574/1‑149 (MQ=255)
                              atgatgGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCgatg                                                                                                                                >  5:139791/1‑139 (MQ=255)
                              atgatgGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCgatg                                                                                                                                <  6:139791/139‑1 (MQ=255)
                                    gCGATCTTTCTCGTCGAGATGGGCTTGGACGGTCTGGTGGTGGAGAAACTCGATTATTTGATCGACTACGCGCTTGCGATTGGTATTGATGACGTGGGGCGGCGCACCATGGTGTCGGGCGGCTTGGTCGAGTATGTTGTGGCGcact                                                                                                                 <  4:294883/148‑2 (MQ=255)
                                    gCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGcaca                                                                                                                 <  7:256293/148‑1 (MQ=255)
                                      gAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTa                                                                                                              >  3:26628/1‑149 (MQ=255)
                                       aaCTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTAc                                                                                                             <  2:330603/149‑1 (MQ=255)
                                        aCTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACg                                                                                                            >  7:268470/1‑149 (MQ=255)
                                            tctcGGCCAGATGGGCAGTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCg                                                                                                         <  5:274864/148‑1 (MQ=255)
                                                ggCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGca                                                                                                                   <  3:334006/134‑1 (MQ=255)
                                                ggCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGca                                                                                                                   >  4:334006/1‑134 (MQ=255)
                                                ggCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTg                                                                                                       <  8:266726/146‑1 (MQ=255)
                                                    aGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCAcgc                                                                                                 >  8:361906/1‑148 (MQ=255)
                                                     gATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGa                                                                                               >  8:139474/1‑149 (MQ=255)
                                                        gggCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGc                                                                                            >  4:66183/1‑149 (MQ=255)
                                                          gCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGt                                                                                          >  4:261747/1‑149 (MQ=255)
                                                                         tggtggAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGaaa                                                                           <  1:305175/149‑1 (MQ=255)
                                                                           gtggAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAgg                                                                         >  8:191432/1‑149 (MQ=255)
                                                                            tggAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGTTTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGc                                                                        <  4:161587/149‑1 (MQ=255)
                                                                              gagaAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCg                                                                      >  3:320203/1‑149 (MQ=255)
                                                                                   aCTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGcc                                                                       >  4:226718/1‑143 (MQ=255)
                                                                                   aCTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGcc                                                                       <  3:226718/143‑1 (MQ=255)
                                                                                       gATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGc                                                             >  2:391837/1‑149 (MQ=255)
                                                                                       gATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGc                                                             >  7:202803/1‑149 (MQ=255)
                                                                                            gTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGata                                                        >  7:308043/1‑149 (MQ=255)
                                                                                             ttGATCGACTACCCGCGTGCGATTGGTATTGATGCCGAGGCGCTGCGCCCCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGatat                                                       >  3:306286/1‑149 (MQ=255)
                                                                                                aTCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCTGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATAtt                                                      >  8:123289/1‑147 (MQ=255)
                                                                                                aTCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCTGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATAtt                                                      >  8:117637/1‑147 (MQ=255)
                                                                                                aTCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCTGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCa                                                    <  7:117637/149‑1 (MQ=255)
                                                                                                aTCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCTGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCa                                                    <  7:123289/149‑1 (MQ=255)
                                                                                                aTCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCa                                                    <  3:237050/149‑1 (MQ=255)
                                                                                                 tCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAg                                                   <  7:139474/149‑1 (MQ=255)
                                                                                                   gACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGc                                                  <  3:261747/148‑1 (MQ=255)
                                                                                                   gACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGc                                                  <  6:18752/148‑1 (MQ=255)
                                                                                                   gACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGAGGATGTTCTGCCGCACCCTACGCCGTGGCACGCGATGCGGTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGcc                                                 <  1:276471/149‑1 (MQ=255)
                                                                                                       aCCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCa                                                                                   >  8:109991/1‑111 (MQ=255)
                                                                                                       aCCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCa                                                                                   <  7:109991/111‑1 (MQ=255)
                                                                                                        cccGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATAtt                                                      <  4:267755/139‑1 (MQ=255)
                                                                                                        cccGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATAtt                                                      >  3:267755/1‑139 (MQ=255)
                                                                                                        cccGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAAt                                              <  7:191432/147‑1 (MQ=255)
                                                                                                         ccGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGAc                                           <  5:325831/149‑1 (MQ=255)
                                                                                                         ccGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGAc                                           <  4:42171/149‑1 (MQ=255)
                                                                                                                    tGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAAttt                                 <  8:200422/148‑1 (MQ=255)
                                                                                                                      gTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTtgg                               <  4:320203/148‑1 (MQ=255)
                                                                                                                            atgaCGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCcg                        >  8:99297/1‑149 (MQ=255)
                                                                                                                              gaCGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTc                                                     >  4:258841/1‑118 (MQ=255)
                                                                                                                              gaCGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTc                                                     <  3:258841/118‑1 (MQ=255)
                                                                                                                              gaCGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCcgcg                      >  6:343541/1‑149 (MQ=255)
                                                                                                                                cGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGcc                    <  4:26628/149‑1 (MQ=255)
                                                                                                                                 gAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAGAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCg                   >  8:404096/1‑149 (MQ=255)
                                                                                                                                 gAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCg                   >  3:184040/1‑149 (MQ=255)
                                                                                                                                 gAGGCGCTGCGCACCAGGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCg                   >  5:310365/1‑149 (MQ=255)
                                                                                                                                   ggCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCGt                  >  5:163873/1‑148 (MQ=255)
                                                                                                                                      gcTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCGTAACg              >  4:301612/1‑149 (MQ=255)
                                                                                                                                        tGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCGTAACGcc            >  8:201522/1‑149 (MQ=255)
                                                                                                                                          cgcACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCGTAACGCCtt          >  2:310319/1‑149 (MQ=255)
                                                                                                                                            cACCATGCAGTCGGTCTGCCTGGTCGATGATGTTCTGCCTCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCAGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCGTAACGCCTTTa        <  6:244772/149‑1 (MQ=255)
                                                                                                                                               cacgCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCGTAACGCCTTTATTc     <  6:310365/146‑1 (MQ=255)
                                                                                                                                               cATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACTCTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCGTAACGCCTTTAtt      >  3:112984/1‑148 (MQ=255)
                                                                                                                                               cATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATAtt                                                      >  2:292025/1‑100 (MQ=255)
                                                                                                                                               cATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATAtt                                                      <  1:292025/100‑1 (MQ=255)
                                                                                                                                                aTGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCGTAACGCCTTTAtt      >  4:357869/1‑147 (MQ=255)
                                                                                                                                                   cAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCGTAACGCCTTTATTCAGc  <  8:268470/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTCGATAAGTTGATCGACTACCCGCGTGCGATTGGTATTGATGACGAGGCGCTGCGCACCATGCAGTCGGTCGGCCTGGTCGATGATGTTCTGCCGCACACTACGCCGTGGCACGCGATGCGTTTTCTCACCCCGAAAGGCCGCTGTTTTGCTGATATTCAGCCAATGACCGATGAATTTGGCTGGCCGCGCCGTAACGCCTTTATTCAGC  >  NC_000913/368665‑368959

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: