Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,446,538 G→C 6.0% L357L (CTC→CTG nicC ← 6‑hydroxynicotinate 3‑monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,446,5380GC 6.0% 369.0 / 12.5 116L357L (CTC→CTGnicC6‑hydroxynicotinate 3‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (51/58);  new base C (5/2);  total (56/60)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

Only 100 of 123 total aligned reads displayed.

CTTGCCGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGCCATG  >  NC_002947/4446389‑4446685
                                                                                                                                                     |                                                                                                                                                   
cTTGCCGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTag                                                                                                                                                     >  2:1629287/1‑150 (MQ=255)
 ttGCCGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTaga                                                                                                                                                    <  2:971909/150‑1 (MQ=255)
 ttGCCGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTaga                                                                                                                                                    <  1:1153775/150‑1 (MQ=255)
      gccATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAg                                                                                                                                               <  1:737255/150‑1 (MQ=255)
      gccATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAg                                                                                                                                               >  1:571067/1‑150 (MQ=255)
        cATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGgt                                                                                                                                             >  2:808157/1‑150 (MQ=255)
        cATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGgt                                                                                                                                             <  2:1180384/150‑1 (MQ=255)
        aaTAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGgt                                                                                                                                             <  2:1180323/149‑1 (MQ=255)
         aTAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGgtg                                                                                                                                            >  2:55235/1‑150 (MQ=255)
         aTAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTCTACAGCCAggaa                                                                                                                                            <  2:649625/150‑3 (MQ=255)
                  gTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTg                                                                                                                                   <  1:837492/150‑1 (MQ=255)
                       gCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATAc                                                                                                                              <  1:2097376/150‑1 (MQ=255)
                         tGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTaga                                                                                                                                                    >  1:1252547/1‑126 (MQ=255)
                         tGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTaga                                                                                                                                                    <  2:1252547/126‑1 (MQ=255)
                                   cATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTg                                                                                                                  >  1:352308/1‑150 (MQ=255)
                                   cATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTg                                                                                                                  <  2:157493/150‑1 (MQ=255)
                                     ttGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCg                                                                                                                >  2:1220895/1‑150 (MQ=255)
                                        gcgcTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCGGCTGCACATAGAGACCGAACCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATGCCGATTGCACCTGCGATg                                                                                                             >  2:1263518/1‑150 (MQ=255)
                                         cgcTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGc                                                                                                            <  1:412022/150‑1 (MQ=255)
                                         cgcTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGGAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGTGTTGGATACCGATTGCACCTGCGATGc                                                                                                            <  1:1620257/150‑1 (MQ=255)
                                         cgcTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCAGAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGc                                                                                                            <  1:996454/150‑1 (MQ=255)
                                          gcTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGcc                                                                                                           >  2:1680452/1‑150 (MQ=255)
                                                gTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGTCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCTTTGGATACCGATTGCGCCTGCGATGCCCGctct                                                                                                     <  2:492025/150‑1 (MQ=255)
                                                   accaccAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTg                                                                                                  <  2:830910/150‑1 (MQ=255)
                                                     caccaGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCg                                                                                                >  2:877353/1‑150 (MQ=255)
                                                     caccaGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAAGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGct                                                                                                >  2:288610/1‑149 (MQ=255)
                                                        caGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTaga                                                                                                                                                    >  2:573020/1‑95 (MQ=255)
                                                        caGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTaga                                                                                                                                                    <  1:573020/95‑1 (MQ=255)
                                                        caGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGc                                                                                                            <  1:1495793/135‑1 (MQ=255)
                                                        caGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGc                                                                                                            >  2:1495793/1‑135 (MQ=255)
                                                        caGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGtt                                                                                             <  1:1755415/150‑1 (MQ=255)
                                                           cgacgaAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGc                                                                                          <  2:2082188/150‑1 (MQ=255)
                                                                 aGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATa                                                                                    <  1:1078005/150‑1 (MQ=255)
                                                                 aGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATa                                                                                    >  2:2321781/1‑150 (MQ=255)
                                                                   gTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACa                                                                                  >  1:553530/1‑150 (MQ=255)
                                                                      cTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATAc                                                                                   <  1:1895415/146‑1 (MQ=255)
                                                                         aTGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACTGCTTTCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGa                                                                            >  1:1452696/1‑150 (MQ=255)
                                                                         aTGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGa                                                                            >  1:208928/1‑150 (MQ=255)
                                                                         aTGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGa                                                                            >  1:1129551/1‑150 (MQ=255)
                                                                             cgccTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTg                                                                                            <  1:1430586/130‑1 (MQ=255)
                                                                             cgccTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTg                                                                                            >  2:1430586/1‑130 (MQ=255)
                                                                              gccTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGcc                                                                       <  1:1680452/150‑1 (MQ=255)
                                                                                  ccccGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGTTCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGAGCTCGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGc                                                                   <  1:1023722/150‑1 (MQ=255)
                                                                                   cccGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCATGCg                                                                  <  1:2413549/150‑1 (MQ=255)
                                                                                    ccGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCgg                                                                 >  1:609043/1‑150 (MQ=255)
                                                                                     cGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCggt                                                                <  1:2185740/150‑1 (MQ=255)
                                                                                      gCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGTTATTCTTGGTTGTAGAGCCAGGTGTTGGCGTTAGAGACCGATTGTATGTGGCAGATCGGCTCTTTGCGGTTGGCTTAATACAACGCGAACTCCGTGCggtg                                                               >  2:1146299/1‑150 (MQ=255)
                                                                                      gCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCggtg                                                               <  1:72280/150‑1 (MQ=255)
                                                                                          tCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCg                                                           <  2:650004/150‑1 (MQ=255)
                                                                                           cAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCTCGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGGCgg                                                          >  1:2132867/1‑150 (MQ=255)
                                                                                            aaGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATCCACCGAGAGCGCCGAGCggcggttgtc                                                         >  2:2352511/1‑146 (MQ=255)
                                                                                             aGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGAc                                                        <  1:1629287/150‑1 (MQ=255)
                                                                                              gctgctGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACTCGAACGCCGTGCGGTGGTCGGACa                                                       <  2:1622280/150‑1 (MQ=255)
                                                                                              gctgctGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTAATACAACGCGAACGCCGTGCGGTGGTCGGACa                                                       >  1:1722466/1‑150 (MQ=255)
                                                                                                 gctgcCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGcc                                                    <  2:553530/150‑1 (MQ=255)
                                                                                                   tgcCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGTCAGGTGTTGGTTTGGCATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCg                                                  <  1:877353/150‑1 (MQ=255)
                                                                                                   tgcCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCg                                                  <  2:1080730/150‑1 (MQ=255)
                                                                                                   tgcCCATAGAGGGCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGTACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCg                                                  <  2:4088/150‑1 (MQ=255)
                                                                                                     cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGt                                                <  1:1202513/150‑1 (MQ=255)
                                                                                                     cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGt                                                >  1:465633/1‑150 (MQ=255)
                                                                                                     cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCAGACACCCCGGt                                                >  1:1150614/1‑150 (MQ=255)
                                                                                                     cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTAATAGAAAGCGAACGCCGTGCTGTGGGCGGACAGCCCGGt                                                >  2:1840621/1‑150 (MQ=255)
                                                                                                     cccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTACATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGt                                                >  1:1930386/1‑150 (MQ=255)
                                                                                                      ccATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGtc                                               >  1:219031/1‑150 (MQ=255)
                                                                                                        aTAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGtctc                                             >  1:1544972/1‑150 (MQ=255)
                                                                                                          agagGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCt                                           <  2:465633/150‑1 (MQ=255)
                                                                                                               tCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTgcagg                                      <  1:2352511/150‑1 (MQ=255)
                                                                                                                  tAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGcc                                                    >  1:1343574/1‑133 (MQ=255)
                                                                                                                  tAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGcc                                                    <  2:1343574/133‑1 (MQ=255)
                                                                                                                    gCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCg                                 >  1:2399286/1‑150 (MQ=255)
                                                                                                                     ccGTAAACCCAGGCTGGCTCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCgg                                <  1:1426027/150‑1 (MQ=255)
                                                                                                                     ccGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTACATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCgg                                <  2:1930386/150‑1 (MQ=255)
                                                                                                                     ccGTAAACCCAGGCTGGATCTTCCTGGCTGAAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCgg                                <  2:2132867/150‑1 (MQ=255)
                                                                                                                       gTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGt                              <  2:2107723/150‑1 (MQ=255)
                                                                                                                         aaaCCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGACAGGCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCa                            >  1:1097001/1‑150 (MQ=255)
                                                                                                                             ccAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCAt                        >  2:1246206/1‑150 (MQ=255)
                                                                                                                              cAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCtt                                                                                        >  2:1126764/1‑85 (MQ=255)
                                                                                                                              cAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCtt                                                                                        <  1:1126764/85‑1 (MQ=255)
                                                                                                                              cAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATg                       >  1:1359408/1‑150 (MQ=255)
                                                                                                                               aGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATgg                      <  2:113987/150‑1 (MQ=255)
                                                                                                                                 gCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTAGGGATACTCGGGCTCCTGCAGGCAGCGGGTCAGCATGGcc                    >  2:788862/1‑150 (MQ=255)
                                                                                                                                  cgggATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCcg                   >  1:1403473/3‑150 (MQ=255)
                                                                                                                                  cTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCACTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGGCAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCcg                   >  1:1708467/1‑150 (MQ=255)
                                                                                                                                   tGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCg                                                           <  1:1218190/109‑1 (MQ=255)
                                                                                                                                         ttCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTc            >  2:942604/1‑150 (MQ=255)
                                                                                                                                         ttCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCGGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTc            >  2:1047074/1‑150 (MQ=255)
                                                                                                                                           ccTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAAc                                                                          <  1:2374838/86‑1 (MQ=255)
                                                                                                                                           ccTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAAc                                                                          >  2:2374838/1‑86 (MQ=255)
                                                                                                                                            cTGGCTCTACAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCgct         >  2:475861/1‑148 (MQ=255)
                                                                                                                                             tGGCTCTACAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCCGGCAGCGGGTCAGCATGGCCGCGTCCTCGATg        >  2:1072311/1‑150 (MQ=255)
                                                                                                                                              ggCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATgg       >  1:133146/1‑150 (MQ=255)
                                                                                                                                              ggCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATgg       >  1:132383/1‑150 (MQ=255)
                                                                                                                                               gCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGtc                                               >  1:410855/1‑109 (MQ=255)
                                                                                                                                               gCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGtc                                               <  2:410855/109‑1 (MQ=255)
                                                                                                                                                cTGTAGAGCCAGGTGTTGGCGTTGGATACAGATTGCACAAGTGATGTCCTCTCTTTTCTGTTGGCTTCATACACTAGTAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGcc     <  1:533758/150‑1 (MQ=255)
                                                                                                                                                  gTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCa                         <  1:1508116/128‑1 (MQ=255)
                                                                                                                                                  gTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCa                         >  2:1508116/1‑128 (MQ=255)
                                                                                                                                                  gTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGCCAt   <  2:571067/150‑1 (MQ=255)
                                                                                                                                                  gTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGCCAt   >  1:2170588/1‑150 (MQ=255)
                                                                                                                                                   tAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGCCATg  <  2:208928/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTTGCCGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTCGATGGCCATG  >  NC_002947/4446389‑4446685

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